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1 Introduction 



1.1 Overview 

This review will be concerned with our knowledge of extended matter under the 
governance of strong interaction, in short: QCD matter. Strictly speaking, the 
hadrons are representing the first layer of extended QCD architecture. In fact 
we encounter the characteristic phenomena of confinement as distances grow to 
the scale of 1 fm (i.e. hadron size): loss of the chiral symmetry property of 
the elementary QCD Lagrangian via non-perturbative generation of "massive" 
quark and gluon condensates, that replace the bare QCD vacuum [i\. However, 
given such first experiences of transition from short range perturbative QCD 
phenomena (jet physics etc.), toward extended, non perturbative QCD hadron 
structure, we shall proceed here to systems with dimensions far exceeding the 
force range: matter in the interior of heavy nuclei, or in neutron stars, and pri- 
mordial matter in the cosmological era from electro- weak decoupling (10~^^ s) 
to hadron formation (0.5 •10~^s). This primordial matter, prior to hadroniza- 
tion, should be deconfined in its QCD sector, forming a plasma (i.e. color 
conducting) state of quarks and gluons [2\: the Quark Gluon Plasma (QGP). 

In order to recreate matter at the corresponding high energy density in the ter- 
restial laboratory one collides heavy nuclei (also called "heavy ions") at ultrarel- 
ativistic energies. Quantum Chromodynamics predicts OElll] a phase transfor- 
mation to occur between deconfined quarks and confined hadrons. At near-zero 
net baryon density (corresponding to big bang conditions) non-perturbative 
Lattice-QCD places this transition at an energy density of about 1 GeV/fm^, 
and at a critical temperature, Tcrit ~ 170 MeV [H El El El E] (see the article 
on Lattice QCD in this Volume). The ultimate goal of the physics with ultra- 
relativistic heavy ions is to locate this transition, elaborate its properties, and 
gain insight into the detailed nature of the deconfined QGP phase that should 
exist above. What is meant by the term " ultrarelativistic" is defined by the 
requirement that the reaction dynamics reaches or exceeds the critical density 
e « 1 GeV/ fm? . Required beam energies turn out [8] to be y^ > 10 GeV , and 
various experimental programs have been carried out or are being prepared 
at the CERN SPS (up to about 20 GeV), at the BNL RHIC cohider (up to 
200 GeV) and finally reaching up to 5.5 TeV at the LHC of CERN. 

QCD confinement-deconfinement is of course not limited to the domain that 
is relevant to cosmological expansion dynamics, at very small excess of baryon 
over anti-baryon number density and, thus, near zero baryo-chemical potential 
^B- In fact, modern QCD suggests [9l[T0l[TT] a detailed phase diagram of QCD 
matter and its states, in the plane of T and baryo-chemical potential ^b- For a 
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map of the QCD matter phase diagram we are thus employing the terminology 
of the grand canonical Gibbs ensemble that describes an extended volume V of 
partonic or hadronic matter at temperature T. In it, total particle number is 
not conserved at relativistic energy, due to particle production-annihilation pro- 
cesses occurring at the microscopic level. However, the probability distributions 
(partition functions) describing the relative particle species abundances have to 
respect the presence of certain, to be conserved net quantum numbers (i), no- 
tably non-zero net baryon number and zero net strangeness and charm. Their 
global conservation is achieved by a thermodynamic trick, adding to the system 
Lagrangian a so-called Lagrange multiplier term, for each of such quantum num- 
ber conservation tasks. This procedure enters a "chemical potential" /ij that 
modifies the partition function via an extra term exp(— /ij/T) occuring in the 
phase space integral (see chapter [3] for detail). It modifies the canonical "pun- 
ishment factor" (exp (—E/T)), where E is the total particle energy in vacuum, 
to arrive at an analogous grand canonical factor for the extended medium, of 
exp {—E/T — Hi/T). This concept is of prime importance for a description of 
the state of matter created in heavy ion collisions, where net-baryon number 
(valence quarks) carrying objects are considered — extended "fireballs" of QCD 
matter. The same applies to the matter in the interior of neutron stars. The 
corresponding conservation of net baryon number is introduced into the grand 
canonical statistical model of QCD matter via the " baryo-chemical potential" 

We employ this terminology to draw a phase diagram of QCD matter in fig- 
ure ll.lt ill the variables T and fis ■ Note that fiB is high at low energies of 
collisions creating a matter fireball. In a head-on collision of two mass 200 
nuclei at y/s = 15GeV the fireball contains about equal numbers of newly cre- 
ated quark-antiquark pairs (of zero net baryon number), and of initial valence 
quarks. The accomodation of the latter, into created hadronic species, thus 
requires a formidable redistribution task of net baryon number, refiecting in a 
high value of fis- Conversely, at LHC energy (5.5TeV for Pb+Pb collisions), 
the initial valence quarks constitute a mere 5% fraction of the total quark den- 
sity, correspondingly requiring a small value of /x^. In the extreme, big bang 
matter evolves toward hadronization (at r=170 MeV) featuring a quark over 
antiquark density excess of 10^^ only, resulting in ^b ~ 0. 

Note that the limits of existence of the hadronic phase are not only reached by 
temperature increase, to the so-called Hagedorn value Th (which coincides with 
Tcrit at fiB -^ 0), but also by density increase to q > {5 — 10) qq: "cold com- 
pression" beyond the nuclear matter ground state baryon density ^o of about 
0.16 B/fm^. We are talking about the deep interior sections of neutron stars 
or about neutron star mergers |121 [T3| I14j . A sketch of the present view of 
the QCD phase diagram (Oj ^0\ [H] is given in Fig. 11.11 It is dominated by 
the parton-hadron phase transition line that interpolates smoothly between the 
extremes of predominant matter heating (high T, low fis) and predominant 
matter compression (T — > 0, hb > 1 GeV). Onward from the latter conditions, 
the transition is expected to be of first order }15J until the critical point of QCD 
matter is reached at 200 < fiB (E) < 500 AleV. The relatively large position 
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uncertainty reflects the preliminary character of Lattice QCD calculations at 
finite ^b [9l[T0l[Tl]. Onward from the critical point, E, the phase transforma- 
tion at lower fis is a cross-over [11 j. 
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Figure 1.1: Sketch of the QCD matter phase diagram in the plane of tempera- 
ture T and baryo-chemical potential ij,b- The parton-hadron phase 
transition line from lattice QCD [HI [9l 1101 111) ends in a critical point 
E. A cross-over transition occurs at smaller /x^. Also shown are 
the points of hadro-chemical freeze-out from the grand canonical 
statistical model. 



We note, however, that these estimates represent a major recent advance of lat- 
tice theory which was, for two decades, believed to be restricted to the fj-B = 
situation. Onward from the critical point, toward lower ^ub, the phase transfor- 
mation should acquire the properties of a rapid cross-over [TC], thus also includ- 
ing the case of primordial cosmological expansion. This would finally rule out 
former ideas, based on the picture of a violent first order "explosive" cosmolog- 
ical hadronization phase transition, that might have caused non-homogeneous 
conditions, prevailing during early nucleo-synthesis [T7], and fluctuations of 
global matter distribution density that could have served as seedlings of galac- 
tic cluster formation [TH]. However, it needs to be stressed that the conjectured 
order of phase transformation, occuring along the parton - hadron phase bound- 
ary line, has not been unambiguously confirmed by experiment, as of now. 



On the other hand, the position of the QCD phase boundary at low fiB has, in 
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fact, been located by the hadronization points in the T, fiB plane that are also 
illustrated in Fig. II. li They are obtained from statistical model analysis [19] 
of the various hadron multiplicities created in nucleus- nucleus collisions, which 
results in a [T, /i^] determination at each incident energy, which ranges from 
SIS via AGS and SPS to RHIC energies, i.e. 3 < ^ < 200 GeV. Toward low 
Hb these hadronic freeze-out points merge with the lattice QCD parton-hadron 
coexistence line: hadron formation coincides with hadronic species freeze-out. 
These points also indicate the ^b domain of the phase diagram which is acces- 
sible to relativistic nuclear collisions. The domain at fiB > 1-5 GeV which is 
predicted to be in a further new phase of QCD featuring color-flavor locking 
and color superconductivity [20] will probably be accessible only to astrophys- 
ical observation. 

One may wonder how states and phases of matter in thermodynamical equi- 
librium - as implied by a description in grand canonical variables - can be 
sampled via the dynamical evolution of relativistic nuclear collisions. Employ- 
ing heavy nuclei, A ss 200, as projectiles/targets or in colliding beams (RHIC, 
LHC), transverse dimensions of the primordial interaction volume do not ex- 
ceed about 8 fm, and strong interaction ceases after about 20 fm/c. We shall 
devote an entire later section to the aspects of equilibrium (12. 5p but note, for 
now, that the time and dimension scale of primordial perturbative QCD inter- 
action at the microscopic partonic level amounts to subfractions of 1 fm/c, the 
latter scale, however, being representative of non perturbative processes (con- 
finement, "string" formation etc.). The A+A fireball size thus exceeds, by far, 
the elementary non perturbative scale. An equilibrium quark gluon plasma rep- 
resents an extended non-perturbative QCD object, and the question whether 
its relaxation time scale can be provided by the expansion time scale of an A+A 
collision, needs careful examination. Reassuringly, however, the hadrons that 
are supposedly created from such a preceding non-perturbative QGP phase at 
top SPS and RHIC energy, do in fact exhibit perfect hadrochemical equilib- 
rium, the derived [T, hb] values p^ thus legitimately appearing in the phase 
diagram. Fig. 11.11 

In the present review we will order the physics observables to be treated, in 
sequence of their origin from successive stages that characterize the overall dy- 
namical evolution of a relativistic nucleus-nucleus collision. In rough outline 
this evolution can be seen to proceed in three major steps. An initial period of 
matter compression and heating occurs in the course of interpenetration of the 
projectile and target baryon density distributions. Inelastic processes occuring 
at the microscopic level convert initial beam longitudinal energy to new internal 
and transverse degrees of freedom, by breaking up the initial baryon structure 
functions. Their partons thus acquire virtual mass, populating transverse phase 
space in the course of inelastic perturbative QCD shower multiplication. This 
stage should be far from thermal equilibrium, initially. However, in step two, 
inelastic interaction between the two arising parton fields (opposing each other 
in longitudinal phase space) should lead to a pile-up of partonic energy density 
centered at mid-rapidity (the longitudinal coordinate of the overall center of 
mass). Due to this mutual stopping down of the initial target and projectile 
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parton fragmentation showers, and from the concurrent decrease of parton vir- 
tuahty (with decreasing average square momentum transfer Q^) there results 
a slowdown of the time scales governing the dynamical evolution. Equilibrium 
could be approached here, the system "lands" on the T, fj, plane of Fig. II. 1^ 
at temperatures of about 300 and 200 MeV at top RHIC and top SPS energy, 
respectively. The third step, system expansion and decay, thus occurs from 
well above the QCD parton-hadron boundary line. Hadrons and hadronic reso- 
nances then form, which decouple swiftly from further inelastic transmutation 
so that their yield ratios become stationary ("frozen-out"). A final expansion 
period dilutes the system to a degree such that strong interaction ceases all 
together. 

In order to verify in detail this qualitative overall model, and to ascertain the 
existence (and to study the properties) of the different states of QCD that 
are populated in sequence, one seeks observable physics quantities that convey 
information imprinted during distinct stages of the dynamical evolution, and 
"freezing-out" without significant obliteration by subsequent stages. Ordered 
in sequence of their formation in the course of the dynamics, the most relevant 
such observables are briefly characterized below: 

1. Suppression of J/^ and Y production by Debye-screening in the QGP. 
These vector mesons result from primordial, pQCD production of cc and 
bb pairs that would hadronize unimpeded in elementary collisions but 
are broken up if immersed into a npQCD deconfined QGP, at certain 
characteristic temperature thresholds. 

2. Suppression of dijets which arise from primordial qq pair production frag- 
menting into partonic showers (jets) in vacuum but being attenuated by 
QGP-medium induced gluonic bremsstrahlung: Jet quenching in A+A 
collisions. 

a) A variant of this: any primordial hard parton suffers a high, spe- 
cific loss of energy when traversing a deconfined medium: High pT 
suppression in A+A collisions. 

3. Hydrodynamic collective motion develops with the onset of (local) thermal 
equilibrium. It is created by partonic pressure gradients that reflect the 
initial collisional impact geometry via non-isotropies in particle emission 
called "directed" and "elliptic" flow. The latter reveals properties of the 
QGP, seen here as an ideal partonic fluid. 

a) Radial hydrodynamical expansion flow ("Hubble expansion") is a 
variant of the above that occurs in central, head on collisions with 
cylinder symmetry, as a consequence of an isentropic expansion. It 
should be sensitive to the mixed phase conditions characteristic of a 
first order parton-hadron phase transition. 

4. Hadronic "chemical" freeze-out fixes the abundance ratios of the hadronic 
species into an equilibrium distribution. Occuring very close to, or at 
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hadronization, it reveals the dynamical evolution path in the [T, /u^] plane 
and determines the critical temperature and density of QCD. The yield 
distributions in A+A collisions show a dramatic strangeness enhancement 
effect, characteristic of an extended QCD medium. 

5. Fluctuations, from one collision event to another (and even within a single 
given event) can be quantified in A+A collisions due to the high charged 
hadron multiplicity density (of up to 600 per rapidity unit at top RHIC 
energy) . Such event-by-event (ebye) fluctuations of pion rapidity density 
and mean transverse momentum (event "temperature"), as well as event- 
wise fluctuations of the strange to non-strange hadron abundance ratio 
(may) reflect the existence and position of the conjectured critical point 
of QCD (Fig.in]). 

6. Two particle Bose-Einstein-Correlations are the analog of the Hanbury- 
Brown, Twiss (HBT) effect of quantum optics. They result from the last 
interaction experienced by hadrons, i.e. from the global decoupling stage. 
Owing to a near isentropic hadronic expansion they reveal information on 
the overall space-time-development of the "fireball" evolution. 

In an overall view the first group of observables (1 to 2a) is anchored in es- 
tablished pQCD physics that is well known from theoretical and experimental 
analysis of elementary collisions (e^e~ annihilation, pp and pp data). In fact, 
the first generation of high Q^ baryon collisions, occuring at the microscopic 
level in A+A collisions, should closely resemble such processes. However, their 
primary partonic products do not escape into pQCD vacuum but get atten- 
uated by interaction with the concurrently developing extended high density 
medium , thus serving as diagnostic tracer probes of that state. The remaining 
observables capture snapshots of the bulk matter medium itself. After initial 
equilibration we may confront elliptic flow data with QCD during the corre- 
sponding partonic phase of the dynamical evolution employing thermodynamic 
|21j and hydrodynamic [22] models of a high temperature parton plasma. The 
hydro-model stays applicable well into the hadronic phase. Hadron formation 
(confinement) occurs in between these phases (at about 5 microseconds time in 
the cosmological evolution). In fact relativistic nuclear collision data may help 
to finally pin down the mechanism(s) of this fascinating QCD process |23ll24t [25] 
as we can vary the conditions of its occurence, along the parton-hadron phase 
separation line of Fig. II. H by proper choice of collisional energy ^/s, and sys- 
tem size A, while maintaining the overall conditions of an extended imbedding 
medium of high energy density within which various patterns [HI \T0\ \TT\ [T5| [TB] 
of the hadronization phase transition may establish. The remaining physics ob- 
servables (3a, 5 and 6 above) essentially provide for auxiliary information about 
the bulk matter system as it traverses (and emerges from) the hadronization 
stage, with special emphasis placed on manifestations of the conjectured critical 
point. 

The present review will briefly cover each of the above physics observables in a 
separate chapter, beginning with the phenomena of confinement and hadroniza- 
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tion (chapter [3]) , then to turn to the preceding primordial dynamics, e.g. to 
elhptical flow (chapter H]), high pT and jet quenching (chapter [5|) and quarko- 
nium suppression (chapter [6]) as weh as in- medium £»-meson "melting". We then 
turn to the late period, with correlation and fluctuation studies (chapter [7]) . We 
conclude (chapter [S]) with a summary, including an outlook to the future of the 
research field. 

However, before turning to such specific observables we shall continue this intro- 
ductory chapter, with a look at the origin, and earlier development of the ideas 
that have shaped this field of research (section [L2]) . Then we turn to a detailed 
description of the overall dynamical evolution of relativistic nucleus- nucleus col- 
lisions, and to the typical overall patterns governing the final distributions in 
transverse and longitudinal (rapidity) phase space (chapter [2]). The aspects 
of an approach toward equilibrium, at various stages of the dynamical evolu- 
tion (which are of key importance toward the intended elucidation of the QCD 
matter phase diagram), will be considered, in particular. 



1.2 History 

The search for the phase diagram of strongly interacting matter arose in the 
1960s, from a coincidence of ideas developing - at first fairly independently - in 
nuclear and astrophysics. In fact, the nuclear proton-neutron matter, a quan- 
tum liquid at T = and energy density e = 0.15 GeV/fm^, represents the 
ground state of extended QCD matter. Of course, QCD was unknown dur- 
ing the development of traditional nuclear physics, and the extended matter 
aspects of nuclei - such as compressibility or the equation of state, in general 
- did not receive much attention until the advent, in the 1960s, of relativistic 
nuclear mean field theory, notably s-matrix theory by Brueckner [2S] and the 
(T-model of Walecka |27) . These theories developed the novel view of "infinite 
nuclear matter" structure, based on in-medium properties of the constituent 
baryons that share parts of their vacuum mass and surface structure with the 
surrounding, continuous field of relativistic scalar and vector mesons. Most 
importantly, in the light of subsequent development, these theories allowed for 
a generalization away from ground state density and zero temperature. Such 
developments turned out to be of key relevance for acute nuclear astrophysics 
problems: the dynamics of type II super-novae and the stability of neutron 
stars, which both required the relation of pressure to density and temperature 
of hadronic matter, i.e. the hadronic matter equation of state (EOS). H. A. 
Bethe, G. Brown and their collaborators f7E\ postulated that the final stages 
of supernova collaps should evolve through the density interval 0.1 < g/go < 5 
where go = 0.16 (baryons per /m^) is the nuclear matter ground state den- 
sity, and a similar domain was expected for the neutron star density variation 
from surface to interior [29]. It was clear that, at the highest thus considered 
densities the EOS might soften due to strange hadron production caused by 
increasing Fermi energy. However the field theoretical models permitted no 
reliable extrapolation to such high densities (which, in retrospect, are perhaps 
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not reached in supernova dynamics [30]), and the experimental information 
concerning the EOS from study of the giant monopole resonance - a collective 
density oscillation also called "breathing mode" - covered only the parabolic 
minimum occuring in the related function of energy vs. density at T = 0, 
= 00- 

The situation changed around 1970 due to the prediction made by W. Greiner 
and his collaborators [3T] that nucleus-nucleus collisions, at relatively mod- 
est relativistic energies, would result in shock compression. This mechanism 
promised to reach matter densities far beyond those of a mere superposition 
(i.e. q/qo < 27) of initial target and projectile densities. Coinciding in time, 
the newly developed Bevalac accelerator at LBL Berkeley offered projectiles up 
to ^^Ar, at just the required energies, WOMeV < Eiab/nudeon < 2GeV. The 
field of "relativistic heavy ion physics" was born. The topic was confronted, 
at first, with experimental methods available both from nuclear and particle 
physics. It was shown that particle production (here still restricted to pions and 
kaons) could indeed be linked to the equation of state [32] and that, even more 
spectacularly, the entire "fireball" of mutually stopped hadrons developed decay 
modes very closely resembling the initial predictions of hydro dynamical shock 
flow modes |33j which directly link primordial pressure gradients with collective 
velocity fields of matter streaming out, again governed by the nuclear/hadronic 
matter EOS. Actually, both these statements do, in fact, apply (mutatis mutan- 
dis) up to the present ultra-relativistic energies (see chapters [21 [3l |1|) . However 
it turned out soon that the equation of state at low or even zero temperature 
(as required in supernova and neutron star studies) could only be obtained in 
a semi-empirical manner [3^. The reason: compression can, in such collisions, 
be only accomplished along with temperature and entropy increase. In an ideal 
baryon gas exp{s/A) oc T^''^/q, i.e. T^/"^ will grow faster than g in a non- 
isentropic compression. Thus the reaction dynamics will be sensitive to various 
isothermes of the ground state EOS P = f{g, T = 0) , staying at T » 0, 
throughout, and, moreover, not at constant T. Thus a relativistic dynamical 
mean field model is required in order to interactively deduce the T = EOS 
from data [31]. The EOS result thus remains model dependent. 



The ideas concerning creation of a quark gluon plasma arose almost concurrent 
with the heavy ion shock compression proposal. In 1974 T. D. Lee formulated 
the idea that the non-perturbative vacuum condensates could be "melted down 
. . by distributing high energy or high nucleon density over a relatively large 
volume" [35]. Collins and Perry [36] realized that the asymptotic freedom 
property of QCD implies the existence of an ultra-hot form of matter with 
deconfined quarks and gluons, an idea that gained wide recognition when S. 
Weinberg |37J proposed an asymptotic freedom phase at the beginning of "The 
first Three minutes" . In fact, this idea of deconfinement by asymptotic freedom 
(with implied temperature of several GeV) was correct, but somewhat besides 
the point, as everybody expected, likewise, that deconfinement sets in right 
above the limiting hadron temperature of R. Hagedorn [38], Th ~ 160 MeV. A 
medium existing down to that temperature would, however, feature an average 
momentum square transfer Q^ < 1 GeV^, i.e. be far into the non perturbative 
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domain, and very far from asymptotic freedom. Right above the hadron to 
parton transition the "quark gluon plasma" (as it was named by E. Shuryak [39]) 
is not a weakly coupled ideal pQCD gas as soon became obvious by Lattice QCD 
calculations for extended matter [40j. Seen in retrospect one obviously can not 
defend a picture of point like quarks (with "current" masses) at Q^ < 0.2GeV^ 
where size scales of 0.5 to 1 fm must play a dominating role. 

An analytic QCD description of deconfinement does not exist. For heavy 
quarkonia, cc (J/^) and bb (Y) deconfinement in partonic matter, Matsui and 
Satz proposed [51] a Debye screening mechanism, caused by the high spatial 
density of free color carriers, that removes the confining long range potential 
as T increases toward about 2 Tc, an effect reproduced by modern lattice QCD 
[42j . However, light hadron deconfinement can not be understood with a non- 
relativistic potential model. Such critical remarks not withstanding, we shall 
demonstrate in chapters [3] to [H] that the very existence, and also crucial prop- 
erties of the QGP can in fact be inferred from experiment, and be confronted 
with corresponding predictions of recent lattice QCD theory. 

Our present level of an initial understanding of the phase diagram of QCD 
matter (Fig. II. ip , is the result of a steady development of both experiment and 
theory, that began about three decades ago, deriving initial momentum from the 
Bevalac physics at LBL which motivated - along with the developing formulation 
of the quark gluon plasma research goals - a succession of experimental facilities 
progressing toward higher ^/s. Beginning with the AGS at BNL (^^Si and ^^'^Au 
beams with ^/s < 5 GeV), the next steps were taken at the CERN SPS {^/s 
from 6 to 20 GeV; '^^O, ^^5, "^^^Pb beams), and at the Relativistic Heavy Ion 
Collider RHIC (the first facility constructed explicitly for nuclear collisions) 
which offers beams of ^^Cu and ^^"^Au at 20 < ^/s < 200 GeV. A final, gigantic 
step in energy will be taken 2008 with the CERN Large Hadron Collider: '^^^Pb 
beams at y/s = 5.5 TeV. 
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2 Bulk Hadron Production in A+A 
Collisions 



In this section we take an overall look at bulk hadron production in nucleus- 
nucleus collisions. In view of the high total cm. energies involved at e.g. top 
SPS (£;*^ w 3.3 TeV) and top RHIC (38 TeV) energies, in central Pb+Pb 
(SPS) and Au+Au (RHIC) collisions, one can expect an extraordinarily high 
spatial density of produced particles. Thus, as an overall idea of analysis, one 
will try to relate the observed flow of energy into transverse and longitudinal 
phase space and particle species to the high energy density contained in the 
primordial interaction volume, thus to infer about its contained matter. The 
typical experimental patterns of such collisions, both in collider mode at RHIC 
and in a fixed target configuration at the SPS, are illustrated in Fig. 12. II which 
shows a fractional view of the total distribution of charged particles (about 4000 
and 1600, respectively) within the tracking volume of the STAR and NA49 ex- 
periments. 

Most of these tracks correspond to "thermal" pions {pT up to 2 GeV) and, in 
general, such thermal hadrons make up for about 95% of the observed multi- 
plicity: the bulk of hadron production. Their distributions in phase space will 
be illustrated in the subsections below. This will lead to a first insight into 
the overall reaction dynamics, and also set the stage for consideration of the 
rare signals, imbedded in this thermal bulk production: correlations, jets, heavy 
flavors, fluctuations, which are the subject of later chapters. 



2.1 Particle Multiplicity and Transverse Energy Density 

Particle production can be assessed globally by the total created transverse en- 
ergy, the overall result of the collisional creation of transverse momentum px or 



transverse mass {mx = \ Pr + ^)) at the microscopic level. Fig. l2.2l shows the 
distribution of total transverse energy Et = X] -^(©i) ' sinQ resulting from a 

i 

calorimetric measurement of energy flow into calorimeter cells centered at angle 
0j relative to the beam [13], for '^^5+^^'^Au collisions at y/s = 20 GeV , and for 
208p6+208p]3 collisions at Vs = 17.3 GeV. 

The shape is characteristic of the impact parameter probability distribution 
(for equal size spheres in the Pb+Pb case). The turnoff at Et = 520 GeV 
indicates the point where geometry runs out of steam, i.e. where 6 — > 0, a 
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2 Bulk Hadron Production in A+A Collisions 




Figure 2.1: Charged particle tracks in central Au+Au and Pb+Pb collision 
events, in collider geometry (top) from RHIC STAR TPC track- 
ing at y^ = 200 GeV, and in fixed target geometry (bottom) from 
NA49 at the SPS, ^/s = 17.3 GeV. 
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2.1 Particle Multiplicity and Transverse Energy Density 
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Figure 2.2: Minimum bias distribution of total transverse energy in Pb+Pb 
collisions at ^/s = 17.3 GeV, and S+Au collisions at ^/s = 20 GeV, 
in the rapidity interval 2.1 < y < 3.4, from [43J. 



configuration generally referred to as a "central collision". The adjacent shoul- 
der results from genuine event by event fluctuations of the actual number of 
participant nucleons from target and projectile (recall the diffuse Woods-Saxon 
nuclear density profiles), and from experimental factors like calorimeter reso- 
lution and limited acceptance. The latter covers 1.3 units of pseudo-rapidity 
and contains mid-rapidity ijmid = 2.9. Re-normalizing [43j to A r/ = 1 leads 
to dEx / dr][mid) = 400 GeV , in agreement with the corresponding WA80 re- 
sult |44| . Also, the total transverse energy of central Pb-|-Pb collisions at 
y/s = 17.3 GeV turns out to be about 1.2 TeV . As the definition of a cen- 
tral collision, indicated in Fig. 12.21 can be shown [42j to correspond to an 
average nucleon participant number of Npart = 370 one finds an average total 
transverse energy per nucleon pair, of Et/ (0.5 Npart) = 6.5 GeV . After proper 
consideration of the baryon pair rest mass (not contained in the calorimetric 
Et response but in the corresponding ^/s) one concludes [43j that the observed 
total Et corresponds to about 0.6 EJp""^ , the maximal Et derived from a situa- 
tion of " complete stopping" in which the incident y/s gets fully transformed into 
internal excitation of a single, ideal isotropic fireball located at mid-rapidity. 
The remaining fraction of EJp'^^ thus stays in longitudinal motion, reflecting 
the onset, at SPS energy, of a transition from a central flreball to a longitudi- 
nally extended "fire-tube", i.e. a cylindrical volume of high primordial energy 
density. In the limit of much higher y^ one may extrapolate to the idealization 
of a boost invariant primordial interaction volume, introduced by Bjorken |j45j. 

We shall show below (section 12. 2p that the charged particle rapidity distribu- 
tions, from top SPS to top RHIC energies, do in fact substantiate a development 
toward a boost-invariant situation. One may thus employ the Bjorken model 
for an estimate of the primordial spatial energy density e, related to the energy 
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density in rapidity space via the relation [l5 



<ro) 



1 1 dEr 
ttR'^ tq dy 



(2.1) 



where the initially produced collision volume is considered as a cylinder of length 
dz = Tody and transverse radius R oc A^'^. Inserting for vri?^ the longitudinally 
projected overlap area of Pb nuclei colliding near head-on ("centrally"), and 
assuming that the evolution of primordial pQCD shower multiplication (i.e. the 
energy transformation into internal degrees of freedom) proceeds at a time scale 
To < Ifm/c, the above average transverse energy density, of dEx/dy = AOOGeV 
at top SPS energy [iSl E] leads to the estimate 



e(ro = Ifm) = 3.0 ± 0.6 GeV/fm^ 



(2.2) 



thus exceeding, by far, the estimate of the critical energy density eo obtained 
from lattice QCD (see below), of about 1.0 GeV/ fm^. Increasing the collision 
energy to ^/s = 200 GeV for Au+Au at RHIC, and keeping the same formation 
time. To = 1 fm/c (a conservative estimate as we shall show in section [23]) , the 
Bjorken estimate grows to e ~ 6.0 it 1 GeV/ fm? . This statement is based on 
the increase of charged particle multiplicity density at mid-rapidity with -y/s, 
as illustrated in Fig. 12.31 From top SPS to top RHIC energy \iQ\ the density 
per participant nucleon pair almost doubles. However, at ^/s = 200 GeV the 
formation or thermalization time tq, employed in the Bjorken model [l5], was 
argued [17] to be shorter by a factor of about 4. We will return to such estimates 
of To in section 12.51 but note, for now, that the above choice of tq = 1 fm/c 
represents a conservative upper limit at RHIC energy. 
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Figure 2.3: Charged hadron rapidity density at mid-rapidity vs. ^/s, compiled 
from e^e~, pp, pp and A+A collisions 
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2.1 Particle Multiplicity and Transverse Energy Density 

These Bjorken-estimates of spatial transverse energy density are confronted in 
Fig. 12.41 with lattice QCD results obtained for three dynamical light quark 
flavors 05, and for zero baryo-chemical potential (as is realistic for RHIC 
energy and beyond but still remains a fair approximation at top SPS energy 
where hb ~ 250 MeV). The energy density of an ideal, relativistic parton gas 
scales with the fourth power of the temperature, 

e = gT^ (2.3) 

where g is related to the number of degrees of freedom. For an ideal gluon 
gas, g = 16 7r^/30; in an interacting system the effective g is smaller. The 
results of Fig. 12.41 show, in fact, that the Stefan-Boltzmann limit esB is not 
reached, due to non perturbative effects, even at four times the critical tem- 
perature Tc = 170 MeV. The density e/T^ = g is seen to ascend steeply, 
within the interval Tc =b 25 MeV. At Tc the critical QCD energy density 
e = 0.6 — 1.0 GeV/fm^. Relating the thermal energy density with the Bjorken 
estimates discussed above, one arrives at an estimate of the initial temperatures 
reached in nucleus-nucleus collisions, thus implying thermal partonic equilib- 
rium to be accomplished at time scale tq (see section [23]) . For the SPS, RHIC 
and LHC energy domains this gives an initial temperature in the range 190 
< T^P^ < 220 MeV, 220 < T^^^"^ < 400 MeV (assuming [37] that tq de- 
creases to about 0.3 fm/c here) and T^^^ > 600 MeV, respectively. From 
such estimates one tends to conclude that the immediate vicinity of the phase 
transformation is sampled at SPS energy, whereas the dynamical evolution at 
RHIC and LHC energies dives deeply into the " quark-gluon-plasma" domain 
of QCD. We shall return to a more critical discussion of such ascertations in 
section I 



One further aspect of the mid-rapidity charged particle densities per participant 
pair requires attention: the comparison with data from elementary collisions. 
Fig. 12.31 shows a compilation of pp, pp and e~^e~' data covering the range from 
ISR to LEP and Tevatron energies. 

The data from e~^e~ represent dNch/dy, the rapidity density along the event 
thrust axis, calculated assuming the pion mass [39] (the difference between 
dN/dy and dN/dr] can be ignored here). Remarkably, they superimpose with 
the central A+A collision data, whereas pp and pp show similar slope but 
amount to only about 60% of the AA and e~^e~ values. This difference between 
e~^e~ annihilation to hadrons, and pp or pp hadro-production has been ascribed 
[50| to the characteristic leading particle effect of minimum bias hadron-hadron 
collisions which is absent in e^e~. It thus appears to be reduced in AA col- 
lisions due to subsequent interaction of the leading parton with the oncoming 
thickness of the remaining target/projectile density distribution. This naturally 
leads to the scaling of total particle production with Npart that is illustrated 
in Fig. 12.51 for three RHIC energies and minimum bias Au+Au collisions; the 
close agreement with e~^e~ annihilation data is obvious again. One might con- 
clude that, analogously, the participating nucleons get "annihilated" at high 
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Figure 2.4: Lattice QCD results at zero baryon potential for energy density 
e/T^ versus T/Tc with three light quark flavors, compared to the 
Stefan-Boltzmann- limit esB of an ideal quark-gluon gas [H] . 
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Figure 2.5: The total number of charged hadrons per participant pair shown as 
a function of Npart in Au+ Au collisions at three RHIC energies [53] . 



^i, their net quantum number content being spread out over phase space (as 
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we shall show in the next section). 



2.2 Rapidity Distributions 
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Figure 2.6: Left panel: charged particle pseudo-rapidity distribution in pp col- 
lisions at \fs = 540 GeV [51] • Right panel: same in RHIC Au+Au 
collisions at ^/s = 130 GeV at different centralities [SI]- Closed lines 
represent fits with the color glass condensate model 



2.2 Rapidity Distributions 



Particle production number in A+A collisions depends globally on ^/s and colli- 
sion centrality, and differentially on px and rapidity y, for each particle species 
i. Integrating over -px results in the rapidity distribution dNi/dy. Particle 

rapidit}o, y = sinh~^ pi/Mx (where Mt = \ nrfl +p^), requires mass identi- 
fication. If that is unknown one employs pseudo-rapidity (77 = —In [tan(G/2)]) 
instead. This is also chosen if the joint rapidity distribution of several unre- 
solved particle species is considered: notably the charged hadron distribution. 
We show two examples in Fig. 12. 6[ The left panel illustrates charged particle 
production in pp collisions studied by UAl at y/s = 540 GeV [^. Whereas 
the minimum bias distribution (dots) exhibits the required symmetry about 
the center of mass coordinate, ?y = 0, the rapidity distribution corresponding to 
events in which a W boson was produced (histogram) features, both, a higher 
average charged particle yield, and an asymmetric shape. The former effect can 



^The rapidity variable represents a compact (logarithmic) description of longitudinal phase 
space. It is based on longitudinal particle velocity (derived from piong and m), y = 
1/2 In ((1 + f3L)/{l — Pl))- The rapidity distribution dN/dy is shape invariant under lon- 
gitudinal Lorentz transformation, and centered at "mid-rapidity" ymid = ycM, for all 
produced particle species; see figures [231 to [2^ 
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be seen to reflect the expectation that the W production rate increases with 
the "centrahty" of pp cohisions, involving more primordial partons as the col- 
lisional overlap of the partonic density profiles gets larger, thus also increasing 
the overall, softer hadro-production rate. The asymmetry should result from a 
detector bias favoring W identification at negative rapidity: the transverse W 
energy, of about 100 GeV would locally deplete the energy store available for 
associated soft production. If correct, this interpretation suggests that the wide 
rapidity gap between target and projectile, arising at such high ^/s, of width 
Ay Ki 2 In {2jcm), makes it possible to define local sub-intervals of rapidity 
within which the species composition of produced particles varies. 

The right panel of Fig. [2T6l shows charged particle pseudo-rapidity density distri- 
butions for Au-|-Au collisions at ^/s = 130 GeV measured by RHIC experiment 
PHOBOS [52] at three different collision centralities, from "central" (the 6% 
highest charged particle multiplicity events) to semi-peripheral (the correspond- 
ing 35-45% cut). We will turn to centrality selection in more detail below. Let 
us first remark that the slight dip at mid-rapidity and, moreover, the distribu- 
tion shape in general, are common to pp and Au-|-Au. This is also the case for 
e+e~ annihilation as is shown in Fig. 12.71 which compares the ALEPH rapid- 
ity distribution along the mean pj- ("thrust") axis of jet production in e'^e~ 
at ^ = 200 GeV 09] with the scaled PHOBOS-RHIC distribution of central 
Au-|-Au at the same ^/s |53] . Note that the mid-rapidity values contained in 
Figs. 12.61 and 12.71 have been employed already in Fig. 12.31 which showed the 
overall a/s dependence of mid-rapidity charged particle production. What we 
concluded there was a perfect scaling of A-l-A with e~^e~ data at ^/s > 20 GeV 
and a 40% suppression of the corresponding pp, pp yields. We see here that 
this observation holds, semi-quantitatively, for the entire rapidity distributions. 
These are not ideally boost invariant at the energies considered here but one sees 
in dNch/dy a relatively smooth "plateau" region extending over | y |< 1.5 — 2.5. 

The production spectrum of charged hadrons is, by far, dominated by soft 
pions {pt < 1 GeV/c) which contribute about 85% of the total yield, both in 
elementary and nuclear collisions. The evolution of the vr^ rapidity distribution 
with ^/s is illustrated in Fig. 12.81 for central Au-|-Au and Pb-|-Pb collisions from 
AGS via SPS to RHIC energy, 2.7 < ./s < 200 GeV [54J. 



At lower y^ the distributions are well described by single Gaussian fits 
with cr{y) nearly linearly proportional to the total rapidity gap Ay oc luy/s as 
shown in the right hand panel of Fig. 12. 8[ Also illustrated is the prediction of 
the schematic hydrodynamical model proposed by Landau [55] , 



,2 ^.i^rVs 



a^(xln{-^) (2.4) 



'2mp 



which pictures hadron production in high ^/s pp collisions to proceed via a 
dynamics of initial complete "stopping down" of the reactants matter/energy 
content in a mid-rapidity fireball that would then expand via 1-dimensional 
ideal hydrodynamics. Remarkably, this model that has always been considered 
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2.2 Rapidity Distributions 




Figure 2.7: Pseudo-rapidity distribution of charged hadrons produced in central 
Au+Au collisions at ^/s = 200 GeV compared with e~^e~ data at 
similar energy. The former data normalized by Npart/'^- From ref. 
153]. 



a wildly extremal proposal falls rather close to the lower ^/s data for central 
A+A collisions but, as longitudinal phase space widens approaching boost in- 
variance we expect that the (non-Gaussian) width of the rapidity distribution 
grows linearly with the rapidity gap Ay. LHC data will finally confirm this 
expectation, but Figs. 12.61 to 12.81 clearly show the advent of boost invariance, 
already at ^s = 200 GeV. 

A short didactic aside: At low y/s the total rapidity gap Ay = 2 — 3 does 
closely resemble the total rapidity width obtained for a thermal pion velocity 
distribution at temperature T = 120— IbQMeV , of a single mid-rapidity fireball, 
the y-distribution of which represents the longitudinal component according to 
the relation HH! 



dN , 2^ 2mr2 

— — oi[m T ^ -— 

ay cosny 



2T^ 



cos 



K^y 



exp [—mcoshy/T] 



(2.5) 



where m is the pion mass. Any model of preferentially longitudinal expan- 
sion of the pion emitting source, away from a trivial single central "completely 
stopped" fireball, can be significantly tested only once Ay > 3 which occurs 
upward from SPS energy. The agreement of the Landau model prediction with 
the data in Fig. 12.81 is thus fortuitous, below y/s ~ 10 GeV, as any created 
fireball occupies the entire rapidity gap with pions. 



The Landau model offers an extreme view of the mechanism of " stopping" , by 
which the initial longitudinal energy of the projectile partons or nucleons is 
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Figure 2.8: Left panel: negative pion rapidity distributions in central Au+Au 
and Pb+Pb collisions from AGS via SPS to RHIC energies |54j . 
Right panel: the Gaussian rapidity width of pions versus y^, con- 
fronted by Landau model predictions (solid line) [54J . 



inelastically transferred to produced particles and redistributed in transverse 
and longitudinal phase space, of which we saw the total transverse fraction in 
Fig. 12.21 Obviously €^e~ annihilation to hadrons represents the extreme stop- 
ping situation. Hadronic and nuclear collisions offer the possibility to analyze 
the final distribution in phase space of their non-zero net quantum numbers, 
notably net baryon number. Figure 12.91 shows the net-proton rapidity distri- 
bution (i.e. the proton rapidity distribution subtracted by the antiproton dis- 
tribution) for central Pb-hPb/Au+Au colhsions at AGS (^ = 5.5 GeV), SPS 
(^/s < 17.3 GeV) and RHIC (^ = 200 GeV) [56j. With increasing energy we 
see a central (but non-Gaussian) peak developing into a double-hump structure 
that widens toward RHIC leaving a plateau about mid-rapidity. The RHIC- 
BRAHMS experiment acceptance for p, p identification does unfortunately not 
reach up to the beam fragmentation domain at y^ = 5.4 (nor does any other 
RHIC experiment) but only to y ~ 3.2, with the consequence that the major 
fraction of p"^* is not accounted for. However the mid-rapidity region is by no 
means net baryon free. At SPS energy the NA49 acceptance covers the ma- 
jor part of the total rapidity gap, and we observe in detail a net p distribution 
shifted down from y.p = 2.9 by an average rapidity shift [SB] of {5y) = 1.7. From 
Fig. 12.91 we infer that {5y) can not scale linearly with yp ss ln{2^cM) ~ ln^/s 
for ever - as it does up to top SPS energy where {6y) = 0.58 yp [56j. Because 
extrapolating this relation to ^/s = 200 GeV would result in {6y) = 3.1, and 
with yp ~ 5.4 at this energy we would expect to observe a major fraction of net 
proton yield in the vicinity of y = 2.3 which is not the case. A saturation must 
thus occur in the {5y) vs. ^/s dependence. 

The re-distribution of net baryon density over longitudinal phase space is, of 
course, only partially captured by the net proton yield but a recent study 
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Figure 2.9: Net proton rapidity distributions in central Au+Au/Pb+Pb colli- 
sions at AGS, SPS and RHIC energies [MIEZ]. 



has shown that proper inclusion of neutron □ and hyperon production data at 
SPS and RHIC energy scales up, of course, the dN/dy distributions of Fig. 12.91 
but leaves the peculiarities of their shapes essentially unchanged. As the net 
baryon rapidity density distribution should resemble the final valence quark 
distribution the Landau model is ruled out as the valence quarks are seen to be 
streaming from their initial position at beam rapidity toward mid-rapidity (not 
vice versa). It is remarkable, however, to see that some fraction gets transported 
very far, during the primordial partonic non-equilibrium phase. We shall turn 
to its theoretical description in section 12.41 but note, for now, that pp collisions 
studied at the CERN ISR [58j lead to a qualitatively similar net baryon rapidity 
distribution, albeit characterized by a smaller (Sy). 

The data described above suggest that the stopping mechanism universally re- 
sides in the primordial, first generation of collisions at the microscopic level. 
The rapidity distributions of charged particle multiplicity, transverse energy 
and valence quark exhibit qualitatively similar shapes (which also evolve sim- 
ilarly with ^/s) in pp, pp, e^e~ reactions, on the one hand, and in central or 
semi-peripheral collisions oi A k, 200 nuclei, on the other. Comparing in de- 
tail we formulate a nuclear modification factor for the bulk hadron rapidity 



^Neutrons are not directly measured in the SPS and RHIC experiments but their produc- 
tion rate, relative to protons, reflects in the ratio of tritium to '^He production measured 
by NA49 [5Z|, applying the isospin mirror symmetry of the corresponding nuclear wave 
functions. 
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2 Bulk Hadron Production in A+A Collisions 
distributions, 



~,AA 



dN^^/dy iy) in A + A 



y ~ 0.5 Npart dN^'^/dy in pp ^ ' ' 

where Npart < 2^ is the mean number of "participating nucleons" (which un- 
dergo at least one inelastic collision with another nucleon) which increases with 

collision centrality. For identical nuclei colliding ( A''p™/ ) — ( ^p'Z-t ) ^^'^ thus 
0.5 Npart gives the number of opposing nucleon pairs. R^^ = 1 if each such 
"opposing" pair contributes the same fraction to the total A+A yield as is 
produced in minimum bias pp at similar y/s. From Figs. 12.31 and 12.51 we infer 
that for I r] I < 1, R = 1.5 at top RHIC energy, and for the pseudo-rapidity 
integrated total N'^^ we find R = 1.36, in central Au+Au collisions. AA 
collisions thus provide for a higher stopping power than pp (which is also re- 
flected in the higher rapidity shift {6y) of Fig. 12. 9p . The observation that their 
stopping power resembles the e~^e~ inelasticity suggests a substantially reduced 
leading particle effect in central collisions of heavy nuclei. This might not be 
surprising. In a Glauber-view of successive minimum bias nucleon collisions 
occuring during interpenetration, each participating nucleon is struck u > 3 
times on average, which might saturate the possible inelasticity, removing the 
leading fragment. 

This view naturally leads to the scaling of the total particle production in 
nuclear collisions with Npart, as seen clearly in Fig. 12.51 reminiscent of the 
"wounded nucleon model" [59] but with the scaling factor determined by e+e" 
rather than pp |60j . Overall we conclude from the still rather close similarity 
between nuclear and elementary collisions that the mechanisms of longitudinal 
phase space population occur primordially, during interpenetration which is 
over after 0.15 fm/c at RHIC, and after 1.5 fm/c at SPS energy. I.e. it 
is the primordial non-equilibrium pQCD shower evolution that accounts for 
stopping, and its time extent should be a lower limit to the formation time 
To employed in the Bjorken model [3S], equation 12.11 Equilibration at the 
partonic level might begin at t > tq only (the development toward a quark- 
gluon-plasma phase), but the primordial parton redistribution processes set 
the stage for this phase, and control the relaxation time scales involved in 
equilibration |61j . More about this in section [231 We infer the existence of a 
saturation scale [62] controlling the total inelasticity: with ever higher reactant 
thickness, proportional to A^'^, one does not get a total rapidity or energy 
density proportional to A'^'^ (the number of " successive binary collisions" ) but 
to A^'^^ only [63]. Note that the lines shown in Fig. 12.61 (right panel) refer to such 
a saturation theory: the color glass condensate (CGC) model [S3] developed by 
McLerran and Venugopulan. The success of these models demonstrates that 
" successive binary baryon scattering" is not an appropriate picture at high ^/s. 
One can free the partons from the nucleonic parton density distributions only 
once, and their corresponding transverse areal density sets the stage for the 
ensuing QCD parton shower evolution [62j. Moreover, an additional saturation 
eff'ect appears to modify this evolution at high transverse areal parton density 
(see chapter [27 



26 



2.3 Dependence on system size 

2.3 Dependence on system size 

We have discussed above a first attempt toward a variable {Npart) that scales the 
system size dependence in A+A collisions. Note that one can vary the size either 
by centrally colliding a sequence of nuclei, Ai + Ai, A2 + A2 etc., or by selecting 
different windows in Npart out of minimum bias collision ensembles obtained 
for heavy nuclei for which BNL employs Au and CERN ^^^Ph. The third 
alternative, scattering a relatively light projectile, such as ^■^S, from increasing 
A nuclear targets, has been employed initially both at the AGS and SPS but got 
disfavored in view of numerous disadvantages, of both experimental (the need 
to measure the entire rapidity distribution, i.e. lab momenta from about 0.3- 
100 GeV /c, with uniform efficiency) and theoretical nature (different density 
distributions of projectile and target; occurence of an"effectiv" center of mass, 
different for hard and soft collisions, and depending on impact parameter). 

The determination of Npart is of central interest, and thus we need to look 
at technicalities, briefly. The approximate linear scaling with Npart that we 
observed in the total transverse energy and the total charged particle number 
(Figs. 12.21 12. 5p is a reflection of the primordial redistribution of partons and 
energy. Whereas all observable properties that refer to the system evolution 
at later times, which are of interest as potential signals from the equilibrium, 
QCD plasma "matter" phase, have different specific dependences on Npart, be 
it suppressions (high pT signals, jets, quarkonia production) or enhancements 
(collective hydrodynamic flow, strangeness production). Npart thus emerges as 
a suitable common reference scale. 

Npart captures the number of potentially directly hit nucleons. It is estimated 
from an eikonal straight trajectory Glauber model as applied to the overlap re- 
gion arising, in dependence of impact parameter h, from the superposition along 
beam direction of the two initial Woods-Saxon density distributions of the in- 
teracting nuclei. To account for the dilute surfaces of these distributions (within 
which the intersecting nucleons might not find an interaction partner) each in- 
cident nucleon trajectory gets equipped with a transverse radius that represents 
the total inelastic NN cross section at the corresponding ^/s. The formalism 
is imbedded into a Monte Carlo simulation (for detail see |66j ) starting from 
random microscopic nucleon positions within the transversely projected initial 
Woods-Saxon density profiles. Overlapping cross sectional tubes of target and 
projectile nucleons are counted as a participant nucleon pair. Owing to the 
statistics of nucleon initial position sampling each considered impact parameter 
geometry thus results in a probability distribution of derived Npart- Its width 
a defines the resolution A (6) of impact parameter b determination within this 
scheme via the relation 

1 ^^^^<i{Nparm ^^^^^ 



A{b) ' ' db 

which, at A=200, leads to the expectation to determine b with about 1.5 fm 
resolution [66] , by measuring Npart ■ 
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How to measure Npart^ In fixed target experiments one can calorimetrically 
count all particles with beam momentum per nucleon and superimposed Fermi 
momentum distributions of nucleons, i.e. one looks for particles in the beam 
fragmentation domain yheam i 0.5, pT < 0.25 GeV/c. These are identified as 
spectator nucleons, and N^^rt = ^ — ^fpel- For identical nuclear collision 

systems (N'^artJ — {^parf/j 9,nd thus Npart gets approximated by 2 Np^"!. 
This scheme was employed in the CERN experiments NA49 and WA80, and 
generalized [67] in a way that is illustrated in Fig. 12.101 
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Figure 2.10: (a) Energy spectrum of the forward calorimeter in Pb+Pb colli- 
sions at 158^ GeV; (b) impact parameter and fraction of total 
inelastic cross section related to forward energy from the VENUS 
model [671 [68]. 



The top panel shows the minimum bias distribution of total energy registered 
in a forward calorimeter that covers the beam fragment domain in Pb+Pb 
collisions at lab. energy of 158 GeV per projectile nucleon, ^/s = 17.3 GeV. 
The energy spectrum extends from about 3TeV which corresponds to about 20 
projectile spectators (indicating a "central" collision), to about 32 TeV which 
is close to the total beam energy and thus corresponds to extremely peripheral 
collisions. Note that the shape of this forward energy spectrum is the mirror 
image of the minimum bias transverse energy distribution of Fig. 12.21 both 
recorded by NA49. From both figures we see that the ideal head-on, 6 — > 
collision can not be selected from these (or any other) data, owing to the facts 
that b = carries zero geometrical weight, and that the diffuse Woods-Saxon 
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nuclear density profiles lead to a fluctuation of participant nucleon number at 
given flnite b. Thus the Npart fluctuation at flnite weight impact parameters 
overshadows the genuinely small contribution of near zero impact parameters. 
Selecting "central" collisions, either by an on-line trigger cut on minimal forward 
energy or maximal total transverse energy or charged particle rapidity density, 
or by corresponding off-line selection, one thus faces a compromise between 
event statistics and selectivity for impact parameters near zero. In the example 
of Fig. 12.101 these considerations suggest a cut at about 8TeV which selects the 
5% most inelastic events, from among the overall minimum bias distribution, 
then to be labeled as "central" collisions. This selection corresponds to a soft 
cutoff at 6 < 3 fm. 

The selectivity of this, or of other less stringent cuts on collision centrality 
is then established by comparison to a Glauber or cascade model. The bot- 
tom panel of Fig. 12.101 employs the VENUS hadron/string cascade model [68] 
which starts from a Monte Carlo position sampling of the nucleons imbed- 
ded in Woods-Saxon nuclear density proflles but (unlike in a Glauber scheme 
with straight trajectory overlap projection) following the cascade of inelastic 
hadron/string multiplication, again by Monte Carlo sampling. It reproduces 
the forward energy data reasonably well and one can thus read off the average 
impact parameter and participant nucleon number corresponding to any desired 
cut on the percent fraction of the total minimum bias cross section. Moreover, 
it is clear that this procedure can also be based on the total minimum bias 
transverse energy distribution. Fig. 12.21 which is the mirror image of the for- 
ward energy distribution in Fig. 12.101 or on the total, and even the mid-rapidity 
charged particle density (Fig. 12. 5|) . The latter method is employed by the RHIC 
experiments STAR and PHENIX. 

How well this machinery works is illustrated in Fig. 12.111 by RHIC-PHOBOS 
results at y/s = 200 GeV [S2] . The charged particle pseudo-rapidity density dis- 
tributions are shown for central (3-6% highest Nch cut) Cu-|-Cu collisions, with 
{]^part) = 100, and semi-peripheral Au-|-Au collisions selecting the cut window 
(35-40%) such that the same (Npart) emerges. The distributions are nearly 
identical. In extrapolation to Npart = 2 one would expect to flnd agreement 
between min. bias p + p, and "super-peripheral" A-l-A collisions, at least at 
high energy where the nuclear Fermi momentum plays no large role. Fig. 12.121 
shows that this expectation is correct [69j . As it is technically difficult to select 
Npart = 2 from A=200 nuclei colliding, NA49 fragmented the incident SPS Pb 
beam to study ^'^C -|-^^ C and '^^Si -|-^^ Si collisions [67j. These systems are 
isospin symmetric, and Fig. 12.121 thus plots 0.5((7r"*") -|- {Tr~))/ {N\y) including 
p + p where Nw = 2 by definition. We see that the pion multiplicity of A-l-A 
collisions interpolates to the p-|-p data point. 

Note that NA49 employs the term "wounded nucleon" number (A^iy) to count 
the nucleons that underwent at least one inelastic nucleon-nucleon collision. 
This is what the RHIC experiments (that follow a Glauber model) call Npart 
whereas NA49 reserves this term for nucleons that underwent any inelastic colli- 
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Figure 2.11: Charged hadron pseudo- Figure 2.12: Charged pion multiphcity 
rapidity distributions in normahzed by Nw vs. cen- 

Cu+Cu and Au+Au colh- trahty in p+p, C+C, Si+Si 

sions at ^/s = 200 GeV, with and Pb+Pb colhsions at 

similar Npart ~ 100 [52j. y/^ = 17.3 GeV [671 [69]. 



sion. Thus Nw in Fig.[2l2]has the same definition as Npart in Figs. [231 [231 [2n[2TTl 
We see that a smooth increase joins the p + p data, via the hght A+A central 
collisions, to a saturation setting in with semi-peripheral Pb+Pb collisions, the 
overall, relative increase amounting to about 40% (as we saw in Fig. 12. 3p . 

1/3 

There is nothing like an N '^ increase (the thickness of the reactants) observed 
here, pointing to the saturation mechanism(s) mentioned in the previous sec- 
tion, which are seen from Fig. 12.121 to dampen the initial, fast increase once 
the primordial interaction volume contains about 80 nucleons. In the Glauber 
model view of successive collisions (to which we attach only symbolical sig- 
nificance at high y/s) this volume corresponds to (i^) ~ 3, and within the 
terminology of such models we might thus argue, intuitively, that the initial ge- 
ometrical cross section, attached to the nucleon structure function as a whole, 
has disappeared at (v) « 3, all constituent partons being freed. 



2.4 Gluon Saturation in A+A Collisions 

We will now take a closer look at the saturation phenomena of high energy QCD 
scattering, and apply results obtained for deep inelastic electron-proton reac- 
tions to nuclear collisions, a procedure that relies on a universality of high en- 
ergy hadron scattering. This arises at high y/s, and at relatively low momentum 
transfer squared Q^ (the condition governing bulk charged particle production 
near mid-rapidity at RHIC, where Feynman x ~ 0.01 and Q^ < 5 GeV'^). Uni- 
versality comes about as the transverse resolution becomes higher and higher, 
with Q^, so that within the small area tested by the collision there is no dif- 
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ference whether the partons sampled there belong to the transverse gluon and 
quark density projection of any hadron species, or even of a nucleus. And sat- 
uration arises once the areal transverse parton density exceeds the resolution, 
leading to interfering QCD sub-amplitudes that do not reflect in the total cross 
section in a manner similar to the mere summation of separate, resolved color 
charges [6ll [62l [Ml EH ES EOl EI] • 

The ideas of saturation and universality are motivated by HERA deep inelastic 
scattering (DIS) data [72] on the gluon distribution function shown in Fig. 12.131 
(left side). The gluon rapidity density, xG{x,Q^) = — ^^^ rises rapidly as 
a function of decreasing fractional momentum, x, or increasing resolution, Q^. 
The origin of this rise in the gluon density is, ultimately, the non-abelian nature 
of QCD. Due to the intrinsic non-linearity of QCD [701 171] . gluon showers 
generate more gluon showers, producing an avalanche toward small x. As a 
consequence of this exponential growth the spatial density of gluons (per unit 
transverse area per unit rapidity) of any hadron or nucleus must increase as x 
decreases |65j . This follows because the transverse size, as seen via the total 
cross section, rises more slowly toward higher energy than the number of gluons. 
This is illustrated in Fig. 12.131 (right side). In a head-on view of a hadronic 
projectile more and more partons (mostly gluons) appear as x decreases. This 
picture reflects a representation of the hadron in the "infinite momentum frame" 
where it has a large light-cone longitudinal momentum P^ ^ M. In this 
frame one can describe the hadron wave function as a collection of constituents 
carrying a fraction p'^ = xP~^, < x < 1, of the total longitudinal momentum 
[73j ("light cone quantization" method [71]). In DIS at large sqrts and Q^ one 
measures the quark distributions dNq/dx at small x, deriving from this the 
gluon distributions xG{x, Q^) of Fig. 12.131 

It is useful |75j to consider the rapidity distribution implied by the parton 
distributions, in this picture. Defining y = yhadron — ln{l/x) as the rapidity of 
the potentially struck parton, the invariant rapidity distribution results as 

dN/dy = X dN/dx = xG{x, Q^). (2.8) 

At high Q^ the measured quark and gluon structure functions are thus simply 
related to the number of partons per unit rapidity, resolved in the hadronic 
wave function. 

The above textbook level |74l [75] recapitulation leads, however, to an impor- 
tant application: the dN/dy distribution of constituent partons of a hadron (or 
nucleus), determined by the DIS experiments, is similar to the rapidity distri- 
bution of produced particles in hadron-hadron or A+A collisions as we expect 
the initial gluon rapidity density to be represented in the finally observed, pro- 
duced hadrons, at high ^/s. Due to the longitudinal boost invariance of the 
rapidity distribution, we can apply the above conclusions to hadron-hadron or 
A+A collisions at high ^/s, by replacing the infinite momentum frame hadron 
rapidity by the center of mass frame projectile rapidity, yproj, while retaining 
the result that the rapidity density of potentially interacting partons grows with 
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Figure 2.13: (left) The HERA data for the gluon distribution function as a 
function of fractional momentum x and square momentum transfer 
Q^ [?2|. (right) Saturation of gluons in a hadron; a head on view 
as X decreases j75j. 



increasing distance from yproj like 

Ay = yproj-y = ln{l/x). (2.9) 

At RHIC energy, ^/s = 200 GeV, Ay at mid-rapidity thus corresponds to 
X < 10^^ (well into the domain of growing structure function gluon density, 
Fig. I2.13p . and the two intersecting partonic transverse density distributions 
thus attempt to resolve each other given the densely packed situation that is 
depicted in the lower circle of Fig. 12.131 (right panel). At given Q^ (which is 
modest, Q^ < 5GeV^, for bulk hadron production at mid-rapidity) the packing 
density at mid-rapidity will increase toward higher ^/s as 



Ay 



midrap 



ln{y/s/M), i.e. l/x ^ ^/M 



(2.10) 



thus sampling smaller x domains in Fig. 12.131 according to equation 12.91 It will 
further increase in proceeding from hadronic to nuclear reaction partners A+A. 
Will it be in proportion to A^'^? We know from the previous sections (12.21 and 
12. 3p that this is not the case, the data indicating an increase with A^'^^. This 
observation is, in fact caused by the parton saturation effect, to which we turn 
now. 

For given transverse resolution Q^ and increasing l/x the parton density of 
Fig. 12.131 becomes so large that one can not neglect their mutual interactions 
any longer. One expects such interactions to produce " shadowing" , a decrease 
of the scattering cross section relative to incoherent independent scattering 
[701 [7T] . As an effect of such shadowed interactions there occurs [75] a saturation 
[SH Ea [Ml El ES EOl EH ES] of the cross section at each given Q^ slowing 
the increase with l/x to become logarithmic once l/x exceeds a certain critical 
value Xs{Q'^). Conversely, for fixed x, saturation occurs for transverse momenta 
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below some critical Q'^{x), 

where dN/dy is the x-dependent gluon density (at y = yproj — ln{l/x)). Q^ 
is called the saturation scale. In equation 12.111 vri?^ is the hadron area (in 
transverse projection), and OgNc is the color charge squared of a single gluon. 
More intuitively, Ql{x) defines an inversely proportional resolution area Fs{x) 
and at each x we have to choose Fs{x) such that the ratio of total area ttR'^ 
to Fs{x) (the number of resolved areal pixels) equals the number of single 
gluon charge sources featured by the total hadron area. As a consequence 
the saturation scale Q'l{x) defines a critical areal resolution, with two different 
types of QCD scattering theory defined, at each x, for Q^ > Q^ and Q^ < Q^, 
respectively [621 [651 [75] . 

As one expects a soft transition between such theories, to occur along the 
transition line implied by Q'^{x), the two types of QCD scattering are best 
studied with processes featuring typical Q^ well above, or below Q'^{x). Jet 
production at ^/s > 200 GeV in pp or AA collisions with typical Q^ above 
about lO^GeV^, clearly falls into the former class, to be described e.g. by QCD 
DGLAP evolution of partonic showers [76]. The accronym DGLAP refers to 
the inventors of the perturbative QCD evolution of parton scattering with the 
"runing" strong coupling constant OsiQ'^), Dokshitzer, Gribov, Levine, Altarelli 
and Parisi. On the other hand, mid-rapidity bulk hadron production at the 
upcoming CERN LHC facihty {y/s = UTeV for pp, and 5.5TeV for A+A), with 
typical Q^ < 5 GeV^ at x < 10"'^, will present a clear case for QCD saturation 
physics, as formulated e.g. in the "Color Glass Condensate (CGC)" formalism 
developed by McLerran, Venugopalan and collaborators [Ml ESl [75l [77] . This 
model develops a classical gluon field theory for the limiting case of a high areal 
occupation number density, i.e for the conceivable limit of the situation depicted 
in Fig. 12.131 (right hand panel) where the amalgamating small x gluons would 
overlap completely, within any finite resolution area at modest Q^. Classical 
field theory captures, by construction, the effects of color charge coherence, 
absent in DGLAP parton cascade evolution theories [75]. This model appears 
to work well already at y^ as "low" as at RHIC, as far as small Q^ bulk 
charged particle production is concerned. We have illustrated this by the CGC 
model fits [64| to the PHOBOS charged particle rapidity distributions, shown 
in Fig. [231 



Conversely, QCD processes falling in the transition region between such limiting 
conditions, such that typical Q^ ~ Ql{x), should present observables that are 
functions of the ratio between the transferred momentum Q^ and the appropri- 
ate saturation scale, expressed by Q'^{x). As Q^ defines the effective transverse 
sampling area, and Q'^{x) the characteristic areal size at which saturation is ex- 
pected to set in, a characteristic behaviour of cross sections, namely that they 
are universal functions of Q^/Q^, is called "geometric scaling''' . The HERA ep 
scattering data obey this scaling law closely [75] , and the idea arises to apply the 
universality principle that we mentioned above: at small enough x, all hadrons 
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or nuclei are similar, their specific properties only coming in via the appropriate 
saturation scales Ql{x, h) or Ql{x, A). Knowing the latter for RHIC conditions 
we will understand the systematics of charged particle production illustrated in 
the previous chapter, and thus also be able to extrapolate toward LHC condi- 
tions in pp and AA collisions. 
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Figure 2.14: (top) Geometric scaling of the virtual photo-absorption cross sec- 
tion a'^P on protons; (middle) cross sections for nuclei normalized 
according to equation I2.13t (bottom) the ratio of a^^ to a fit of 
a"''^ (see [63j for data reference). 



All data for the virtual photo-absorption cross section a'^P{x, Q"^) in deep inelas- 
tic ep scattering with x < 0.01 (which is also the RHIC mid-rapidity x-domain) 
have been found [TS] to lie on a single curve when plotted against Q"^ /Ql, with 



Ql{x) ~ (^)^ IGeV^ 

X 



(2.12) 



with A ~ 0.3 and xq — 10~^. This scaling f79] with r = Q'^ /Ql is shown in 
Fig. 12.141 (top panel) to interpolate all data. A chain of arguments, proposed by 
Armesto, Salgado and Wiedemann [63] connects a fit to these data with photo- 
absorption data for (virtual) photon-A interactions [80j via the geometrical 
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scaling ansatz 

a^^JTA) a-^nrp = ta) ._„. 

7,R\ ~ -kRI ^ ''' 

assuming that the scale in the nucleus grows with the ratio of the transverse 
parton densities, raised to the power 1/5 (a free parameter), 



(^\ .r.^rJl^ 



QiA = QiA^^ '^A = r, -^ . (2.14) 



Fig. 12.141 (middle and bottom panels) shows their fit to the nuclear photo- 
absorbtion data which fixes 6 = 0.79 and irRp = 1.57 fm? (see ref. p3l for 
detail). The essential step in transforming these findings to the case of A+A 
collisions is then taken by the empirical ansatz 

—— {at y ~ 0) oc QIAx)ttR\ (2.15) 

ay 

by which the mid-rapidity parton (gluon) density dN/dy in equation 12.111 gets 
related to the charged particle mid-rapidity density at y ss [TOl [H] , measured 
in nucleus-nucleus collisions. Replacing, further, the total nucleon number 2A 
in a collision of identical nuclei of mass A by the number Npart of participating 
nucleons, the final result is |63] 



1 (JN^^ 

{aty^O) = No{V~s)'Kart (2.16) 



where the exponent a = (1 — 6)/35 = 0.089, and A^o = 0.47. The exponent a 
is far smaller than 1/3, a value that represents the thickness of the reactants, 
and would be our naive guess in a picture of "successive" independent nucleon 
participant collisions, whose average number (z^) oc (Npartf^)^ ■ The obser- 
vational fact (see Fig. 1.13) that a < 1/3 for mid-rapidity low Q^ bulk hadron 
production in A+A collisions illustrates the importance of the QCD satura- 
tion effect. This is shown [63j in Fig. 12.151 where equation 12.161 is applied to 
the RHIC PHOBOS data for mid-rapidity charged particle rapidity density per 
participant pair, in Au+Au collisions at ^/s = 19.6, 130 and 200 GeV [82], also 
including a prediction for LHC energy. Note that the factorization of energy 
and centrality dependence, implied by the RHIC data [52j, is well captured by 
equation l2.11l and the resulting fits in Fig. 12.151 Furthermore, the steeper slope, 
predicted for Npart < 60 (not covered by the employed data set), interpolates 
to the corresponding pp and pp data, at Npart = 2. It resembles the pattern 
observed in the NA49 data (Fig. I2.12p for small Npart collisions of light A+A 
systems, at y/s = 17 — 20 GeV , and may be seen, to reflect the onset of QCD 
saturation. Finally we note that the conclusions of the above, partially heuris- 
tic approach [63] , represented by equations 12.131 to 12.161 have been backed up 
by the CGC theory of McLerran and Venugopulan [6-1] [65] I75j . predictions of 
which we have illustrated in Fig. 12. 6i 
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Figure 2.15: Saturation model fit |63j applied to RHIC charged hadron multi- 
plicity data at mid-rapidity normalized by number of participant 
pairs, at various energies [82]. Also shown is an extrapolation to 
pp data and a prediction for minimum bias Pb+Pb collisions at 
LHC energy, ^ = 5500 GeV. 



Bulk hadron production in AA collisions at high ^/s can be related, via the as- 
sumption of universality of high energy QCD scattering, to the phenomenon of 
geometric scaling first observed in HERA deep inelastic ep cross sections. The 
underlying feature is a QCD saturation effect arising from the diverging areal 
parton density, as confronted with the limited areal resolution Q^, inherent in 
the considered scattering process. The "saturation scale" Q'^{x,A) captures 
the condition that a single partonic charge source within the transverse par- 
tonic density profile can just be resolved by a sufficiently high Q^. Bulk hadron 
production in A-|-A collisions falls below this scale. 



36 



2.5 Transverse phase space: equilibrium and the QGP state 

2.5 Transverse phase space: equilibrium and the QGP 
state 



At RHIC energy, ^/s = 200 GeV, the Au+Au collision reactants are longitudi- 
nally contracted discs. At a nuclear radius R ~ A^'^ fm and Lorentz 7 ~ 100 
their primordial interpenetration phase ends at time tq < 0.15 fm/c. This 
time scale is absent in e^e~ annihilation at similar ^/s where tq ~ 0.1 fni/c 
marks the end of the primordial pQCD partonic shower evolution [83j during 
which the initially created qq pair, of "virtually" Q = y/s/2 each, multiplies 
in the course of the QCD DGLAP evolution in perturbative vacuum, giving 
rise to daughter partons of far lower virtuality, of a few GeV. In A+A colli- 
sions this shower era should last longer, due to the interpenetrational spread of 
primordial collision time. It should be over by about 0.25 fm/c. The shower 
partons in e~^e~ annihilation are localized within back to back cone geometry 
reflecting the directions of the primordial quark pair. The eventually observed 
"jet" signal, created by an initial Q^ of 10^ GeV^, is established by then. Upon 
a slow-down of the dynamical evolution time scale to r ~ 1 fmjc the shower 
partons fragment further, acquiring transverse momentum and yet lower vir- 
tuality, then to enter a non perturbative QCD phase of color neutralization 
during which hadron-like singlet parton clusters are formed. Their net initial 
pQCD virtuality, in pQCD vacuum, is recast in terms of non-perturbative vac- 
uum hadron mass. The evolution ends with on-shell, observed jet-hadrons after 
about 3 fmjc of overall reaction time. 

Remarkably, even in this, somehow most elementary process of QCD evolution, 
an aspect of equilibrium formation is observed, not in the narrowly focussed final 
dijet momentum topology but in the relative production rates of the various 
created hadronic species. This so-called "hadro chemical" equilibrium among 
the hadronic species is documented in Fig. 12.161 The hadron multiplicities per 
e"'"e~ annihilation event at \fs = 91.2 GeV [3S] are confronted with a Hage- 
dorn [38] canonical statistical Gibbs ensemble prediction [8l] which reveals that 
the apparent species equilibrium was fixed at a temperature of T = 165 MeV, 
which turns out to be the universal hadronization temperature of all elemen- 
tary and nuclear collisions at high y/s (Hagedorns limiting temperature of the 
hadronic phase of matter). We shall return to this topic in chapter [3] but note, 
for now, that reactions with as few as 20 charged particles exhibit such statisti- 
cal equilibrium properties, a pre-requisite for application of thermodynamic or 
hydrodynamic concepts. 

What happens with parton (and hadron) dynamics in A+A collisions after 
tq? There will not be a QCD evolution in vacuum (which would be over after 
3 fm/c) as the transverse radius of the interacting system is large. It may grow 
to about twice the nuclear radius, i.e. to about 15 fm before interactions cease; 
i.e. the system needs about 15 fmjc to decouple. This simple fact is the key to 
our expectation that the expansive evolution of the initial high energy density 
deposited in a cylinder of considerable diameter (about 10 /?ti), may create 
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Figure 2.16: Hadron multiplicities in LEP e^e^ annihilation at y/s = 91.2GeV 
confronted with the predictions of the canonical statistical 
hadronization model 



certain equilibrium properties that allow us to treat the contained particles and 
energy in terms of thermodynamic phases of matter, such as a partonic QGP 
liquid, or a hadronic liquid or gas, etc.. Such that the expansion dynamics 
makes contact to the phase diagram illustrated in Fig. 11.11 This expectation 
turns out to be justified as we shall describe in chapters [3] and [H What results 
for the evolution after tq in a central A+A collision is sketched in Fig. 12.171 by 
means of a schematic 2-dimensional light cone diagram, which is entered by the 
two reactant nuclei along z = t trajectories where z is the beam direction and 
Lorentz contraction has been taken to an extreme, such that there occurs an 
idealized t = z = interaction "point". Toward positive t the light cone proper 
time profiles of progressing parton-hadron matter evolution are illustrated. The 
first profile illustrated here corresponds to the end of formation time tq. From 
our above discussion of the e~'"e~ annihilation process one obtains a first esti- 
mate. To > 0.25 fra/c (including interpenetration time of 0.15 fm/c at RHIC) 
which refers to processes of very high Q^ > 10^ GeV"^, far above the saturation 
scale Qg discussed in the previous section. The latter scale has to be taken into 
account for low pT hadron production. 

It is the specific resolution scale Q^ of a QCD sub-process, as enveloped in 
the overall collision dynamics of two slabs of given transverse partonic struc- 
ture function density, that determines which fraction of the constituent partons 
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Figure 2.17: Schematic light cone diagram of the evolution of a high energy 
heavy ion collision, indicating a formation phase tq (see text). 



enters interaction. In the simple case of extremely high Q^ processes the an- 
swer is that all constituents are resolved. However, at modest Q^ (dominating 
bulk hadron production) the characteristic QCD saturation scale Q'^{x) gains 
prominence, defined such that processes with Q^ < Q^ do not exploit the initial 
transverse parton densities at the level of independent single constituent color 
field sources (see equation l2.1ip . For such processes the proper formation time 
scale, To, is of order of the inverse saturation momentum [61], 1/Qs ~ 0.2 fm/c 
at ^/s = 200 GeV. The first profile of the time evolution, sketched in Fig. \2.n\ 
should correspond to proper time t = tq = 0.25 fm/c at RHIC energy. At 
top SPS energy, ^/s = 17.3 GeV, we can not refer to such detailed QCD con- 
siderations. A pragmatic approach suggests to take the interpenetration time, 
at 7 ~ 8.5, for guidance concerning the formation time, which thus results as 
To ss 1.5 fm/c. 

In summary of the above considerations we assume that the initial partonic 
color sources, as contained in the structure functions (Fig. 12.13]) . are spread out 
in longitudinal phase space after light cone proper time t = tq « 0.2 fm/c, at 
top RHIC energy, and after tq « 1.5 fm/c at top SPS energy. No significant 
transverse expansion has occured at this early stage, in a central collision oi Ak, 
200 nuclei with transverse diameter of about 12 /ttt,. The Bjorken estimate |45] 
of initial energy density e (equation 12. ip refers to exactly this condition, after 
formation time tq. In order to account for the finite longitudinal source size and 
interpenetration time, at RHIC, we finally put the average tq ~ 0.3 fm, at y/s = 
200 GeV, indicating the "initialization time" after which all partons that have 
been resolved from the structure functions are engaged in shower multiplication. 
As is apparent from Fig. 12.171 this time scale is Lorentz dilated for partons 
with a large longitudinal momentum, or rapidity. This means that the slow 
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particles are produced first toward the center of the colhsion region, and the 
fast (large rapidity) particles are produced later, away from the collision region. 
This Bjorken "inside-out" correlation [15] between coordinate- and momentum- 
space is similar to the Hubble expansion pattern in cosmology: more distant 
galaxies have higher outward velocities. This means that the matter created in 
A-|-A collisions at high y/s is also born expanding, however with the difference 
that the Hubble flow is initially one dimensional along the collision axis. This 
pattern will continue, at y/s = 200 GeV, until the system begins to feel the 
effects of finite size in the transverse direction which will occur at some time to 
in the vicinity of 1 fm/c. However, the tight correlation between position and 
momentum initially imprinted on the system will survive all further expansive 
evolution of the initial "firetube" , and is well recovered in the expansion pattern 
of the finally released hadrons of modest pT as we shall show when discussing 
radial flow and pion pair Bose-Einstein momentum correlation (see section 12.61 
and chapter [Tj) . 

In order to proceed to a more quantitative description of the primordial dynam- 
ics (that occurs onward from tq for as long the time period of predominantly 
longitudinal expansion might extend) we return to the Bjorken estimate of en- 
ergy density, corresponding to this picture [35], as implied by equation 12.1^ 
which we now recast as 



dy J 



'EDi-nRlt,)-^ (2.17) 



where the first term is the (average) total hadron multiplicity per unit rapid- 
ity which , multiplied with the average hadron transverse energy, equals the 
total transverse energy recorded in the calorimetric study shown in Fig. 12.21 
as employed in equation 12.11 The quantity Ra is, strictly speaking, not the 
radius parameter of the spherical Woods-Saxon nuclear density profile but the 
rms of the reactant overlap profiles as projected onto the transverse plane (and 
thus slightly smaller than Ra ~ A^'^ I'm)- Employing A^'^ here (as is done 
throughout) leads to a conservative estimate of e, a minor concern. However, 
the basic assumption in equation 12.171 is to identify the primordial transverse 
energy " radiation" , of an interactional cylindric source of radius Ra and length 
to (where To < to < 1 fm/c, not Lorentz dilated at midrapidity) , with the 
finally emerging bulk hadronic transverse energy. We justify this assumption 
by the two observations, made above, that 

1. the bulk hadron multiplicity density per unit rapidity -j-^ resembles the 
parton density, primordially released at saturation scale tq (Figs. [T6ll2.15p 
at ^ = 200 GeV, and that 

2. the global emission pattern of bulk hadrons (in rapidity and pt) closely 
refiects the initial correlation between coordinate and momentum space, 
characteristic of a primordial period of a predominantly longitudinal ex- 
pansion, as implied in the Bjorken model. 

Both these observations are surprising, at first sight. The Bjorken model was 
conceived for elementary hadron collisions where the expansion proceeds into 
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vacuum, i.e. directly toward observation. Fig. 12.171 proposes that, to the con- 
trary, primordially produced partons have to transform through further, succes- 
sive stages of partonic and hadronic matter, at decreasing but stih substantial 
energy density, in central A+A collisions. The very fact of high energy den- 
sity, with implied short mean free path of the constituent particles, invites a 
hydrodynamic description of the expansive evolution. With initial conditions 
fixed between tq and Iq, an ensuing 3-dimensional hydrodynamic expansion 
would preserve the primordial Bjorken-type correlation between position and 
momentum space, up to lower density conditions and, thus, close to emission of 
the eventually observed hadrons. We thus feel justified to employ equation 12.11 
or 12.171 for the initial conditions at RHIC, obtaining |6H I84j 



6 GeV/fm^ < e < 20 GeV/fm^ (2.18) 

for the interval 0.3 fm/c < tQ < 1 fm/c, in central Au+Au collisions at y ~ 
and ^/s = 200 GeV. The energy density at top SPS energy, ^/s = 17.3 GeV, 
can similarly be estimated [l3l Hlj to amount to about 3 GeV/ fw? at a io 
of 1 fm/c but we can not identify conditions at tq < to in this case as the 
mere interpenetration of two Pb nuclei takes 1.4 fm/c. Thus the commonly 
accepted io = 1 fm/c may lead to a high estimate. An application of the 
parton-hadron transport model of Ellis and Geiger j85t [86] to this collision 
finds e = 3.3 GeV/ fm? at t = 1 fm/c. A primordial energy density of about 
3 GeV/fm^ is twenty times po ~ 0.15 GeV/fm^, the average energy density 
of ground state nuclear matter, and it also exceeds, by far, the critical QCD 
energy density, of 0.6 < ec < 1 GeV/fm'^ according to lattice QCD |48j . The 
initial dynamics thus clearly proceeds in a deconfined QCD system also at top 
SPS energy, and similarly so with strikingly higher energy density, at RHIC, 
where time scales below 1 fm/c can be resolved. 

However, in order now to clarify the key question as to whether, and when 
conditions of partonic dynamical equilibrium may arise under such initial con- 
ditions, we need estimates both of the proper relaxation time scale (which will, 
obviously, depend on energy density and related collision frequency), and of 
the expansion time scale as governed by the overall evolution of the collision 
volume. Only if T{relax.) < T{expans.) one may conclude that the "decon- 
fined partonic system" can be identified with a "deconfined QGP state of QCD 
matter" as described e.g. by lattice QCD, and implied in the phase diagram of 
QCD matter suggested in Fig. II. li 



For guidance concerning the overall time-order of the system evolution we con- 
sider information [87] obtained from Bose-Einstein correlation analysis of pion 
pair emission in momentum space (for detail see chapter [7]) . Note that pions 
should be emitted at any stage of the evolution, after formation time, from the 
surface regions of the evolving "fire-tube". Bulk emission of pions occurs, of 
course, after hadronization (the latest stages illustrated in the evolution sketch 
given in Fig. I2.17|) . The dynamical pion source expansion models by Heinz |88] 
and Sinyukov [89j elaborate a Gaussian emission time profile, with mean Tf 
(the decoupling time) and width At (the duration of emission) . 
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20 

X (fm/c) 



Figure 2.18: Time profile of pion decoupling rate from the fireball in a cen- 
tral Pb+Pb collision, with r = the end of the formation phase. 
Bose-Einstein correlation of vr^vr^pairs yields an average Gaussian 
decoupling profile with tj = 8 fm/c and duration of emission pa- 
rameter At = 4 fm/c [871188]. 



Fig. 12.181 shows an application of this analysis to central Pb+Pb collision 
negative pion pair correlation data obtained by NA49 at top SPS energy, 
^/s = 17.3 GeV [90], where r/ ss 8 fm/c and Ar ss 4 fm/c (note that r = 
in Fig. 12.181 corresponds, not to interaction time t = but to i ~ 1.4 fm/c^ 
the end of the interpenetration phase). We see, first of all, that the overall 
dynamical evolution of a central Pb+Pb collision at \fs = 17.3 GeV is ending 
at about 15 fm,/c\ the proper time defines the position of the last, decoupling 
profile illustrated in Fig. 12.171 for the SPS collisions considered here. While 
the details of Fig. 12.181 will turn out to be relevant to our later discussion of 
hadronization (chapter [3]) and hadronic expansion (chapter H]), we are concerned 
here with the average proper time at which the partonic phase ends. After 
consideration of the duration widths of these latter expansion phases [861 EZ] 
one arrives at an estimate for the average time, spent before hadronization, of 
Ai = 3 — 4 fm/c, again in agreement with the parton cascade model mentioned 
above ^86j. This model also leads to the conclusion that parton thermal equi- 
librium is, at least, closely approached locally in these central Pb+Pb collisions 
as far as mid-rapidity hadron production is concerned (at forward-backward 
rapidity the cascade re-scattering processes do not suffice, however). 

This finding agrees with earlier predictions of Treiax = 1 — 2 fm,/c at top SPS 
energy [9l]. However we note that all such calculations employ perturbative 
QCD methods, implying the paradoxical consequence that equilibrium is closely 
approached only at the end of the partonic phase, at such low ^/s, i.e. in a 
QGP state at about T = 200 MeV which is, by definition, of non-perturbative 
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nature. We shall return to the question of partonic equilibrium attainment at 
SPS energy in the discussion of the hadronization process in nuclear collisions 
(chapter [3]) . 

Equilibrium conditions should set in earlier at top RHIC energy. As transverse 
partonic expansion should set in after the proper time interval 0.3 fm/c <to< 
Ifm/c (which is now resolved by the early dynamics, unlike at top SPS energy), 
we take guidance from the Bjorken estimate of primordial energy density which 
is based on transverse energy production data. Conservatively interpreting the 
result in equation 12.181 we conclude that e is about four times higher than at 
i/s = 17.3 GeV in the above proper time interval. As the binary partonic 
collision frequency scales with the square of with the square of the density p 
(related to the energy density e via the relation e = (E) p = Tp), and is inversely 
proportional to the relaxation time Tj-e^aa; we expect 

Trelax OC (l/pf « (T/ef (2.19) 

which implies that TreiaxiRHIC) ss 0.25 Treiax{SPS) w 0.5 fm/c if we employ 
the estimate T{RHIC) = 2T{SPS). This crude estimate is, however, confirmed 
by the parton transport model of Molar and Gyulassy |92| . 

Partonic equilibration at -y/s = 200 GeV should thus set in at a time scale 
commensurate to the (slightly smaller) formation time scale, at which the to be 
participant partons are resolved from the initial nucleon structure functions and 
enter shower multiplication. Extrapolating to the conditions expected at LHC 
energy (^/s = 5.5 TeV for A+A collisions), where the initial parton density of 
the structure functions in Fig. 12.131 is even higher (x ~ 10~^ at mid-rapidity), 
and so is the initial energy density, we may expect conditions at which the 
resolved partons are almost " born into equilibrium" . 

Early dynamical local equilibrium at RHIC is required to understand the obser- 
vations concerning elliptic flow, with which we shall deal, in detail, in chapterfU 
This term refers to a collective anisotropic azimuthal emission pattern of bulk 
hadrons in semi-peripheral collisions, a hydrodynamical phenomenon that orig- 
inates from the initial geometrical non-isotropy of the primordial interaction 
zone |93[ I94| . A detailed hydrodynamic model analysis of the corresponding 
elliptic flow signal at RHIC [95j leads to the conclusion that local equilibrium 
(a prerequisite to the hydrodynamic description) sets in at to ~ O.dfm/c. This 
conclusion agrees with the estimate via equation 12.191 above, based on Bjorken 
energy density and corresponding parton collisions frequency. 

We note that the concept of a hydrodynamic evolution appears to be, almost 
necessarily ingrained in the physics of a system born into (Hubble-type) expan- 
sion, with a primordial correlation between coordinate and momentum space, 
and at extreme initial parton density at which the partonic mean free path 
length A is close to the overall spatial resolution resulting from the saturation 
scale, i.e. A « l/Qs- 

The above considerations suggest that a quark-gluon plasma state should be 
created early in the expansion dynamics at ^/s = 200 GeV, at about T = 
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300 MeV, that expands hydro dynamically until hadronization is reached, at 
T ~ 165 — 170 MeV. Its manifestations will be considered in chapters [3] to [6j 
At the lower SPS energy, up to 17.3 GeV, we can conclude, with some caution, 
that a deconfined hadronic matter system should exist at T ~ 200 MeV, in 
the closer vicinity of the hadronization transition. It may closely resemble the 
QGP state of lattice QCD, near Tc. 

2.6 Bulk hadron transverse spectra and radial expansion 
flow 

In this chapter we analyze bulk hadron transverse momentum spectra obtained 
at SPS and RHIC energy, confronting the data with predictions of the hydro- 
dynamical model of collective expansion matter flow that we have suggested in 
the previous section, to arise, almost necessarily, from the primordial Hubble- 
type coupling between coordinate and momentum space that prevails at the 
onset of the dynamical evolution in A+A collisions at high ^/s. As all hadronic 
transverse momentum spectra initially follow an approximately exponential fall- 
off (see below) the bulk hadronic output is represented by thermal transverse 
spectra at pT < 2GeV/c. We shall turn to high px information in later sections. 

Furthermore we shall focus here on mid-rapidity production in near central 
A-|-A collisions, because hydrodynamic models refer to an initialization period 
characterized by Bjorken-type longitudinal boost invariance, which we have seen 
in Figs. 12.61 and 12.81 to be restricted to a relatively narrow interval centered at 
mid-rapidity. Central collisions are selected to exploit the azimuthal symmetry 
of emission, in an ideal impact parameter b ^ geometry. We thus select 
the predominant, relevant hydrodynamic "radial flow" expansion mode, from 
among other, azimuthaly oriented (directed) flow patterns that arise once this 
cylindrical symmetry (with respect to the beam direction) is broken in finite 
impact parameter geometries. 

In order to define, quantitatively, the flow phenomena mentioned above, we 
rewrite the invariant cross section for production of hadron species i in terms 
of transverse momentum, rapidity, impact parameter b and azimuthal emission 
angle ipp (relative to the reaction plane), 

^£^, = 2^ S [' + ^"- <^- ""'^^ + '^'' ^'^' ""•(^^'' + -1 

(2.20) 
where we have expanded the dependence on ipp into a Fourier series. Due to 
reflection symmetry with respect to the reaction plane in collisions of identical 
nuclei, only cosine terms appear. Restricting to mid-rapidity production all 
odd harmonics vanish, in particular the "directed flow" coefficient v\, and we 
have dropped the y-dependence in the flow coefficients vl and V2- The latter 
quantifies the amount of " elliptic flow" , to which we turn in chapter HI In 
the following, we will restrict to central collisions which we shall idealize as 
near-zero impact parameter processes governed by cylinder symmetry, whence 
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all azimuthal dependence (expressed by the v\, V2, ■■■ terms) vanishes, and the 
invariant cross section reduces to the first term in equation 12.201 which by 
definition also corresponds to all measurements in which the orientation of the 
reaction plane is not observed. 

Typical transverse momentum spectra of the latter type are shown in Fig. 12.191 
for charged hadron production in Au+Au collisions at ^/s = 200GeV, exhibiting 
mid-rapidity data at various collision centralities ^7\. We observe a clear- 
cut transition, from bulk hadron emission at pT ^ 2 GeV/c featuring a near- 
exponential cross section (i.e. a thermal spectrum), to a high pT power-law 
spectral pattern. Within the context of our previous discussion (section 12. 4p 
we tentatively identify the low pT region with the QCD physics near saturation 
scale. Hadron production at px -^ 10 GeV/c should, on the other hand, be 
the consequence of primordial leading parton fragmentation originating from 
"hard", high Q^ perturbative QCD processes. 

We thus identify bulk hadron production at low pT as the emergence of the ini- 
tial parton saturation conditions that give rise to high energy density and small 
equilibration time scale, leading to a hydrodynamical bulk matter expansion 
evolution. Conversely, the initially produced hard partons, from high Q^ pro- 
cesses, are not thermalized into the bulk but traverse it, as tracers, while being 
attenuated by medium-induced rescattering and gluon radiation, the combined 
effects being reflected in the high pT inclusive hadron yield, and in jet corre- 
lations of hadron emission. We shall turn to the latter physics observables in 
chapter while staying here with low pT physics, related to hydrodynamical 
expansion modes, focusing on radially symmetric expansion. 

In order to infer from the spectral shapes of the hadronic species about the ex- 
pansion mechanism, we first transform to the transverse mass variable, mj- = 
{p^ -l-m^)^/^, via 

1 dNi 1 dNi , , 

— = (2.21) 

2iT pTdprdy 2tt mTdrriTdy 

because it has been shown in p-|-p collisions [98j near RHIC energy that the 
rriT distributions of various hadronic species exhibit a universal pattern {"rriT 
scaling") at low tjit- 

V- T\ = ^^ expi-myT) (2.22) 

with a universal inverse slope parameter T and a species dependent normal- 
ization factor A. Hagedorn showed [99] that this scaling is characteristic of an 
adiabatic expansion of a fireball at temperature T. We recall that, on the other 
hand, an ideal hydrodynamical expansion is isentropic. 

Fig. 12.201 shows the ^/s dependence of the average transverse kinetic energy 
(^mlp'j — 771* for pions, kaons and protons observed at mid-rapidity in central 
Au-|-Au/Pb-|-Pb collisions [54J- Similarly, the inverse slope parameter T re- 
sulting from a fit of equation 12.221 to K^ and K^ transverse mass spectra (at 
Pt ^ 2 GeV/c) is shown in Fig. 12. 2H both for nuclear and p+p collisions [100) . 
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Figure 2.19: Transverse momentum spectra of charged hadrons in Au+Au col- 
lisions at y^ = 200 GeV, in dependence of collision centrality [97] 
(offset as indicated), featuring transition from exponential to power 
law shape. 



We see, first of all, that niT scaling does not apply in A+A collisions, and that 
the kaon inverse slope parameter, T ~ 230 MeV over the SPS energy regime, 
can not be identified with the fireball temperature at hadron formation which 
is Th ~ 165 MeV from Fig. 11.11 The latter is seen, however, to be well repre- 
sented by the p-|-p spectral data exhibited in the left panel of Fig. I2.21i There 
is, thus, not only thermal energy present in A-|-A transverse expansion, but also 
hydrodynamical radial flow. 



We note that the indications in Figs. 12.201 and 12. 2H of a plateau in both (ttzt) 
and T, extending over the domain of SPS energies, 6 < ^/s < 17 GeV, have not 
yet been explained by any fundamental expansive evolution model, including hy- 
drodynamics. Within the framework of the latter model, this is a consequence 
of the initialization problem [UB] which requires a detailed modeling, both of 
primordial energy density vs. equilibration time scale, and of the appropriate 
partonic matter equation of state (EOS) which relates expansion pressure to 
energy density. At top RHIC energy, this initialization of hydro-flow occurs, 
both, at a time scale to ^ 0.6 fm/c which is far smaller than the time scale of 
eventual bulk hadronization (t ~ 3 fm/c), and at a primordial energy density 
far in excess of the critical QCD confinement density. After initialization, the 
partonic plasma phase thus dominates the overall expansive evolution, over a 
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Figure 2.20: The average transverse kinetic energy {mj-) — rriQ for pions, kaons 
and protons vs. y/s in central Au+Au/Pb+Pb collisions |54j . 
Open symbols represent negative hadrons. 



time interval far exceeding the formation and relaxation time scale. 

Thus, at RHIC energy, parton transport [92] and relativistic hydrodynamic 
[95l l96] models establish a well developed expansion mode that survives the 
subsequent stages of hadronization and hadronic expansion. This is reflected in 
their success in describing elliptic flow. On the other hand, the hydro dynamical 
model far overestimates elliptic flow at SPS energy [96] at which, as we have 
shown in sectioiESl the initialization period may be not well separated from 
the confinement (hadronization) stage. Thus, whereas the expansion evolution 
at ^/s = 200 GeV (occuring at near-zero baryo-chemical potential in Fig. II. ip 
"races" across the parton-hadron phase boundary with fully established flow 
patterns, near fis = where lattice QCD predicts the phase transformation 
to be merely a soft cross-over [1^, the dynamics at ^/s = 10 — 20 GeV may 
originate from only slightly above, or even at the phase boundary, thus sam- 
pling the domain 200 < fiB ^ 500 MeV where the equation of state might 
exhibit a "softest point" [96]. The hydrodynamic model thus faces formidable 
uncertainties regarding initialization at SPS energy. 



The plateau in Figs. 12.201 12.211 may be the consequence of the fact that not 
much flow is generated in, or transmitted from the partonic phase, at SPS 
energies, because it is initialized close to the phase boundary ^lOOj where the 
expected critical point [HI dU] (Fig. II. 1|) . and the corresponding adjacent first 
order phase transition might focus [lOlj or stall |96j the expansion trajectory, 
such that the observed radial flow stems almost exclusively from the hadronic 
expansion phase. The SPS plateau, which we shall subsequently encounter in 
other bulk hadron variables (elliptic flow, HBT radii) might thus emerge as a 
consequence of the critical point or, in general, of the flatness of the parton- 
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Figure 2.21: The inverse slope parameter T of equation 12.221 for X^and K^ 
transverse mass spectra at pr < 2 GeV/c and mid-rapidity in cen- 
tral A-|-A, and in minimum bias p-|-p collisions [lOOj . 



hadron coexistence line. RHIC dynamics, on the other hand, originates from 
far above this line. 



Hadronic expansion is known to proceed isentropically |102j : commensurate to 
expansive volume increase the momentum space volume must decrease, from 
a random isotropic thermal distribution to a restricted momentum orientation 
preferentially perpendicular to the fireball surface, i.e. radial. The initial ther- 
mal energy, implied by the hadron formation temperature Th = 165 MeV, 
will thus fall down to a residual Tp at hadronic decoupling from the flow field 
("thermal freeze-out") plus a radial transverse kinetic energy term ttij (/3r) 
where m, is the mass of the considered hadron species and {Pt) the average 
radial velocity. We thus expect [103 j for the slope of equation 12.221 



T = TF + mi {Pt) ,Pt<2 GeV/c 



and 






mi 



(2.23) 



(2.24) 



the latter expression valid at pT larger than hadron mass scale (T then is the 
"blue-shifted temperature" at decoupling |1U4| and (vt) the average transverse 
velocity). The assumption that radial flow mostly originates from the hadronic 
expansion phase is underlined by the proportionality of flow energy to hadron 
mass (equation 12. 23p . 

Fig. 12. 22] illustrates this proportionality, by a recent compilation |103j of RHIC 
results for central Au-|-Au collisions at y/s = 200 GeV, and SPS results for cen- 
tral Pb-|-Pb collisions at top SPS energy, ^/s = 17.3 GeV. At the latter energy 
the slope parameter of the $ meson is seen to be close to that of the similar mass 
baryons p and A, emphasizing the occurence of rrii scaling as opposed to valence 
quark number scaling that we will encounter in RHIC elliptic flow data [M]. As 
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Figure 2.22: Hadron slope parameters T at mid-rapidity as a function of mass. 
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200 GeV (circles); from [ T03] . 



is obvious from Fig. 12.221 the multi-strange liyperons and charmonia exhibit a 
slope saturation which is usually explained |103j as a consequence of their small 
total cross sections of rescattering from other hadrons, leading to an early de- 
coupling from the bulk hadron radial flow field, such that {Pt)q < {Pt)p- 

According to our observations with equation 12.231 a hydrodynamical ansatz 
for the transverse mass spectrum of hadrons should thus contain the variables 
"true temperature" Tp at decoupling from the flow field, and its average veloc- 
ity {Pt), common to all hadrons. This is the case for the blast wave model |104] 
developed as an approximation to the full hydrodynamic formalism [96] . assum- 
ing a common decoupling or "freeze-out" from flow, for all hadronic species, and 
a boost-invariant longitudinal expansion: 



dN,, 



mTdniTdy 



Ai niT Ki 



ruTCOshp 
Tf 



In 



pj-sinhp 
Tf 



(2.25) 



where p = tanh^^j3T- In an extended version of this model a function is in- 
cluded that describes the radial profile of the transverse velocity field, /3t(^) = 
j^max j-jji^ instead of employing a fixed /3t at decoupling [106]. Fig. 12.231 shows 
[5^ the resulting energy dependence of Tf and (/?t), for the same set of data 
as implied already in Figs. 12.20] and YI7IT\ The "true" decoupling temperature 
rises steeply at the AGS and less so at SPS energy (as does (/3t)), to a value of 
about 95MeV at top SPS energy, which is considerably lower than the chemical 
freeze-out temperature, Th = 165 MeV, at which the hadronic species relative 
yield composition of the hadronic phase becomes stationary (see chapter O and 
Fig. ll.ip . Chemical decoupling thus occurs early, near the parton-hadron phase 
boundary, whereas hadronic radial flow ceases after significant further expan- 
sion and cooling, whence the surface radial velocity (its average value given by 
{(3t) in Fig. I2.23P approaches (3t ~ 0.65. Both data sets again exhibit an indi- 
cation of saturation, over the interval toward top SPS energy: the SPS plateau. 
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Figure 2.23: Hadron decoupling temperature Tf, and average radial flow veloc- 
ity {/St) extracted from blast wave model (see equation r2.25p fits 
of niT spectra vs. -y/s 



This supports our above conjecture that radial flow is, predominantly, a conse- 
quence of isentropic bulk hadronic expansion in this energy domain, which sets 
in at Th- At RHIC energy, both parameters exhibit a further rise, suggesting 
that primordial partonic flow begins to contribute significantly to radial flow. 
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3 Hadronization and hadronic 
freeze-out in A+A collisions 



Within the course of the global expansion of the primordial reaction volume 
the local flow " cells" will hit the parton-hadron phase boundary as their energy 
density approaches ecrU ~ 1 GeV/ fm? . Hadronization will thus occur, not at 
an instant over the entire interaction volume, but within a finite overall time 
interval [86] that results from the spread of proper time at which individual cells, 
or coherent clusters of such cells (as developed during expansion) arrive at the 
phase boundary. However, irrespective of such a local-temporal occurence, the 
hadronization process (which is governed by non perturbative QCD at the low 
Q^ corresponding to bulk hadronization) universally results in a novel, global 
equilibrium property that concerns the relative abundance of produced hadrons 
and resonances. This so-called "hadrochemical equilibrium state" is directly 
observable, in contrast to the stages of primordial parton equilibration that are 
only indirectly assessed, via dynamical model studies. 

This equilibrium population of species occurs both in elementary and nuclear 
collisions |107| . We have seen in Fig. 12.161 a first illustration, by e^e~ anni- 
hilation data at ^/s = 91.2 GeV LEP energy, that are well reproduced by the 
partition functions of the statistical hadronization model (SHM) in its canonical 
form [M]. The derived hadronization temperature, Th = 165MeV, turns out to 
be universal to all elementary and nuclear collision processes at ^/s > 20 GeV, 
and it agrees with the limiting temperature predicted by Hagedorn [38j to occur 
in any multi-hadronic equilibrium system once the energy density approaches 
about 0.6 GeV I fw? . Thus, the upper limit of hadronic equilibrium density 
corresponds, closely, to the lower limit, tcrit = 0.6 — 1.0 GeV/ fm'^ of partonic 
equilibrium matter, according to lattice QCD [38]. In elementary collisions 
only about 20 partons or hadrons participate: there should be no chance to 
approach thermodynamic equilibrium of species by rescattering cascades, nei- 
ther in the partonic nor in the hadronic phase. The fact that, nevertheless, the 
hadron formation temperature Th coincides with the Hagedorn limiting tem- 
perature and with the QCD confinement temperature, is a consequence of the 
non-perturbative QCD hadronization process itself [H^, which "gives birth" to 
hadrons/resonances in canonical equilibrium, at high ^/s, as we shall see below. 
This process also governs A+A collisions but, as it occurs here under conditions 
of high energy density extended over considerable volume, the SHM descrip- 
tion now requires a grand canonical ensemble, with important consequences for 
production of strange hadrons (strangeness enhancement). 

The grand canonical order of hadron/resonance production in central A+A col- 
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lisions, and its characteristic strangeness enhancement shows that a state of ex- 
tended matter that is quantum mechanically coherent must exist at hadroniza- 
tion |87ll88t[l07j . Whether or not it also reflects partonic equilibrium properties 
(including flavor equilibrium), that would allow us to claim the direct obser- 
vation of a quark gluon plasma state near Tc, can not be decided on the basis 
of this observation alone, as the hadronization process somehow generates, by 
itself, the observed hadronic equilibrium. This conclusion, however, is still the 
subject of controversy |107j . 
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Figure 3.1: Total hadron multiplicities in central Pb-|-Pb collisions at ^/s = 
17.3 GeV |1UUJ versus prediction of the grand canonical statistical 
hadronization model 1191. 



Two typical examples of grand canonical SHM application are illustrated in 
Figs. l3.l1 and l3.2l the first showing total hadron multiplicities in central Pb-|-Pb 
collisions at ^/s = 17.3 GeV by NA49 [lOOj confronted with SHM predictions 
by Becattini et al. [19]. This plot is similar to Fig. 12.161 in which e~^e~ an- 
nihilation to hadrons is confronted with a SHM prediction derived from the 
canonical ensemble [84] , Central Au+Au collision data at ^/s = 200 GeV 
from several RHIC experiments are compared to grand canonical model predic- 
tions by Braun-Munzinger et al. |108j in Fig. 13.21 The key model parameters, 
Th and the baryo-chemical potential /ig result as 159 MeV (160 MeV), and 
247 MeV (20 MeV) at ./s = 17.3 (200) GeV, respectively The universality of 
the hadronization temperature is obvious from comparison of the present values 
with the results of the canonical procedure employed in e^e~ annihilation to 
hadrons at ^ = 91.2 GeV (Fig. \TTE\i . and in canonical SHM fits [TOgJ to p+p 
collision data at ^/s = 27.4 GeV where Th = 159 and 169 MeV, respectively. 



52 



o 

■« 1 

DC 



10 



10 



10 



o 



A 



Vsnn=200 GeV 



o 



STAR 

PHENIX 

BRAHMS 

- T=160.5, ^t,=20MeV 

- T=155, ^j,=26MeV 



o 



p 



(? 



XL 



^. 



"^ 



n* K* p A E Q. n^ n n n n n p p K K 



A A" 
A p 



Figure 3.2: Hadron multiplicity ratios at mid-rapidity in central Au+Au colli- 
sions at ^ = 200 GeV from RHIC experiments STAR, PHENIX 
and BRAHMS, compared to predictions of the grand canonical sta- 
tistical model IIT08]. 



Figs. l3.11 and l3.2l illustrate two different approaches employed in grand canonical 
SHM application, the former addressing the values of the hadronic multiplici- 
ties as obtained in approximate full 47r acceptance (within limitations implied 
by detector performance), the latter employing a set of multiplicity ratios ob- 
tained in the vicinity of mid-rapidity as invited, at RHIC energy, by the limited 
acceptance of the STAR and PHENIX experiments. The latter approach is 
appropriate, clearly, in the limit of boost-invariant rapidity distributions where 
hadron production ratios would not depend on the choice of the observational 
rapidity interval. We have shown in section 12.21 that such conditions do, in 
fact, set in at top RHIC energy, as referred to in Fig. I3.2i However, at low -^i 
the y-distributions are far from boost-invariant, and the total rapidity gap Ay 
may become comparable, in the extreme case, to the natural rapidity widths 
of hadrons emitted in the idealized situation of a single, isotropically decaying 
fireball positioned at mid-rapidity. Its rapidity spectra, equation 12.51 resemble 
Gaussians with widths Fj ~ 2.35 (T/mj)^'^ for hadron masses m,. Clearly, 
the particle ratios {dNi/dy)/{dNj/dy) then depend strongly on the position 
of the rapidity interval dy: away from y = heavy hadrons will be strongly 
suppressed, and particle yields in narrow rapidity intervals are useless for a 
statistical model analysis unless it is known a priori that the radiator is a single 
stationary spherical fireball [llOj . This is not the case toward top SPS energy 
(see Fig. l2.9p . due to significant primordial longitudinal expansion of the hadron 
emitting source. Given such conditions, the total multiplicity per collision event 
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(the invariant yield divided by the total overall inelastic cross section) should 
be employed in the SHM analysis, as is exemplified in Fig. I3.1[ 

3.1 Hadronic freeze-out from expansion flow 

The hadronic multiplicities result from integration of the invariant triple differ- 
ential cross section over pT and y. Instrumental, experiment-specific conditions 
tend to result in incomplete pT and/or y acceptances. It is important to ascer- 
tain that the effects of hydrodynamic transverse and longitudinal flow do not 
blast a significant part of the total hadron yield to outside the acceptance, and 
that they, more generally, do not change the relative hadron yield composition, 
thus basically affecting the SHM analysis. To see that hadronization incorpo- 
rates only the internal energy in the co- moving frame [llUj . we first assume 
that hadrochemical freeze-out occurs on a sharp hypersurface S, and write the 
total yield of particle species i as 

Ni = 1^1 p^dV^(x) Mx,p) = I (iV^(x)if (x) (3.1) 

where d^a is the outward normal vector on the surface, and 

Jiix) =9^ [ tfVe(/)(5(/ _ mf) pf^^exp [p • u{x) - /i,]/r ± 1)-^ (3.2) 



is the grand canonical number current density of species i, fii the chemical 
potential, u{x) the local flow velocity, and gi the degeneracy factor. In thermal 
equilibrium it is given by 

jf(x) = pi{x)u^{x) with 

Pi{x) = u^(x)jf (x) = / (fp2e{p°)5{p'^ - mj) p ■ u(x) fi{p ■ u{x); T; pi) 



= fd^p'fi{Epr,T,fii)=pi{T,p,). (3.3) 

Here Ep' is the energy in the local rest frame at point x. The total particle 
yield of species i is therefore 

N, = p,{T, Pi) [ (fa^{x)u^'{x) = pi{T, p,) ^^(u'^) (3.4) 

where only the total comoving volume Vs of the freeze-out hypersurface S 
depends on the flow profile u^. V \s thus a common total volume factor at 
hadronization (to be determined separately), and the flow pattern drops out 
from the yield distribution over species in 47r acceptance |110 j. For nuclear 
collisions at SPS energies and below one thus should perform a SHM analysis 
of the total, 47r-integrated hadronic multiplicities, as was done in Fig. 13.11 



We note that the derivation above illustrates the termination problem of the 
hydrodynamic description of A-|-A collisions, the validity of which depends on 
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conditions of a short mean free path, A < 1 fm. A precise argumentation sug- 
gests that two different free paths are relevant here, concerning hadron occu- 
pation number and hadron spectral freeze-out, respectively. As hadrochemical 
freeze-out occurs in the immediate vicinity of Tc (and Th ~ 160 — IGSMeT^ from 
Figs. 13. H 13. 2p . the hadron species distribution stays constant throughout the 
ensuing hadronic phase, i.e the "chemical" mean free path abruptly becomes 
infinite at Th, whereas elastic and resonant rescattering may well extend far 
into the hadronic phase, and so does collective pressure and flow. In fact we 
have seen in section 12.61 that the decoupling from flow occurs at Tp as low as 
90 — 100 MeV (Fig. 12.230 . Thus the hydrodynamic evolution of high ^/s colli- 
sions has to be, somehow artificially, stopped at the parton-hadron boundary in 
order to get the correct hadron multiplicities Ni, of equations 13.11 to 13.41 which 
then stay frozen-out during the subsequent hadronic expansion. 

The equations (|3.m3.4p demonstrate the application of the Cooper-Frye pre- 
scription [lllj for termination of the hydrodynamic evolution. The hyper- 
surface S describes the space-time location at which individual flow cells arrive 
at the freeze-out conditions, e = Cc and T = Tc, of hadronization. At this point, 
the resulting hadron/resonance spectra (for species i) are then given by the 
Cooper-Prye formula 

^^ = 1 ^ = T^T^s I fi{p-u{x),x)p-(fa{x), (3.5) 

where p^ fid^cr^j, is the local flux of particle i with momentum p through the 
surface S. For the phase space distribution / in this formula one takes the 
local equilibrium distribution at hadronic species freeze-out from the grand 
canonical SHM 

fi{E, x) = [exp{{E, - ^li{x))/T] ± 1]^^ (3.6) 

boosted with the local flow velocity u^{x) to the global reference frame by the 
substitution E ^>- p ■ u{x). Fixing T = T^ (taken e.g. from lattice QCD) the 
hadron multiplicities Ni then follow from equation \'6A\ and one compares to 
experiment, as in Figs. 13. H 13.21 In order now to follow the further evolution, 
throughout the hadronic rescattering phase, and to finally compare predictions 
of equation 13.51 to the observed flow data as represented by the various Fourier- 
terms of equation 12.201 one has to re-initialize (with hadronic EOS) the expan- 
sion from Ti{Tc = 165 MeV) until final decoupling [SH], at T ~ 100 MeV , thus 
describing e.g. radial and elliptic flow. 

Alternatively, one might end the hydrodynamic description at T = Tc and 
match the thus obtained phase space distribution of equation 13.51 to a micro- 
scopic hadron transport model of the hadronic expansion phase |95 tlll2| . This 
procedure is illustrated in Fig. 13.31 by an UrQMD [113J calculation of Bass and 
Dumitru [114j for central Au+Au collisions at top RHIC energy. We select here 
the results concerning the survival of the hadronic multiplicities Ni throughout 
the dynamics of the hadronic expansion phase, which we have postulated above, 
based on the equality of the hadronization temperatures, Th ~ 160 MeV , ob- 
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served in e~^e~ annihilation (Fig. l2.16|) . where no hadronic expansion phase ex- 
ists, and in central collisions of A ~ 200 nuclei (Figs. 13. 1113. 2p . In fact, Fig. 13.31 
shows that the {Ni} observed at the end of the hadronic cascade evolution 
agree, closely, with the initial {N^} as derived from a Cooper-Frye procedure 
(equation I3.4|) directly at hadronization. On the other hand, pT spectra and 
radial flow observables change, drastically, during the hadronic cascade expan- 
sion phase. 



Au+Auat RHIC 



u 

>> 



-a 



10^ 


^ 








■ Hydro -|- had. decays 


5 




■• 


■• 




• Hydro + UrQMD 


2 

10^ 






2 










10 
5 










2 
1 








■• % 


5 










2 

-in 1 


S.A. Bass et al. PRC60:021902 (1999) 




7r 


K 


K 


p a-p A a-A S" S+ iT a* 

hadron-spezies 



Figure 3.3: Modification of mid-rapidity hadron multiplicities in central Au-|-Au 
collisions at y/s = 200 GeV after chemical freeze-out at T = Tc- 
Squares show a hydrodynamic model prediction at T = Tc (without 
further interaction); circles show the result of an attached UrQMD 
hadronic cascade expansion calculation [ll4 J. 



The hadronic multiplicity distribution {Ni}, arising from the hadronization 
process at high ^/s, freezes-out instantaneously also in A+A collisions, and is 
thus preserved throughout the (isentropic) hadronic expansion phase. It is thus 
directly measurable and, moreover, its hadrochemical equilibrium features lend 
themselves to an analysis within the framework of Hagedorn-type statistical, 
grand canonical models. As we shall show below, the outcome of this analysis 
is contained in a [Thi^ib] parameter pair that reflects the conditions of QCD 
matter prevailing at hadronization, at each considered ^/s. In fact, the [T, /i] 
points resulting from the SHM analysis exhibited in Figs. 13. li 13.21 (at y/s = 17.3 
and 200 Ge^, respectively) have been shown in the QCD matter phase diagram 
of Fig. 11.11 to approach, closely, the parton-hadron phase coexistence line pre- 
dicted by lattice QCD. Thus, Th ~ Tc at high ^/s: hadrochemical freeze-out 
occurs in the immediate vicinity of QCD hadronization, thus providing for a 
location of the QCD phase boundary. 
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3.2 Grand canonical strangeness enhancement 

3.2 Grand canonical strangeness enhancement 

The statistical model analysis [HI 11071 I108J of the hadronization species distri- 
bution Ni in A+A collisions is based on the grand canonical partition function 
for species i, 

1 7 =^j^ r ^^ n 7^ 

" * evr^T io E,[k)exp{{E,{k)-iJi,)/T}±l ^ ' ^ 

where Ef = k"^ + m?, and /Xj = hbBi + /x^S'^ + /i/I| is the total chemical 
potential for baryon number S, strangeness S and isospin 3-component I3. Its 
role in equation 13.71 is to enforce, on average over the entire hadron source 
volume, the conservation of these quantum numbers. In fact, making use of 
overall strangeness neutrality (^^ NiSi = 0) as well as of conserved baryon 
number (participant Z+N) and isospin (participant (N-Z)/Z) one can reduce 
/Xj to a single effective potential /i^. Hadronic freeze-out is thus captured in 
three parameters, T, V and /i^. The density of hadron/resonance species i then 
results as 

rij = — — InZi (3-8) 



V5^ 



which gives 



''*- ^"^- (2vr)2 io exp{{E,{k)-fii)/T}±l- ^^'^^ 



We see that the common freeze-out volume parameter is canceled if one consid- 
ers hadron multiplicity ratios, Ni/Nj, as was done in Fig. 13.21 Integration over 
momentum yields the one-particle function 



yTgi 2 V^ (±1)"+^ fnmi\ /nfii 
m- > Kn PT.n 



n=l 

where K2 is the modified Bessel function. At high T the effects of Bose or Fermi 
statistics (represented by the ±1 term in the denominators of equation 13.71 and 
equation 13. 9p may be ignored , finally leading to the Boltzmann approximation 

' = ^^ "^i ^2 [-y) exp \^—j (3.11) 

which is the first term of equation I3.1U[ This approximation is employed 
throughout the SHM analysis. It describes the primary yield of hadron species 
i, directly at hadronization. The abundance of hadronic resonance states is 
obtained convoluting equation 13.91 with a relativistic Breit-Wigner distribution 
[T9j . Finally, the overall multiplicity, to be compared to the data, is determined 
as the sum of the primary multiplicity equation 13.111 and the contributions 
arising from the unresolved decay of heavier hadrons and resonances: 



57 



3 Hadronization and hadronic freeze-out in A+A collisions 

After having exposed the formal gear of grand canonical ensemble analysis 
we note that equation 13.111 permits a simple, first orientation concerning the 
relation of T to /ib in A+A collisions by considering, e.g., the antiproton to 
proton production ratio. From equation 13.111 we infer the simple expression 



N{p)/N{p) = exp{-2fiB/T). 



(3.13) 



Taking the mid-rapidity value 0.8 for p/p (from Fig. 13. 2p at top RHIC energy, 
and assuming that hadronization occurs directly at the QCD phase boundary, 
and hence T ^ Tc ^ 165 MeV, we get /u^ — 18 MeV from equation 13.131 in 
close agreement with the result, ^b = 20 MeV, obtained [108J from the full 
SHM analysis. Equation 13.131 illustrates the role played by ^b in the grand 
canonical ensemble. It logarithmically depends on the ratio of newly created 
quark-antiquark pairs (the latter represented by the p yield), to the total num- 
ber of quarks including the net baryon number-carrying valence quarks (repre- 
sented by the p yield) . 
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Figure 3.4: The ratio of total K~^ to total vr"*" multiplicity as a function of y/s, 
in central Au+Au and Pb+Pb collisions and in p+p minimum bias 
collisions |100j . 



The most outstanding property of the hadronic multiplicities observed in cen- 
tral A+A collisions is the enhancement of all strange hadron species, by factors 
ranging from about 2 to 20, as compared to the corresponding production rates 
in elementary hadron-hadron (and e"'"e~ annihilation) reactions at the same 
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y/s. I.e. the nuclear collision modifies the relative strangeness output by a 
"nuclear modification factor", R^^ = N^^ /{).b Npart • -A^f^, which depends on 
^/s and Npart and features a hierarchy with regard to the strangeness number 
s = 1, 2, 3 of the considered species, Rf=i < R'^=2 < ^3^=3- These properties are 



illustrated in Figs. [T^ and 13.51 The former shows the ratio of total K~^ to pos- 
itive pion multiplicities in central Au+Au/Pb+Pb collisions, from lower AGS 
to top RHIC energies, in comparison to corresponding ratios from minimum 
bias p+p collisions [lOOj . We have chosen this ratio, instead of (K^) /Npart, 
because it reflects, rather directly, the "Wroblewski ratio" of produced strange 
to non-strange quarks [107], contained in the produced hadrons. 



2((s> + (^» 



(u) + (d) + (n) + (d) 



0.2 in pp 
0.45 in AA. 



(3.14) 



The low value of A^ in pp (and all other elementary) collisions reflects a quark 
population far away from u, d, s flavor equilibrium, indicating strangeness sup- 
pression [109J. 
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Figure 3.5: The nuclear modification factors Rg^i 2 3 for hyperon and anti- 
hyperon production in nucleus-nucleus collisions at ^/s = 17.3GeV, 
relative to the p+p reference at the same energy scaled by N]y{= 
Npart)- The NA49 data refer to total 4tt yields ^TW\ . Closed hues 
represent the inverse strangeness suppression factors from ref . [119j , 
at this energy. 



The so-called strangeness enhancement property of A+A collisions (obvious 
from Figs. l3.4l and l3.5p is, thus, seen as the removal of strangeness suppression; it 
is also referred to as a strangeness saturation, in SHM analysis [1071 1108| , for the 
reason that A^ ~ 0.45 corresponds to the grand canonical limit of strangeness 
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production, implicit in the analysis illustrated in Figs. l3.l1 and l3.2l The average 
Rf^i at -y/i > 10 GeV thus is about 2.2, both in the data of Fig. [33] and in the 
statistical model. It increases (Fig. reffig:Figure29) toward about 10 in s = 3 
production of fl hyperons. 

In order to provide for a first guidance concerning the above facts and ter- 
minology regarding strangeness production we propose an extremely naive ar- 
gument, based on the empirical fact of a universal hadronization temperature 
(Figs. \Tm EH E2D at high ^. Noting that (s) = {s) and (u) « («) ^ 
(d) ss (^d'j in a QGP system at fi^ near zero, and T = 165 MeV, just prior to 
hadronization. As ss {s) / (u) ss exp{{mu — 'ms)/T} = 0.45 at pT ^ if we 
take current quark masses, nis — rriu ~ 130 MeV . I.e. the value of As in A+A 
collisions at high ^/s resembles that of a grand canonical QGP at /i^ — > 0, as 
was indeed shown in a 3 flavor lattice QCD calculation |115j at T ~ Tc- On 
the contrary, a p+p collision features no QGP but a small fireball volume, at 
T K, Tc, within which local strangeness neutrality, (s) = (s) has to be strictly 
enforced, implying a canonical treatment |109j . In our naive model the expo- 
nential penalty factor thus contains twice the strangeness quark mass in the 
exponent. As in pp collisions ~ exp{2{mu — ms)/T} ~ 0.2, in agreement with 
the observations concerning strangeness suppression, which are thus referred to 
as canonical suppression. In a further extension of our toy model, now ignoring 
the u, d masses in comparison to rus ~ 135 MeV , we can estimate the hierarchy 
of hyperon enhancement in A+A collisions, 

Rj^ ex N^'^/NPP-Q.f, Npart ~ exp{{-sms + 2sms)/T} = 2.2, 5.1, 11.6 (3.15) 

for s = 1,2,3, respectively. Fig. 13.51 shows that these estimates correspond well 
with the data |116j for R derived in Att acceptance for A, Handfi as well as for 
their antiparticles, from central Pb+Pb collisions at ^/s = 17.3 GeV. The p+p 
reference data, and C+C, Si+Si central collisions (obtained by fragmentation 
of the SPS Pb beam) refer to separate NA49 measurements at the same energy. 

The above, qualitative considerations suggest that the relative strangeness yields 
reflect a transition concerning the fireball volume (that is formed in the course 
of a preceding dynamical evolution) once it enters hadronization. Within the 
small volumes, featured by elementary collisions (see section [3^ . phase space 
is severely reduced by the requirement of local quantum number conservation 
[1091 I117J including, in particular, local strangeness neutrality. These con- 
straints are seen to be removed in A+A collisions, in which extended volumes 
of high primordial energy density are formed. Entering the hadronization stage, 
after an evolution of expansive cooling, these extended volumes will decay to 
hadrons under conditions of global quantum mechanical coherence, resulting in 
quantum number conservation occuring, non-locally, and on average over the 
entire decaying volume. This large coherent volume decay mode removes the 
restrictions, implied by local quantum number balancing. In the above, naive 
model we have thus assumed that the hadronization of an Omega hyperon in 
A+A collisions faces the phase space penalty factor of only three s quarks to 
be gathered, the corresponding three s quarks being taken care of elsewhere in 
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the extended volume by global strangeness conservation. In the framework of 
the SHM this situation is represented by the grand canonical ensemble (equa- 
tions 13.91 I3.1ip : the global chemical potential ^b expresses quantum number 
conservation on average. Strict, local conservation is represented by the canon- 
ical ensemble. 

The grand canonical (GC) situation can be shown to be the large collision 
volume limit (with high multiplicities {Ni}) of the canonical (C) formulation 
[list I119J . with a continuous transition concerning the degree of canonical 
strangeness suppression |119j . To see this one starts from a system that is 
already in the GC limit with respect to baryon number and charge conserva- 
tion whereas strangeness is treated canonically. Restricting to s = 1 and — 1 
the GC strange particle densities can be written (from equation 13. lip as 

-fiLi = %±i A±i (3.16) 

with 

Zs=±i = -^ mj, K2 (^) expi^{BsfiB + QsA^q)/t| (3.17) 

and a "fugacity factor" X^ = exp {^is/T). The canonical strange particle 
density can be written as [119) 

nf = nf'=^ • (A,) (3.18) 

with an effective fugacity factor 

A, = -^ Mil (3.19) 

y/blb-l Io{x) 

where S±i = Yls=±i ^s=±i is the sum over all created hadrons and resonances 
with s = ±1, the In{x) are modified Bessel functions, and x = 2^/SlS-l is 
proportional to the total fireball volume V. In the limit j; ~ 1/ — > 00 the 
suppression factor /i(x)/Io(x) -^ 1, and the ratio S±l/^/SlS-l corresponds 
exactly to the fugacity As in the GC formulation (see equation I3.16p . Thus 
the C and GC formulations are equivalent in this limit, and the canonical 
strangeness suppression effect disappears. Upon generalization to the complete 
strange hadron spectrum, with s = ibl,ib2,ib3, the strangeness suppression 
factor results [119] as 

r/(s) = /,(x)//o(x). (3.20) 

In particular for small x (volume), r]{s) — > (x/2)*, and one expects that the 
larger the strangeness content of the particle the smaller the suppression factor, 
and hence the larger the enhancement in going from elementary to central A+A 
collisions. This explains the hierarchy addressed in equation l3.15l and apparent 
from the data shown in Fig. 13.51 In fact, the curves shown in this figure represent 
the results obtained from equation I3.20[ for s = 1,2,3 hyperon production at 
^/s = 17.3 GeV |119j . They are seen to be in qualitative agreement with the 
data. However the scarcity of data, existing at top SPS energy for total hyperon 
yields, obtained in 47r acceptance (recall the arguments in section \3A\i both for 
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A+A and p+p collisions does not yet permit to cover the SHM strangeness 
saturation curves in detail, for s > 1. 

This saturation is seen in Fig. l3.5[ to set in already at modest system sizes, but 
sequentially so, for ascending hyperon strangeness. Note that SHM saturation 
is sequentially approached, from equation I3.2UI with increasing fireball volume 
V. In order to make contact to the experimental size scaling with centrality, e.g. 
Npart , the model of ref . |119j , which is illustrated in Fig. \'S.5\ has converted the 
genuine volume scale to the Npart scale by assuming a universal eigenvolume of 
7 fm^ per participant nucleon. I.e. Npart = 10 really means a coherent fireball 
volume of 70 fm^, in Fig. 13.51 Within this definition, saturation of s = 1,2,3 
sets in at fireball volumes at hadronization of about 60, 240 and 600 /m^, 
respectively: this is the real message of the SHM curves in Fig. 13.51 
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Figure 3.6: The mid-rapidity K^ to vr^ ratio vs. Npart in minimum bias Au+Au 
collisions at ^/s = 200GeV, compared to the percolation model |120j 
(solid line); a prediction of which for Cu+Cu at similar energy is 
given by the long dashed line (see text for detail). 



The above direct translation of coherent fireball volume to participant number is 
problematic [120j as it assumes that all participating nucleons enter into a single 
primordially coherent fireball. This is, however, not the case [120j particularly in 
the relative small scattering systems that cover the initial, steep increase of ??(s), 
where several local high density clusters are formed, each containing a fraction 
of Npart- This is revealed by a percolation model [12U) of cluster overlapp 
attached to a Glauber calculation of the collision/energy density. At each Npart 
an average cluster volume distribution results which can be transformed by 
equation 13.201 to an average {??(s,y)} distribution whose weighted mean is 
the appropriate effective canonical suppression factor corresponding to Npart- 
On the latter scale, the SHM suppression curve thus shifts to higher Npart, 
as is shown in Fig. 13.61 for the K^ /n^ ratio vs. Npart, measured at mid- 
rapidity by PHENIX in Au-|-Au collisions at y/s = 200GeV , which is reproduced 
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by the percolation model [120J. Also included is a prediction for Cu+Cu at 
this energy which rises more steeply on the common Npart scale because the 
collision and energy density reached in central Cu+Cu collisions, at Npart ~ 100, 
exceeds that in peripheral Au+Au collisions (at the same Npart) which share 
a more prominent contribution from the dilute surface regions of the nuclear 
density profile. We note, finally, that this deviation from universal Npart scaling 
does not contradict the observations of a perfect such scaling as far as overall 
charged particle multiplicity densities are concerned (recall Fig. 12. lip which are 
dominated by pions, not subject to size dependent canonical suppression. 



3.3 Origin of hadro-chemical equilibrium 

The statistical hadronization model (SHM) is not a model of the QCD con- 
finement process leading to hadrons, which occurs once the dynamical cooling 
evolution of the system arrives at Tc- At this stage the partonic reaction volume, 
small in elementary collisions but extended in A+A collisions, will decay (by 
whatever elementary QCD process) to on-shell hadrons and resonances. This 
coherent quantum mechanical decay results in a de-coherent quasi-classical, pri- 
mordial on-shell hadron-resonance population which, at the instant of its for- 
mation, lends itself to a quasi-classical Gibbs ensemble description. Its detailed 
modalities (canonical for small decaying systems, grand canonical for extended 
fireballs in A+A collisions) , and its derived parameters [T, /i^] merely recast 
the conditions, prevailing at hadronization. The success of SHM analysis thus 
implies that the QCD hadronization process ends in statistical equilibrium con- 
cerning the hadron-resonance species population. 

In order to identify mechanisms in QCD hadronization that introduce the hadro- 
chemical equilibrium we refer to jet hadronization in e~^e~ annihilation reac- 
tions, which we showed in Fig. 12. 16) to be well described by the canonical SHM. 
In di-jet formation at LEP energy, y^ = 92 GeV, we find a charged particle 
multiplicity of about 10 per jet, and we estimate that, likewise, about 10 pri- 
mordial partons participate on either side of the back-to-back di-jet [85]. There 
is thus no chance for either a partonic or hadronic, extensive rescattering toward 
chemical equilibrium. However, in the jet hadronization models developed by 
Amati and Veneziano |83], Webber p^l] and Ellis and Geiger [85] the period of 
QCD DGLAP parton shower evolution (and of perturbative QCD, in general) 
ends with local color neutralization, by formation of spatial partonic singlet 
clusters. This QCD "color pre-confinement" [83] process reminds of a coales- 
cence mechanism, in which the momenta and the initial virtual masses of the 
individual clustering partons get converted to internal, invariant virtual mass 
of color neutral, spatially extended objects. Their mass spectrum [121] extends 
from about 0.5 to 10 GeV. This cluster mass distribution, shown in Fig. 13. 7^ 
represents the first stochastic element in this hadronization model. 

The clusters are then re- interpreted within non-perturbative QCD: their inter- 
nal, initially perturbative QCD vacuum energy gets replaced by non-perturbative 
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Figure 3.7: Invariant mass spectrum of color neutralization clusters in the 
Veneziano- Webber hadronization model |83 t I121| . 



quark and gluon condensates, making the clusters appear like hadronic reso- 
nances. Their subsequent quantum mechanical decay to on-shell hadrons is 
governed by the phase space weights given by the hadron and resonance spec- 
trum [85lll21j . I.e. the clusters decay under "phase space dominance" [85], the 
outcome being a micro-canonical or a canonical hadron and resonance ensem- 
ble [HI 1107] . The apparent hadro-chemical equilibrium thus is the consequence 
of QCD color neutralization to clusters, and their quantum mechanical decay 
under local quantum number conservation and phase space weights. We note 
that the alternative description of hadronization, by string decay |122| . contains 
a quantum mechanical tunneling mechanism, leading to a similar phase space 
dominance [T23J . 

Hadronization in e^e^ annihilation thus occurs from local clusters (or strings), 
isolated in vacuum, of different mass but similar energy density corresponding 
to QCD confinement. These clusters are boosted with respect to each other 
but it was shown |124j that for a Lorentz invariant scalar, such as multiplicity, 
the contributions of each cluster (at similar T) can be represented by a single 
canonical system with volume equal to the sum of clusters. In the fit of Fig. 12.161 
this volume sum amounts to about 45 frn? [84]; the individual cluster volumes 
are thus quite small, of magnitude a few /m^ [85.J. This implies maximum 
canonical strangeness suppression but may, in fact, require a micro-canonical 
treatment of strangeness [109j . implying a further suppression. These MC ef- 
fects are oftentimes included [125] in the canonical partition functions by an 
extra strangeness fugacity parameter 7^ < 1 which suppresses s = 1,2,3 in a 
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hierarchical manner, {Ni{s)) ~ (isY*- The fit of Fig. [2.161 requires 7^ = 0.66, a 
value typical of canonical multiplicity analysis in p + p, p + p and e~^e~ annihi- 
lation collisions |109j at -y/s > 30 GeV. 

The above picture, of hadrochemical equilibrium resulting from the combined 
stochastic features of QCD color neutralization by cluster formation, and sub- 
sequent quantum mechanical decay to the on-shell hadron and resonance spec- 
trum (under phase space governance) lends itself to a straight forward extension 
to A+A collisions. The essential new features, of grand canonical hadronization 
including strangeness enhancement, should result from the fact that extended 
space-time volumes of e > €crit sue formed in the course of primordial partonic 
shower evolution, an overlap effect increasing both with ^/s and with the size of 
the primordial interaction volume. As the volume of the elementary hadroniza- 
tion clusters amounts to several fm^ it is inevitable that the clusters coalesce, 
to form extended " super-cluster" volumes prior to hadronization [120J . As these 
super-clusters develop toward hadronization via non perturbative QCD dynam- 
ics, it is plausible to assume an overall quantum mechanical coherence to arise 
over the entire extended volume, which will thus decay to hadrons under global 
quantum number conservation, the decay products thus modeled by the GC 
ensemble. 

Our expectation that space-time coalescence of individual hadronization clus- 
ters will lead to a global, quantum mechanically coherent extended super-cluster 
volume, that decays under phase space dominance, appears as an analogy to 
the dynamics and quantum mechanics governing low energy nuclear fission from 
a preceding "compound nucleus" |126j . Note that the observation of a smooth 
transition from canonical strangeness suppression to grand canonical satura- 
tion (Figs. 13.51 13. 6|) lends further support to the above picture of a percolative 
growth ^120j of the volume that is about to undergo hadronization. 

An extended, coherent quark gluon plasma state would, of course, represent 
an ideal example of such a volume [127J and, in fact, we could imagine that 
the spatial extension of the plasma state results from a percolative overlap 
of primordial zones of high energy density, which becomes more prominent 
with increasing y^ and Npart- A QGP state preceding hadronization will thus 
lead to all the observed features. However, to be precise: the hadronizing 
QCD system of extended matter decaying quantum coherently, could still be a 
non-equilibrium precursor of the ideal equilibrium QGP, because we have seen 
above that hadrochemical equilibrium also occurs in e^e~ annihilation, where 
no partonic equilibrium exists. It gets established in the course of hadronization, 
irrespective of the degree of equilibrium prevailing in the preceding partonic 
phase. 



3.4 Hadronization vs. rapidity and ^/s 

We have argued in section 13.11 that, at relatively low ^/s, the total rapidity 
gap Ay does not significantly exceed the natural thermal rapidity spreading 
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width Fj ss 2.35 {T/rriiY'^ of ^ single, isotropically decaying fireball, centered 
at mid-rapidity and emitting hadrons of mass nii |110j . However, this procedure 
involves an idealization because in the real Pb+Pb collision the intersecting di- 
lute surface sections of the nuclear density profiles will lead to a significant 
contribution of single-scattering NN collisions, outside the central high density 
fireball. The leading hadron properties of such "corona collisions" result in 
wider proper rapidity distributions, quite different from those of the central 
fireball decay hadrons. Their contribution will thus be prominent near tar- 
get/projectile rapidity, and will feature a canonically suppressed strangeness. 
The one-fireball assumption, although inevitable at small Ay, does not quite 
resemble the physical reality. This may explain the need for an additional 
strangeness suppression factor in the GC one-particle partition function (equa- 
tion 13. 7p that has, unfortunately, also been labeled 7^ but expresses physics 
reasons quite different from the extra suppression factor that reflects micro- 
canonical phase space constraints in elementary collisions. It turns out that all 
GC analysis of central A-l-A collisions at low \/s, and addressed to total 47r 
multiplicities, requires a 7^ of 0.7 - 0.85 |19j; in the fit of Fig. I3.1l ^c = 0.84. 

At RHIC, Ay w 11 ^ Fj, and such difficulties disappear: 7^ « 1 at midrapidity 
and, moreover, the wide gap permits a SHM analysis which is differential in 
y. Fig. 13.81 shows the y-dependence of the ratios ir" /tt^, K^ /K^ and p/p as 
obtained by BRAHMS [HH] in central Au+Au collisions at ^ = 200 GeV . 
The figure shows a dramatic dependence of the p/p ratio, which reflects the 
local baryochemical potential according to equation 13.131 At ycM > 1 the p/p 
ratio drops down steeply, to about 0.2 at y ~ 3.5, thus making close contact to 
the top SPS energy value obtained by NA49 [T^. The K-/K+ ratio follows 
a similar but weaker drop-off pattern, to about 0.65 again matching with the 
top SPS energy value of about 0.6 [130j . The deviation from unity of these 
ratios reflects the rapidity densities of initial valence u, d quarks, relative to the 
densities of newly created light and strange quark-antiquark pairs, i.e. the y 
distribution of the net baryon number density, and of the related baryo-chemical 
potential of the GC ensemble. Thus, in analyzing successive bins of the rapidity 
distributions in Fig. 13.81 the major variation in the GC fit concerns the baryo- 
chemical potential usiy) which increases from about 20 MeV (Fig. 13. 2p at 
mid-rapidity, to about 150 MeV at y > 3 while the hadronization temperature 
stays constant, at T = 160 MeV. This interplay between K~/K~^, p/p and fiB 
is illustrated |128| in the right hand panel of Fig. 13.81 and shown to be well 
accounted for by the GC statistical model [131]. 

These considerations imply that hadronization at RHIC (and LHC) energy 
occurs local in y-space and late in time. The density distribution of net baryon 
number results from the primordial pQCD shower evolution (c.f. section [2^ . 
and is thus fixed at formation time, to ^ 0.6 fm/c at RHIC. Hadronization of 
the bulk partonic matter occurs later, at t > 3 fm/c [86l[95], and transmits the 
local conditions in rapidity space by preserving the local net baryon quantum 
number density. Most importantly we conclude that hadronization occurs, not 
from a single longitudinally boosted fireball but from a succession of "super- 
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interpretation of the correlation between p/p and K^ /K^ in terms 
of baryo-chemical potential ^b variation in the grand canonical sta- 
tistical model. From |128j . 



clusters", of different partonic composition depending on y, and decaying at 
different time due to the Lorentz-boost that increases with y, in an "inside- 
outside" pattern (c.f. Fig. I2.17p . We are thus whitnessing at hadronization a 
Hubble expanding system of local fireballs. The detailed implications of this 
picture have not been analyzed yet. Note that a central RHIC collision thus 
does not correspond to a single hadronization "point" in the [T, ^] plane of 
Fig. 11.11 but samples {T, /i} along the QCD parton-hadron coexistence line 

Throughout this chapter we have discussed hadronic freeze-out at high a/s only 
(top SPS to RHIC energy), because of the proximity of the chemical freeze- 
out parameters [T, |Uf,] to the QCD phase boundary from lattice QCD, which 
suggests an overall picture of hadronization, to occur directly from a partonic 
cluster or super-cluster. Our discussion of the GC statistical hadronization 
model has been explicitly or implicitly based on the assumption that hadronic 
freeze-out coincides with hadronization. However, the GC model has also been 
applied successfully to hadro-chemical freeze-out at ^/s down to a few GeV 
[19t I107t I108| where it is not expected that the dynamical evolution traverses 
the phase boundary at all, but grand canonical multiplicity distributions, and 
their characteristic strangeness enhancement pattern, are observed throughout. 
Toward lower y^, T decreases while /Uf, increases, as is shown in Fig. 13.91 which 
presents a compilation of all reported freeze-out parameters |108j . 

These points have also been included in the phase diagram of Fig. 11.11 which 
shows that they are gradually branching away from the phase separation bound- 
ary line that could recently be predicted by lattice QCD once new methods had 
been developed to extend the theory to finite ^b [lO]. At ^/s > 20 GeV we 
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see that 



edQCD) ^en- ecc 



(3.21) 



where ecc is the freeze-out density inferred from GC analysis [T9 | 1107^1108) . 

In turn, the GC hadronic freeze-out points drop below the lattice QCD coexis- 
tence line at lower ^/s, implying that chemical freeze-out now occurs within the 
hadronic expansion phase. This requires a model of freeze-out, now governed 
by the properties of a high density hadronic medium, upon expansive cooling 
and dilution. Holding on to the model of a quantum mechanical de-coherence 
decay to on-shell hadrons that we discussed in section 13.31 we argue that an 
initial, extended high density hadronic fireball, given sufficient life-time at T 
smaller, but not far below Tc, could also be seen as a quantum mechanically 
coherent super-cluster, as governed by effective mean fields [133] . In such a 
medium hadrons, at T near Tc, acquire effective masses and/or decay widths 
far off their corresponding properties in vacuum: they are off-shell, approach- 
ing conditions of QCD chiral symmetry restoration as T —> Tc |134j . This 
symmetry is inherent in the elementary QCD Langrangian, and "softly" bro- 
ken within the light quark sector by the small non-zero current quark masses, 
but severely broken at T ^ by the high effective constituent quark masses 
that get dressed by non perturbative QCD vacuum condensates. Pictorially 
speaking, hadrons gradually loose this dressing as T —> Tc |135j . introducing a 
change, away from in vacuum properties, in the hadronic mass and width spec- 
trum. Such in-medium chiral restoration effects have, in fact, been observed in 
relativistic A-l-A collisions, by means of reconstructing the in-medium decay of 



68 



3.4 Hadronization vs. rapidity and y/s 

the p vector meson to an observed e"'"e^ pair [136] (see chapter [0]l . 

A dense, high T hadronic system, with mean-field induced off-shell constituents 
is also, clearly, quantum mechanically coherent. At a certain characteristic den- 
sity, e < Ec, and temperature T < Tc, as reached in the course of overall hadronic 
expansion, this extended medium will undergo a decoherence transition to clas- 
sical on-shell hadrons. Its frozen-out hadronic multiplicity distribution should 
be, again, characterized by the phase space weights of a grand canonical en- 
semble at r < Tc. Theoretical studies of such a mean field hadronic expansion 
mode |137j have also shown that such mechanisms play essentially no role at 
yfs > 20 GeV because the expanding system is already in rapid flow once it 
traverses the phase boundary, with an expansion time scale shorter than the 
formation time scale of mean field phenomena. At lower energies, on the other 
hand, the system might not even dive into the deconfined phase but spend 
a comparatively long time in its direct vicinity, at the turning point between 
compression and re-expansion where all dynamical time constants are large, and 
the hadron density is high, such that the inelastic hadronic transmutation rate 
becomes high (particularly in collisions of more than two hadronic reactants, 
with reaction rates ^138j proportional to e"), and sufficiently so for maintaining 
hadronic chemical equilibrium after it is first established at maximum hadron 
density, in low ^/s systems that do not cross the phase boundary at all. 

The GC freeze-out parameters [T,fi] at various ^/s in Fig. 13. 91 permit a smooth 
interpolation in the T, /i plane |139j . which, in turn, allows for GC model predic- 
tions which are continuous in ^/s. Such a curve is shown in Fig. 13.41 compared 
to the 4tt data points for the K^ /tt~^ multiplicity ratio in central collisions 
Au+Au/Pb+Pb, at all ^/s investigated thus far. It exhibits a smooth max- 
imum, due to the interplay of T saturation and /is fall-off to zero, but does 
not account for the sharp peak structure seen in the data at ^/s ~ 7 GeV and 
fiB ~ 480 MeV . This behavior is not a pecularity of the K^ channel only; 
it also is reflected in an unusually high Wroblewski ratio (see equation I3.14p 
obtained at y^ = 7.6 GeV, of As = 0.60 [19j. This sharp strangeness max- 
imum is unexplained as of yet. It implies that hadron formation at this ^/s 
reflects influences that are less prominent above and below, and most attempts 
to understand the effect |14H 11421 I143J are centered at the assumption that at 
this particular y/s the overall bulk dynamics will settle directly at the phase 
boundary where, moreover, finite /i^ lattice theory also expects a QCD critical 
point Uni [TT] . This would cause a softest point to occur in the equation of 
state, i.e. a minimum in the relation of expansion pressure vs. energy den- 
sity, slowing down the dynamical evolution |144t I145| , and thus increasing the 
sensitivity to expansion modes characteristic of a first order phase transition 
|143) . which occurs at /i^ > /ij;^**. Such conditions may modify the K/vr ratio 

(Fig.[33D [US]- 

It thus appears that the interval from top AGS to lower SPS energy, 5 < 
^/s < IQGeV , promises highly interesting information regarding the QCD phase 
diagram (Fig. ll.l| l in the direct vicinity of the parton-hadron coexistence line. In 
particular, the physics of a critical point of QCD matter deserves further study. 
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Observable consequences also comprise so-called "critical fluctuations" |146t 
I147J of multiplicity density, mean transverse momentum and hadron-chemical 
composition |148j . the latter in fact being observed near y/s = 7 GeV in an 
event by event study of the K/tt ratio in central Pb+Pb collisions [l49] . We 
shall return to critical point physics in chapter [71 



70 



4 Elliptic flow 



We have up to now mostly stressed the importance of central collisions and 
mid-rapidity data because they provide for the highest primordial energy den- 
sity and avoid problems caused by emission anisotropy and the presence of 
cold spectator sections of target and projectile nuclei. On the other hand, a 
fundamentally new window of observation is opened by non-central collisions 
as the finite impact parameter breaks cylinder symmetry, defining emission 
anisotropics with respect to the orientation of the impact vector b as we have 
shown in equation 12. 2Ui In a strongly interacting fireball collision dynamics, the 
initial geometric anisotropy of the reaction volume gets transferred to the final 
momentum spectra and thus becomes experimentally accessible. Furthermore, 
the high charged particle multiplicity allows for an event-by-event determina- 
tion of the reaction plane (direction of b) , enabling the study of observables at 
azimuth 99, relative to the known reaction plane. We shall show that this opens 
a window into the very early stages of A-l-A collisions onward from the end of 
nuclear interpenetration, at r ps 2 R{A)/^cM- Our observation thus begins at 
the extreme energy densities prevailing right at formation time (section [231 [23]), 
i.e. concurrent with the initialization phase of relativistic hydrodynamic expan- 
sion. We access the phase diagram of Fig. 11.11 in regions far above the QCD 
phase boundary. 

Before turning to the details of elliptic flow data we wish to illustrate j96j the 
above statements. Fig. 14.11 exhibits the transverse projection of primordial en- 
ergy density, assumed to be proportional to the number density of participant 
nucleons in the overlap volume arising from a Au-|-Au collision at impact pa- 
rameter b = 7 fm. The nuclear density profiles (assumed to be of Woods-Saxon 
type) intersect in an ellipsoidal fireball, with minor axis along the direction 
of b which is positioned at y = 0. The obvious geometrical deformation can 
be quantified by the spatial excentricity (unfortunately also labeled e in the 
literature) 

-(" = wly '"' 

where the averages are taken with respect to the transverse density profiles of 
Fig. 14.11 ex is zero for 6 = 0, reaching a value of about 0.3 in the case b = 7 fm 
illustrated in Fig. 14.11 

Translated into the initialization of the hydrodynamic expansion the density 
anisotropy implies a corresponding pressure anisotropy. The pressure is higher 
in X than in y direction, and thus is the initial acceleration, leading to an 
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Figure 4.1: Transverse projection of primordial binary collision density in an 
Au+Au collision at impact parameter 7 fm, exhibiting participant 
parton spatial excentricity 



increasing momentum anisotropy, 

fdxdy{T^^-Tyy) 

''^ ' Jdxdy{T-- + Tyy) ^ ' 

where Tf^ is the fluid's energy-momentum tensor. Fig. 14.21 shows |96t 1150) 
the time evolution of the spatial and momentum anisotropies for the collision 
considered in Fig. 14.11 implementing two different equations of state which are 
modeled with (without) implication of a first order phase transition in "RHIC" 
("EOSl"). A steep initial rise is observed for ep, in both cases: momentum 
anisotropy builds up during the early partonic phase at RHIC, while the spa- 
tial deformation disappears. I.e. the initial source geometry, which is washed 
out later on, imprints a flow anisotropy which is preserved, and observable as 
"elliptic flow". A flrst order phase transition essentially stalls the buildup of ep 
at about r = 3 fm/c when the system enters the mixed phase, such that the 
emerging signal is almost entirely due to the partonic phase. We have to note 
here that the "ideal fluid" (zero viscosity) hydrodynamics [15U| employed in 
Fig. 14.21 is. at flrst, a mere hypothesis, in view also of the fact that microscopic 
transport models have predicted a signiflcant viscosity, both for a perturbative 
QCD parton gas [92l [TS^lj and for a hadron gas [95l I152J . Proven to be correct 
by the data, the ideal fluid description of elliptic flow tells us that the QGP is 
a non-perturbative liquid [1531 1154] . 

Elliptic flow is quantifled by the coefficient V2 of the second harmonic term in the 
Fourier expansion (see equation I2.20p of the invariant cross section; it depends 
on ^/s, b, y and pT- Fig. 14.31 shows the ^/s dependence of V2 at mid-rapidity and 
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Figure 4.2: Time evolution of the spatial excentricity ex and the momentum 
space anisotropy Ep (equations 14.11 and 14. 2p in the hydrodynamic 
model of an Au+Au collision at 6 = 7/m, occuring at ^/s = 200Gey 
|96j . The dynamics is illustrated with two equations of state. 
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Figure 4.3: Energy dependence of the elliptic flow parameter V2 at mid-rapidity 
and averaged over pT, in Au+Au and Pb+Pb semi-peripheral colli- 
sions |155j . 



averaged over pT, in Au+Au/Pb+Pb semi-peripheral collisions [^ I155J . We 
see that the momentum space anisotropy is relatively small overall, but exhibits 
a steep rise toward top RHIC energy. 

Fig. l4.4l shows the (pseudo)-rapidity dependence of U2 at ^/s = 130 and 200 GeV 
as obtained by PHOBOS [156J for charged particles in minimum bias Au+Au 
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Figure 4.4: Pseudo-rapidity dependence of the pT-averaged elliptic flow coeffi- 
cient V2 for charged hadrons at ^/s = 130 and 200 GeV [156J. 



collisions. It resembles the corresponding charged particle rapidity density dis- 
tribution of Fig. \'2.7\ suggesting that prominent elliptic ffow arises only at the 
highest attainable primordial energy density. 
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Figure 4.5: Transverse momentum dependence of elliptic flow V2 for mesons and 
baryons in Au-|-Au collisions at \fs = 200 GeV . The hydrodynamic 
model [Ml I159J describes the mass dependence at pr < 2 GeV jc. 



That such conditions are just reached at top RHIC energy is shown in Figs. 14.51 
and WR The former combines STAR [ST] and PHENIX [EH] data for the 
Pt dependence of elliptic flow, observed for various identifled hadron species 
TT, i^,p, K^ and A, A in Au-|-Au at 200 GeV. The predicted hydrodynamic flow 
pattern [Ml I159J agrees well with observations in the bulk px < 2 GeV/c do- 
main. Fig. 14.61 (from [155| ) unifles average V2 data from AGS to top RHIC 
energies in a scaled representation [93] where V2 divided by the initial spatial 
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anisotropy tx is plotted versus charged particle mid-rapidity density per unit 
transverse area 5, the latter giving the density weighted transverse surface area 
of primordial overlap, Fig. 14.11 Fig. 14.61 includes the hydrodynamic predictions 
|951 [96l 11501 11591 I161J for various primordial participant or energy densities as 
implied by the quantity {1/ S)dnch/dy [93j. Scaling V2 by ex enhances the ellip- 
tic flow effect of near-central collisions where ex is small, and we see that only 
such collisions at top RHIC energy reach the hydrodynamical ideal flow limit in 
models that include an EOS ansatz which incorporates [96] the effect of a first 
order phase transition, which reduces the primordial flow signal as was shown 
in Fig. 1121 
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Figure 4.6: Elliptic flow V2 scaled by spatial excentricity e as a function of 
charged particle density per unit transverse area S, from AGS to 
top RHIC energy. The hydrodynamic limit is only attained at RHIC 

^551. 



At top RHIC energy, the interval between tg ~ 0.6 fm/c^ and hadronization 
time, tn ^ ^ ffn/c^ is long enough to establish dynamical consequences of an 
early approach toward local equilibrium. The "lucky coincidence" of such a pri- 
mordial resolution of dynamical time scale, with the extreme primordial density, 
offered by semi-central collisions of heavy nuclei, results in an extremely short 
mean free path of the primordial matter constituents, thus inviting a hydro- 
dynamic description of the expansive evolution. Consistent application of this 
model reveals a low viscosity: the primordial matter resembles an ideal fluid, 
quite different from earlier concepts, of a weakly interacting partonic gas plasma 
state (QGP) governed by perturbative QCD screening conditions [MlHI]- 

A further, characteristic scaling property of elliptic flow is derived from the pT 
dependence of V2, observed for the different hadronic species. In Fig. 14.51 one 
observes a hadron mass dependence, the V2 signal of pions and charged kaons 
rising faster with px than that of baryons. Clearly, within a hydrodynamic flow 
velocity field entering hadronization, heavier hadronic species will capture a 
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higher pT, at a given flow velocity. However, unlike in hadronic radial expansion 
flow phenomena (c.f. section 12. 6p it is not the hadronic mass that sets the scale 
for the total px derived, per particle species, from the elliptic flow field, but the 
hadronic valence quark content. This conclusion is elaborated fMj in Fig. 14. 7i 
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Figure 4.7: V2 vs. px (left panel) and transverse kinetic energy KEt = mT — mo 
(middle) for several hadronic species in min. bias Au+Au collisions 
at ^/s = 200 GeV, showing separate meson and baryon branches. 
Scaling (right panel) is obtained by valence quark number n^, di- 
viding V2 and KEt 



The left panel shows measurements of the pT dependence of V2 for several 
hadronic species, in minimum bias Au+Au collisions at ^/s = 200 GeV |161j . 
The middle panel bears out the hydrodynamically expected |162| particle mass 
scaling when V2 is plotted vs. the relativistic transverse kinetic energy KEt = 
rriT — 'm where niT = {p^ + rn?Yi'^ . For KEt > 1 GeV , clear splitting into a 
meson branch (lower V2) and a baryon branch (higher V2) occurs. However, both 
of these branches show good scaling separately. The right panel shows the result 
obtained after scaling both V2 and KEt (i-e. the data in the middle panel) by 
the constituent quark number, Uq = 2 for mesons and n^ = 3 for baryons. The 
resulting perfect, universal scaling is an indication of the inherent quark degrees 
of freedom in the flowing matter as it approaches hadronization. We thus assert 
that the bulk of the elliptic flow signal develops in the pre-hadronization phase. 

The above scaling analysis has been extended to $ meson and Q hyperon pro- 
duction, and also to first PHENIX results |163j concerning elliptic flow of the 
charmed D meson ^94J, with perfect agreement to the observations made in 
Fig. 14.71 of a separation into meson/hadron branches on the KEt scale, which 
merge into a universal V2 scaling once both V2 and KEt per valence quark are 
considered. The observation that the D meson charm quark apparently shares 
in the universal flow pattern is remarkable as its proper relaxation time is, in 
principle, lengthened by a factor M/T |164| . A high partonic rescattering cross 
section a is thus required in the primordial QGP fireball, to reduce the par- 
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tonic mean free path A = Xjuo (where n is the partonic density), such that 
A <SC A^''^ (the overall system size) and, simultaneously, A < 1 /m in order 
to conform with the near-zero mean free path implication of the hydrodynamic 
description of the elliptic flow, which reproduces the data gathered at RHIC 
energy. The presence of a high partonic rescattering cross section was born out 
in a parton transport model study [92] of the steep linear rise of the elliptic 
flow signal with p^ (Fig. 14. 5p . In such a classical Boltzmann equation approach 
the cross sections required for the system to translate the initial spatial into 
momentum space anisotropy, fast enough before the initial deformation gets 
washed out during expansion (c.f. Fig. 14. 2p . turn out to exceed by almost an 
order of magnitude the values expected in a perturbative QCD quark-gluon gas 



The non-perturbative quark-gluon plasma is thus a strongly coupled state (which 
has been labeled sQGP [165J ). At RHIC energy this reduces the partonic mean 
free path to a degree that makes hydrodynamics applicable. A Navier-Stokes 
analysis |166| of RHIC flow data indicates that the viscosity of the QGP must 
be about ten times smaller than expected if the QGP were a weakly interacting 
pQCD Debye screened plasma. This justifies the use of perfect fluid dynamics. 

Considering first attempts to derive a quantitative estimate of the dimensionless 
ratio of (shear) viscosity to entropy,we note that rj/s is a good way to charac- 
terize the intrinsic ability of a substance to relax toward equilibrium |167j . It 
can be estimated from the expression 



rj/s ^ TXf cs (4.3) 

where T is the temperature, Xf the mean free path, and Cs is the sound speed 
derived from the partonic matter EOS. A fit by the perfect fluid "Buda-Lund" 
model [IBS] to the scaled V2 data shown in Fig. 14.71 vields T = 165 ± 3 MeV; 
Cs is estimated as 0.35 ± 0.05 [1621 IM] . and Xj ^ 0.30 fm taken from a parton 
cascade calculation including 2 <-> 3 scattering |169] . The overall result is 



rj/s = 0.09 ± 0.02 (4.4) 

in agreement with former estimates of Teaney and Gavin [170j . This value 
is very small and, in fact, close to the universal lower bound of rj/s = 1/Att 
recently derived field theoretically [l7l]. 

Elliptic flow measurements thus confirm that the quark-gluon matter produced 
as y/s -^ 200 GeV is to a good approximation in local thermal equilibrium 
up to about 3 — 4 fm/c. In addition, the final hadron mass dependence of 
the flow pattern is consistent with a universal scaling appropriate for a nearly 
non- viscous hydrodynamic flow of partons, and the observed V2 signal reflects a 
primordial equation of state that is consistent with first numerical QCD compu- 
tations [153^1154] 1165] of a strongly coupled quark-gluon plasma (sQGP) state. 
First estimates of its proper shear viscosity to entropy ratio, ry/s, are emerg- 
ing from systematic analysis of the elliptic flow signal. At lower y/s precursor 
elliptic flow phenomena are observed, as well, but are more difficult to ana- 
lyze as the crucial, new feature offered by top RHIC energies is missing here: 
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a clear cut separation in time, between primordial formation of local partonic 
equilibrium conditions, and hadronization. At RHIC (and at future LHC) en- 
ergy elliptic flow systematics thus captures the emerging quark-gluon state of 
QCD at energy densities in the vicinity of e = 6 — 15 GeV/ fnfi , at temperature 
T « 300 MeV , and ^b -^ 0, describing it as a strongly coupled, low viscosity 
liquid medium. 
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5 In-medium attenuation of high px 
hadron and jet production 



In the preceding sections we have fohowed the dynamical evolution of bulk 
matter in A+A collisions (at pr ^ 2 GeV which covers about 95% of the 
hadronic output), from initial formation time of partonic matter which reflects 
in charged particle transverse energy and multiplicity density, also giving birth 
to hadrons, and to the elliptic expansion flow signal. 

An alternative approach toward QCD plasma diagnostics exploits the idea 
[411 I172J of implanting partonic products of primordial high Q^ processes into 
the evolving bulk medium, that could serve as "tracers" of the surrounding, 
co-traveling matter. The ideal situation, of being able to scatter well defined 
partons, or electrons, from a plasma fireball, is approximated by employing 
primordially formed charm-anticharm quark pairs [H] , or leading partons from 
primordial di-jet production |172j . Both processes are firmly anchored in per- 
turbative QCD and well studied in elementary collisions (where such partons 
are directly released into vacuum), which thus serve as a reference in an anal- 
ysis that quantifies the in-medium modification of such tracer partons. Not a 
surprise, in view of our above inferences, from elliptic flow, of a high tempera- 
ture, strongly coupled primordial medium: these in-medium modifications are 
quite dramatic, leading to a suppression of J/^ production from primordial cc 
pairs (chapter E]) , and to high px hadron and jet quenching, the subject of this 
chapter. 

5.1 High pt inclusive hadron production quenching 

At top RHIC energy, ^/s = 200 GeV, di-jet production from primordial hard 
pQCD parton-parton collisions (of partons from the initial baryonic structure 
functions) is the source of "leading" partons, with Et up to about 30 GeV. 
They are derived from the inclusive cross section arising if the A-l-A collision is 
considered, first, as an inco/iereni superposition of independent nucleon-nucleon 
collisions, as enveloped within the target-projectile nucleon densities. In this 
framework, the pQCD cross section for producing an Et parton in A+B takes 
the form of " factorization" 11731 



d(^ _S^ f ,.f J., r /„. ^i• /„ \ daab 



dXa I dXbfa/AiXa)fb/Bixb) ^f„^^, (5-1) 

a,b 



dErdy ff J^^ J^^ "' dExdy 



where the f{x) are the parton distributions inside projectile A and target B 
nuclei, and the last term is the pQCD hard scattering cross section. This 
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equation describes the primordial production rate of hard partons, leading to 
the conclusion [172] that, at RHIC energy, all hadrons at pj- > 6 — 10 GeV 
should arise from initial pQCD parton production. 

As partons are effectively frozen during the hard scattering, one can treat each 
nucleus as a collection of free partons. Thus, with regard to high pT production, 
the density of partons within the parton distribution function of an atomic 
number A nucleus should be equivalent to the superposition of A independent 
nucleons A'^: 

fa/A{^,Q')=Afa/N{x,Q^). (5.2) 

Prom equations 15.11 and 15.21 it is clear that the primordial high Q^ inclusive 
parton cross section in a minimum bias A+B reaction scales as A • B times 
the corresponding (N+N or) p+p cross section. Furthermore, as each leading 
parton ends up in an observed high pT hadron h , we thus write the invariant 
hard hadron cross section as 

EdaAB-.h/cfip = A-B-E dapp^h/d^P- (5.3) 

Since nucleus-nucleus experiments usually measure invariant yields Nh for a 
given centrality bin, corresponding to an average impact parameter b, one writes 
instead: 

EdNAB^h{b)/d^p = {TAB{b)) E dapp^^/d^P, (5.4) 

where r^s(6) is the Glauber geometrical overlap function of nuclei A, B at im- 
pact parameter b, which accounts for the average number of participant parton 
collisions at given impact geometry [174j . {Ncoii{b)). One can thus quantify the 
attenuating medium effects, as experienced by the primordially produced tracer 
parton on its way toward hadronization, by the so-called nuclear modification 
factor for hard collisions (analogous to equation 12.61 that refers to soft, bulk 
hadron production): 

D f h\ d^NAs/dydpT , . 

Rab{pt, y, b) = — — ^ , , , ■ (5.5) 

{TAB{b)) d'^fjpp/dydpT 

Obviously, this concept of assessing the in-medium modification of hadron pro- 
duction at high pt requires corresponding p+p collision data, as a reference 
basis. Such data have been, in fact, gathered at top RHIC, and top SPS ener- 
gies, ^/s = 200 and 17.3 GeV , respectively. Alternatively, in situations where 
the relevant reference data are not known, one considers the production ratio 
of hadronic species /i, observed in central relative to peripheral collisions: 

d^Nh{bi)/dydpt {TAB{b2)) ,_, 

^^^^^^'^^ = imMJd^t " -WI^(b^) ^'-'^ 

where 6i << 62 are the average impact parameters corresponding to the em- 
ployed trigger criteria for "central" and "peripheral" A+A collisions, respec- 
tively. This ratio recasts, to a certain extent, the in-medium attenuation anal- 
ysis, offered by Rab, insofar as peripheral collisions approach the limiting con- 
ditions, of a few single nucleon-nucleon collisions occuring in the dilute surface 
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Figure 5.1: The nuclear modification factor Raa for vr'' in min. bias Cu+Cu 
and Au+Au collisions at t/s = 200 GeV, in the range pT > 7GeV/c 
|175j . plotted vs. centrality as measured by participant nucleon 
number Npart- 



sections of the nuclear density profiles, i.e. essentially in vacuum. 

Employing the above analysis schemes, the RHIC experiments have, in fact, 
demonstrated a dramatic in-medium suppression of the high px yield, for vir- 
tually all hadronic species. Fig. 15.11 shows Raa for neutral pions produced in 
min. bias Cu-^Cu and Au-^Au colhsions at ^ = 200 GeV where PHENIX 
extended the pT range up to 18 GeV/c [175j : the nuclear modification factor 
refers to the range pT > 7.0 GeV/c and is shown as a function of centrality, 
reflected by Npart- We infer a drastic suppression, by an Raa ~ 0.25 in near 
central collisions. Fig. 15.21 shows the pT dependence [175) of neutral pion Raa 
in central Au-|-Au collisions {Npart = 350), with a suppression to below 0.2 at 
PT > 4 GeV/c. 

Note that Raa can not reach arbitrarily small values because of the unatten- 
uated contribution of quasi-in- vacuum surface " corona" single nucleon- nucleon 
collisions, closely resembling the p + p ^f n^ + X inclusive yield employed in 
the denominator of Raa , equation 15.51 Even in a central trigger configuration 
this suggests a lower bound, Raa ~ 0.15. Nuclear attenuation of high px pions 
thus appears to be almost "maximal" , implying a situation in which the interior 
sections of the high energy density fireball feature a very high opacity, i.e. they 
are almost "black" to the high pT partons contained in pions. 
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Figure 5.2: Raa for vr*^ production in central Au+Au collisions at ^/s = 200GeV 
|175j, compared to a hydrodynamic calculation for different opaci- 
ties (transport coefficients |195] ) of the plasma |204j . 



We remark here, briefly, on a confusing feature that occurs in Fig. 15.21 (and in 
the majority of other Raa vs. pt plots emerging from the RHIC experiments): 
the pQCD number of collision scaling employed in Raa does not describe pion 
production at low pT as we have demonstrated in chapter [2j The entries at 
Pt < 3 GeV/c are thus besides the point, and so are pQCD guided model 
predictions, as shown here [176j . 

Nuclear modification analysis at RHIC covers, by now, the high pT produc- 
tion of a multitude of mesonic and baryonic species |177j . most remarkably 
even including the charmed D meson which is measured via electrons from 
semi-leptonic heavy flavor decay by PHENIX [175] and STAR [HE]. Fig. lOl 
illustrates the first PHENIX results from central Au+Au at y/s = 200 GeV, 
Raa falling to about 0.3 at 5 GeV/c. Heavy flavor attenuation thus resembles 
that of light quarks, as is also attested by the predictions of in-medium parton 
transport models [179^1180] included in Fig. 15.31 which cast the medium opacity 
into an effective parton transport coefficient q which is seen here to approach a 
value 14 GeV'^/c at high prp, again corresponding to a highly opaque medium. 
We shall describe this approach in more detail below. 



A most fundamental cross-check of the in-medium attenuation picture of color 
charged partons consists in measuring Raa for primordial, "direct" photons. 
Fig. 15.41 shows first PHENIX results [181J for direct photon Raa vs. pt in 
central Au+Au at ^/s = 200 GeV. Photons obey pQCD number of collisions 
scaling, Raa ~ 1! Also included in Fig. l5.4l are the attenuation ratios for neutral 
pions (already shown in Fig. 15. 2p . and for r/ mesons that follow the pattern of 
extreme suppression. In essence, the PHENIX results in Figs. 15.31 and 15.41 wrap 
up all one needs to know for a theoretical analysis of fireball medium opacity, 
for various fiavors, and indicate transparency for photons. 
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Figure 5.3: The nuclear modification factor Raa vs. pt for electrons from semi- 
leptonic decays of heavy flavor (mostly D) mesons in central Au+Au 
collisions at ^/s = 200 GeV |175j ; with calculations of in medium 
energy loss using different attenuation models [1791 I180J . 



We turn, briefly, to Rcp as opposed to Raa analysis, in order to ascertain sim- 
ilar resulting conclusions at RHIC energy. Fig. 15.5! illustrates an Rcp analysis 
of IT, p and charged hadron high px production in Au+Au at ^/s = 200 GeV 
by STAR [182J . the central to peripheral yield ratio referring to a 5%, and a 
60-80% cut of minimum bias data. We showed in equation 15.61 that the Rcp 
measure also refers to a picture of pQCD number of binary collision scaling, 
inappropriate at low pT- Thus ignoring the features of Fig. [531 at px < 3GeV/c 
we conclude that the high pT data again suggest a suppression by about 0.3, 
common to pions and protons, thus approaching the ratio, of about 0.2, ob- 
served in Figs. 15. Il and 15.2] which employ the "ideal" in- vacuum p+p ^ hadron 
+ X reference. 



At top SPS energy, ^/s = 17.3GeV, the experimentally covered px range is fairly 
limited |183j . pr < ^GeV/c. Fig. l5.6l shows NA49 results, Rcp for p and charged 
pions. Contrary to former expectations that such data would be overwhelmed 
by Croonin-enhancement |184j the same systematic behaviour as at RHIC is 
observed, qualitatively: Rcp (baryon)> Rcp (meson) at px > 3GeV/c. Note 
that, again, the data do not approach unity at px ^ because of the employed 
binary scaling, and that the strong rise of the proton signal at px < 2 GeV jc 
is largely the result of strong radial flow in central Pb+Pb collisions, an ef- 
fect much less prominent in pion px spectra. The high px suppression is much 
weaker than at RHIC but it is strong enough that the expected Croonin en- 
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Figure 5.4: Modification factor Raa vs. pt for direct photons in central Au+Au 
at 200 GeV (squares) . Also shown are Raa for tt" (triangles) and rj 
(circles) |181| fitted by the attenuation model |176t I180| . 
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Figure 5.5: Rcp, the ratio of scaled vr, p and charged hadron yields vs. pT in 
central (5%) and peripheral (60-80%) Au+Au collisions at ^/s = 
200 GeV [152] . 



hancement of high pT mesons is not observed. These SPS data, as well as 
first results obtained at the intermediate RHIC energy of y/s = 62.4 GeV |185j 
are reproduced by an attenuation model based on the primordial gluon density 
dNg/dy that scales as the charged particle midrapidity density dNch/dy |176j . 
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and was also employed in Figs. 15.31 and 15.41 
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Figure 5.6: Rcp results from SPS Pb+Pb collisions at y/s = 17.3GeV, for pions 
and protons |T85j . with attenuation model fits jl76j. 



5.2 Energy loss in a QCD medium 

The attenuation model that we have hinted at consists of a gluon radiative en- 
ergy loss theory of the primordially produced leading, high px parton as it tra- 
verses a medium of color charges, by means of emission of gluon bremsstrahlung. 
We expect that the resulting partonic specific energy loss, per unit pathlength 
(i.e. its dE/dx) should reflect characteristic properties of the traversed medium, 
most prominently the spatial density of color charges [176j but also the average 
momentum transfer, per in-medium collision of the considered parton or, more 
general, per unit pathlength at constant density. Most importantly, an aspect 
of non-abelian QCD leads to a characteristic difference from the corresponding 
QED situation: the radiated gluon is itself color-charged, and its emission prob- 
ability is influenced, again, by its subsequent interaction in the medium |186] 
which, in turn, is proportional to medium color charge density and traversed 
pathlength L. Thus, the traversed path-length L in-medium occurs, both in 
the probability to emit a bremsstrahlung gluon, and in its subsequent rescat- 
tering trajectory, also of length L, until the gluon finally decoheres. Quantum 
mechanical coherence thus leads to the conclusion that non-abelian dE/dx is 
not proportional to pathlength L (as in QED) but to L^ [187j . 

This phenomenon occurs at intermediate values of the radiated gluon energy, w, 
in between the limits known as the Bethe-Heitler, and the factorization regimes 
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LOBH ~ A g|. « o; < Ljfact ~ L'^ Qt/^ < E (5.7) 

where A is the in- medium mean free path, q^ the (average) parton transverse 
momentum square, created per coUision, and E the total cm energy of the 
travehng charge. In the BDMPSZ model [Ml HHZl [Mj the properties of the 
medium are encoded in the transport coefficient, defined as the average induced 
transverse momentum squared per unit mean free path, 

g=(gr>/A. (5.8) 

The scale of the radiated gluon energy distribution uj dN/duj is set by the 
characteristic gluon energy [1861 I187J 

LU, = ^qL\ (5.9) 

To see more explicitly how the various properties of the color charged medium 
enter in q we rewrite it as 

Or dqr ^ = pcr (q^) = A"^ (q^) (5.10) 

where p is the color charge density of scattering centers, a the effective binary 
cross section for interaction of the considered leading parton at scale q'^ (which 
may depend on quark fiavor), and (^q^) as above. Obviously, both a and (g|.) 
refer to detailed, intrinsic properties of the QCD medium, globally described by 
density p. The leading parton cross section with in-medium color charges should 
depend on the implied resolution scale, Q^ = (gfi), and can thus be obtained 
from perturbative QCD |186lll87U188j only if Q^ > Q^at' the saturation scale 
that we discussed in section [2j Likewise, (g|n^ itself refiects a medium property, 
the effective range of the color force, which is different in confined and decon- 
fined media. Hadron size limits the force range in the former case, such that 
q is minimal in ground state hadronic matter also, of course, due to the small 
energy density p = po = 0.15 GeV/fm^ [189] . This was, in fact, confirmed by a 
RHIC run with deuteron-gold collisions, in which mid-rapidity hadrons traverse 
sections of cold Au nucleus spectator matter. Fig. 15.71 shows results obtained 
for Riij^ dependence on px, for vr*^ from PHENIX [190j . and for charged hadrons 
from STAR [191]. For comparison, both panels also include the corresponding 
Raa data for central Au-|-Au collisions (all at ^/s = 200 GeV/c), exhibiting 
the typical, drastic high pT quenching of inclusive hadron production, clearly 
absent in d+Au collisions. 

We have shown a first application of the BDMPSZ model, to RHIC inclusive D 
meson production [179j, in Fig. 15.31 Before engaging in further model applica- 
tion we note, first, that equations 15. 7115710] above refer to the idealized conditions 
of an infinitely extended medium of uniform composition. In reality, the fire- 
ball medium expands, at time scale concurrent with the proper time incurred 
in the leading partons propagation over distance L, such that all ingredients in 
q, exhibited in equation 15.101 vary with expansion time [192j . However, before 
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Figure 5.7: Raa vs. pt for d+Au collisions at ^/s = 200 GeV, compared to 
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for vr^ (left) and charged hadrons (right). 



turning to adaption to reality of the infinite matter model, we wish to expose 
its prediction for the final connection of the specific average partonic energy 
loss {AE) and in-medium path length L traversed: 



(A^) 



CO— — duj oc as t-R q L . 
duj 



(5.11) 



This relation |193j represents the eikonal approximation limit of extended medium 
and high leading parton initial energy E > uJc (equation 15. 9p . The average en- 
ergy loss is thus proportional to the appropriate strong coupling constant a^, 
to the Casimir factor corresponding to the leading parton (with value 4/3 for 
quarks, and 3 for gluons), as well as to q and L^. 

In order to verify the non-abelian nature of radiative parton energy loss in a 
partonic QCD medium it would, of course, be most convincing if a direct, ex- 
plicit L^ dependence of (AE) could be demonstrated. Such a demonstration 
is lacking thus far, chiefiy because of the obvious difficulty of simultaneously 
knowing the primordial parton energy E, the transport coefficient q and - in 
finite nuclear collision geometry - its variation along the actual traversed path 
L as the surrounding medium expands with propagation time. Moreover, the 
partonic medium induced energy loss AE of the primordial parton is not di- 
rectly observable. Even if we assume that high pT partons evolve into the 
observed hadrons only after leaving the fireball medium |176j . their ensuing 
"fragmentation" to hadrons (which is known from p + p jet physics) results 
in several hadrons usually comprising a "leading" hadron which transports a 
major fraction (z) = (^E /E^'j of the fragmenting parton energy E^, which, 
in turn, equals E^ (primordial) - AE, with AE sampled from a probability 
distribution with mean {AE) according to equation 15. Ill The observed leading 
hadron energy or transverse momentum is thus subject to sampling, both, z 
from the fragmentation function, and AE from in-medium energy loss. Finally, 
inclusive high pT leading hadron observation in A+A collisions involves an aver- 
age over all potential initial parton production points, within the primordially 
produced density profile. A specific distribution of in medium path lengths 
/(L) arises, for each such production point, which, moreover, depends on a 
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model of space-time fireball expansion. The final inclusive yield thus requires a 
further, weighted volume average over f{L) per production point. Thus, typical 
of an inclusive mode of observation, the "ideal" relationship of equation 15. IH 
between radiative in-medium energy loss AE and traversed path length L gets 
shrouded by double averages, independently occuring at either side of the equa- 
tion [Ml [1791 [EHl [IMl [Ml [196]. A detailed L^ law verification can not be 
expected from inclusive central collision data alone (see next section). 

However, the unmistakeably clear signal of a strong, in- medium high pT par- 
ton quenching effect, gathered at RHIC by Raa measurement for a multitude 
of hadronic species (Figs. [EH \EM 12391 EH \SM EZj), in Au+Au collisions at 
y/s = 200 GeV, has resulted in first estimates of the transport coefficient q, the 
medium - specific quantity entering equation 15. Ill in addition to the geometry - 
specific path length L. In fact, the transport coefficient can, to some extent, be 
analyzed independently, owing to the fact that q (x g from equation I5.1U[ The 
density g falls down rapidly during expansion, but it is initially rather well con- 
strained by the conditions of one-dimensional Bjorken expansion that we have 
described in chapters [2] and [H The major contribution to partonic AE arises 
in the early expansion phase (via a high q), in close analogy to the formation of 
the elliptic flow signal. These two signals are, thus, closely correlated: the pri- 
mordial hydrodynamic expansion phase of bulk matter evolution sets the stage 
for the attenuation, during this stage of QCD matter, of primordially produced 
"tracer" partons, traversing the bulk matter medium as test particles. 

The bias in partonic AE to the primordial expansion period is borne out in 
an expression [1931 I195J which replaces the q coefficient, appropriate to an 
infinitely extended static medium considered in the original BDMPSZ model, 
by an effective, expansion time averaged 

Qeff = ^j dtit-to)q{t) (5.12) 

to be employed in the realistic case of an expanding fireball source. Due to the 
rapid fall-off of g, in q = ga (g|i) from equation 15.101 the integral depends, far 
more strongly, on g(t ~ to) than on total path length L. Furthermore, inserting 
(Jeff into the BDMPSZ formula [1931 I195J for the transverse momentum down- 
ward shift, occuring in leading parton or hadron pT spectra (of power law form 
p:^", see Fig. [2l9]) 

Apr ~ -assjT^qL'^PTJv, (5.13) 

we see that the first order proportionality to L^ is removed. The downward pT 
shift is thus, primarily, a consequence of q^ff which, in turn, is biased to reflect 
the "ideal" transport coefficient q at early evolution time. Within this termi- 
nology,the pT shift (see equation 15.130 determines the experimentally measured 
ratio RAAij>T) which quantifies the effective transport coefficient q^ff for the 
Pt domain considered. It can be related, as a cross check, to the initial gluon 
rapidity density if the collision region expands according to Bjorken scaling 
[1871 [T971: 

,^„.|fl-^. (5.14) 



5.2 Energy loss in a QCD medium 



A typical result of application of the model described above jl95j is shown in 
Fig. 15. 8[ Analogous to Fig. 15.11 Raa for neutral pions and charged hadrons 
is averaged over the range 4.5 < pt ^ 10 GeV/c, and shown as a function 
of centrality (assessed by Npart) in minimum bias Au+Au collision at ^/s = 
200 GeV UMllMllIM]. A path-averaged q^ff of UGeV^/fm is inferred from 
the fit, in close agreement to the value found in |195j . 

A more recent study |199j of the PHENIX Raa data for vr*^ in central Au+Au 
collisions (Fig. 15. 2p is shown in Fig. 15.91 The analysis is carried out in the 
framework of the WHDG model [200j . which replaces the (effective) transport 
coefficient q (employed in the BDMPSZ model [186^ 11871 I188J . and turned 
into the data analysis formalism of |193t I195J ) by the primordial gluon mid- 
rapidity density dN^ /dy, as the fundamental parameter, via equation 15.141 
The initial gluon density, in turn, being related to the charged hadron mid- 
rapidity density [61] . Fig. 15.91 shows that, within the still preliminary statistics 
at pt > 10 GeV/c, the "conservative" estimate of q^ = 0.3 and dN^ /dy = 1000 
does not appear to be the most appropriate choice, the data rather requiring 
1000 < dN^ /dy < 2000. Overall, Fig. 15.91 demonstrates a certain insensitivity 
of the data, to the basic parameters of theoretical high pj- quenching models 
of inclusive hadron production at RHIC energy, that we have already inferred 
from Fig. 15. 3^ concerning choice of q. 

Radiative in-medium energy loss of primordially produced partons, traversing 
the evolving bulk medium as "tracers", must be extremely strong. This is in- 
dicated by the inclusive attenuation ratios Raa, which fall down to about 0.2 
in central collisions thus almost reaching the absolute lower limit of about 0.15 
that arises from the unavoidable fraction of un-attenuated primordial surface 
"corona" nucleon-nucleon interaction products |201j . 
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Figure 5.8: The effective transport coefficient q = 14 GeV"^ / fm in the parton 
quenching model (PQM) of |195j determined from the centrality 
dependence of Raa for vr'' and charged hadrons, averaged over 4.5 
<PT< 10 GeV/c, in Au+Au at ^ = 200 GeV. 
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Figure 5.9: Application of the WDHG transport model |200j based on equa- 
tion EHH to PHENIX Raa data for vr'^ [175], indicating primordial 



1000 < ^ < 2000. 



Not expected from early estimates of medium opacity based on the picture of 
a weakly coupled perturbative QCD medium [172^ 1202] . the interior sections 
of the central Au+Au interaction volume must be almost "black" toward high 
Pt parton transport [1941 I195J at \/i = 200 GeV, also including charm quark 
propagation (Fig. 15.3]) . The remaining signal should thus stem, primarily, from 
the dilute surface sections, and from the finite fraction of partons traversing 
the interior with small, or zero radiative energy loss, as a consequence of the 
finite width of the AE probability distribution }193tll94] . Seen in this light, the 
smooth decrease of Raa with centrality (Fig. 15. 8p should reflect the combined 
effects, of a decreasing surface to volume ratio, an increasing effective q (due to 
interior density increase) that confronts the increasing average geometrical po- 
tential path length (L) (essentially enhancing its effect), and a thus diminishing 
fraction of primordial high pT partons experiencing a small AE. 

Not surprisingly, the ideal non abelian QCD relationship of AE proportional 
to in-medium high px parton path length L^ can, thus, not be established from 
inclusive high px quenching data alone. We shall show in the next section that 
di-jet primordial production can offer a mechanism of higher resolution. The in- 
clusive Raa attenuation data, obtained at RHIC, are seen to establish an unex- 
pected, high opacity of the primordial interaction volume, extending to high px 
parton propagation. The required, high transport coefficient q = lAGeV'^ / fm 
from Fig. 15.81 confirms and extends the picture derived from elliptic flow data 
|167) : at top RHIC energy the plasma is non-perturbatively, strongly coupled, 
a new challenge to lattice QCD [61j. The QGP may be largely composed of 
colored, string- like partonic aggregates [203j . 
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5.3 Di-jet production and attenuation in A+A collisions 

In order to analyze leading parton attenuation in a more constrained situation 
[2U4j , one investigates parton tracer attenuation under further geometrical con- 
straints concerning the in-medium path length L, by means of di-jet analysis, 
and/or by varying the primordial parton density that enters q via equation I5.1UI 
in studies at different y/s while maintaining the observational geometrical con- 
straints. 

We shall concentrate here on di-jet attenuation data obtained in Au+Au col- 
lisions at top RHIC energy, ^/s = 200 GeV . At this relatively modest energy 
the initial pQCD production cross section of leading partons (as described in 
equation IS.ip reaches up to pT = 25 GeV/c. The ensuing DGLAP shower 
multiplication initiates " parton fragmentation" to hadrons [83t [5^ I121| , each 
carrying a momentum fraction zt = PT/PT(primord. parton). The created en- 
semble of hadrons h belonging to the observed hadronic jet can be wrapped up 
by the total fragmentation function 

/dz 
— Ci {Z, \/s)Dpart-.h{z, \/s) (5.15) 

•I 

which summarizes the contributions arising from the different shower partons i. 
Here, Ci are the weight coefficients of the particular process, and -Dpart(j)-*/i are 
the individual fragmentation functions (FFs) for turning parton i into hadron 
h. Similar to the parton distribution functions (PDFs) in equation 15.11 derived 
from deep inelastic electron-parton scattering (DIS) and shown in Fig. 12.13} the 
FFs are semi-empirical non-perturbative QCD functions that have an intuitive 
probabilistic interpretation. They quantify the probability that the primordial 
parton produced at short distance 1/Q fragments into i shower partons, forming 
a jet that includes the hadron h [ 20511206] . 

At Fermilab energy, y/s = 1.8 TeV , the jet spectrum reaches up to Et ~ 
400 GeV , and a typical 100 GeV jet comprises about 10 hadrons which can be 
identified above background by jet-cone reconstruction algorithms |205| . This 
allows for a complete determination of the corresponding fragmentation func- 
tion, and for a rather accurate reconstruction of the pT and Et of the primordial 
parton that initiated the jet. Similar conditions will prevail in jet spectroscopy 
of Pb-|-Pb collisions at LHC energy, ^/s = 5.5 TeV . 

However, at RHIC energy a typical jet at 15 < Et < 25 GeV features a frag- 
mentation function comprised of a few hadrons with Et in the 2-15 GeV range. 
Considering the high background, arising in the lower fraction of this energy 
domain from concurrent, unrelated high pT hadron production processes, a com- 
plete jet-cone analysis can not succeed. The RHIC experiments thus confront 
back-to-back di-jet production with an analysis of the azimuthal correlation be- 
tween high pt hadrons. Defining the observational geometry by selecting a high 
Pt "trigger" hadron observed at azimuthal emission angle ^trig-, the associated 
production of high pT hadrons is inspected as a function of Ac^ = ipass — Vtrig- 
If the trigger has caught a leading jet hadron one expects the hadrons of the 
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balancing back-to-back jet to occur at the side opposite to the trigger, Aip ~ vr. 
The trigger condition thus imposes the definition of a " near-side" and an " away 
side" azimuthal domain. Furthermore, the relatively narrow rapidity accep- 
tance of the STAR and PHENIX experiments (centered at y = 0) selects di-jets 
with axis perpendicular to the beam direction. 

Originating from a uniform distribution of primordial back-to-back di-parton 
production vertices, throughout the primordial reaction volume, the trigger 
selected di-jet partons thus experience an (anti-) correlated average path length 
(L) to arrive at the surface while experiencing medium-specific attenuation, 
with {Ltrig) ~ 2R — (Laway), R being the transverse medium radius. No such 
geometric constraint exists in the study of inclusive high px hadron production. 
We thus expect information different from the inclusive Raa{pt) signal. The 
geometrical selectivity can be even further constrained by fixing the direction of 
the impact parameter (i.e. the reaction plane) in semi-central collisions (recall 
chapter H]), and observing the di-jet correlation signal in dependence of the 
di-jet axis orientation relative to the reaction plane. 

The very first di-hadron correlation measurements confirmed the existence of 
strong in-medium attenuation. Fig. 15.101 shows the azimuthal yield distribu- 
tions, per trigger hadron, as observed by STAR at y/s = 200 GeV pO?] . The 
left panel shows the distribution of hadrons with pt ^ '^ GeV/c relative to a 
trigger hadron with p^*^ > 4 GeV/c. Data for p-|-p, d-|-Au and central Au-|-Au 
are illustrated. At the "near side" (the trigger position is $ = 0) all three 
reactions exhibit a similar narrow distribution of hadrons associated with the 
trigger hadron, typical of a jet cone fragmentation mechanism. Note that the 
associated near-side central Au+Au signal thus exhibits no signs of an atten- 
uation softened fragmentation function, indicating that the trigger imposed 
high pt hadron should predominantly stem from primordial jet production ver- 
tex points located near to the surface of the reaction volume, in azimuthal 
trigger direction. Thus, conversely, the balancing opposite jet has to develop 
while traversing almost the entire transverse diameter of the interaction volume. 
I.e. (Loppos) ~ 2 i? thus emphasizing the expectation that di-jet spectroscopy 
should allow for stricter constraints on path length L in comparison to single 
high pt hadron Raa analysis. In fact, no trigger related away side signal of 
Pt > 2 GeV/c hadrons is observed in Fig. 15.101 for central Au-|-Au collisions, 
whereas p+p and central d+Au collisions exhibit a clear away-side di-jet signal. 

We conclude that the trigger bias, selecting a single near side hadron of pT > 
4 GeV jc in central Au-|-Au collisions, responds to a primordial di-jet of about 
10 GeV per back-to-back parton. After traversal of in medium average path 
length L ^> 2R the fragmentation function of the opposite side parton contains 
on average no hadron at px > 2 GeV/c, indicating that it should have lost a 
fraction (AEt) > 5 GeV. The medium is thus highly opaque, but the total 
disappearance of the opposite side signal can only provide for a lower limit es- 
timate of {AEt), within the trigger conditions employed here. We shall show 
below that the situation changes with more recent RHIC data |208j that extend 
the trigger hadron pj- range toward 20 GeV/c. 
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Figure 5.10: Di-hadron correlation from back-to-back di-jet production in 
Au-|-Au collisions at y/s = 200 GeV. The trigger particle is at 
azimuth <1> = 0, with a pT > ^GeV threshold. The away side peak 
at A$ = vr is observed (left panel) in p+p, d+A but absent in 
central Au+Au. Right panel shows the correlation in Au+Au for 
different orientations of the trigger direction relative to the reaction 
plane [W]. 



However, the right hand panel of Fig. 15.101 shows that an improved constraint 
on partonic in-medium path length can already be obtained by studying the 
di-jet back-to-back production geometry in correlation with the orientation of 
the reaction plane that arises from non-zero impact parameter in semi-central 
A+A collisions. We have seen in Fig. 14.11 that such collisions exhibit an ellip- 
tical primordial transverse density profile, with minor axis along the impact 
vector 6, defining the reaction plane. Di-jets with axis "in-plane" thus traverse 
a shorter in-medium path length as orthogonal " out-of-plane" jets, the differ- 
ence in average path length being quantified by the spatial excentricity e(6), 
equation 14. 1[ Fig. 15.101 (right) shows the di-hadron correlation results in semi- 
peripheral Au+Au collisions, as compared to the in-vacuum p+p reference. At 
the 20-60% centrality window employed here, out of plane jet emission occurs 
along the major axis, the reaction volume diameter still of magnitude 2R, as in 
central collisions. The trigger condition thus again selects opposite side path 
lengths L — > 2 i?, but the energy density should be lower than in central colli- 
sions. Even so, the average opacity along the away side parton path appears to 
be high enough to wipe out the correlation signal. In-plane geometry, however, 
shows a partially attenuated signal (as compared to the global p+p reference) 
at the opposite side, corresponding to path lengths L ^ R. These data thus 
provide for first information concerning the relation of {^E) and average tra- 
versed path length [176| . 



Obviously, one wants to extend the above study to higher di-jet energies, i.e. 
to measurement of di-hadron correlations at hadron trigger pT — > 20 GeV/c, 
conditions sampling the very low primordial cross section, at ^/s = 200GeV, of 
jet production at primordial Et -^ 30 GeV [208. .209]. Fig. 15.111 shows the cor- 
responding jet correlations selected with high pT trigger, 8 < p^^^ < 15 GeV/c, 
and high pT associated hadrons, pT > Q GeV/c, in minimum bias d+Au, semi- 
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Figure 5.11: Di-hadron correlation at high px in central Au+Au collisions at 
^/s = 200 GeV, compared to d+Au and peripheral Au+Au; for 
8 < p^*^ < 15 GeV, and pf?"°^ > 6 GeV. From [208]. 



central Au+Au and central Au+Au at ^/s = 200 GeV. Very clear and almost 
background-free back-to-back jets are observed in all three cases, in sharp con- 
trast with the away side jet total disappearance in central Au+Au at lower jet 
energy, Fig. 15.101 (left panel). The near side trigger associated hadron yield 
decreases only slightly from d+A to central Au+Au, while the away side yield 
drops down by an attenuation factor of about 0.2, the signal thus not being 
completely extinguished. We infer Raa{L -^ 2R) ^ 0.2. 

In order to show how such data can be evaluated in a picture of in-medium 
leading parton attenuation (as finally reflected in leading hadron production 
observed in the above di-hadron correlation data) we briefly consult the pQCD 
factorization |210j prediction for the inclusive production of a high pT hadron 
at central rapidity, in the nuclear collision A + B ^ h + x [ 196j , 
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(5.16) 



where the parton (a, b) distribution functions F in nucleus A, B and the el- 
ementary pQCD cross section for a + b ^> c + x have been already implied 
in equation 15.11 Their spatial integral gets convoluted with the fragmentation 
function D that describes the conversion of the leading parton c to a hadron 
carrying a fraction < Zc < 1 of its transverse momentum, i^ is a factor in- 
troduced as a phenomenological correction for "next to leading order" (NLO) 
QCD effects. Within the (further) approximation that the leading parton c suf- 
fers medium induced gluon bremsstrahlung energy loss but hadronizes outside 
the interaction volume (in vacuum), the in- medium quenching leads, merely, to 
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a re-scaling of the fragmentation function, 



Timed 



jdeP{e)-^^Dl% (y^,Q'), (5.17) 



where the primary parton is imphed to lose an energy fraction e = /\E/Ec with 
probability P{e) [196J. Therefore the leading hadron is a fragment of a parton 
with reduced energy (1 — e)£'c, and accordingly must carry a larger fraction of 
the parton energy, Zc/{1 — e)- If no final state quenching is considered, P{e) 
reduces to 5{e). The entire effect of medium attenuation on the leading parton 
is thus contained in the shift of the fragmentation function. 

An application of this formalism ^196j is shown in Fig. 15.121 The in-medium 
modification of the hadron-triggered fragmentation function (see equation l5.17p 
is evaluated for central Au-|-Au collisions at ^/s = 200Gey. In adaptation to the 
modalities of RHIC di-hadron correlation data, the opposite side fragmentation 
function (for observation of trigger-related hadrons with px > 2 GeV/c) is stud- 
ied in dependence of the trigger selected pT window. Its attenuation is quanti- 
fied by the ratio D (with quenching) to D (without quenching), as a function 
of the fraction zt, of opposite side hadron px to trigger hadron px- Referring 
to the observational conditions implied in Fig. I5.1UI (left) p^^^ « 4 — 6 GeV/c 
and opposite side pT > 2 GeV, the predicted suppression amounts to a fac- 
tor of about 0.35, whereas the data imply a factor smaller than 0.2. Likewise, 
comparing to the data of Fig. 15.111 the (harder) trigger conditions should lead 
to a suppression of about 0.45 from Fig. 15.121 but the observed value is close 
to 0.2. The predictions appear to fall short of the actually observed suppres- 
sion. This calculation employs an ansatz for the transport coefficient q similar 
to equation I5.14[ recurring to the primordial gluon density at tq < 0.6 fm/c 
which, in turn, is estimated by the charged hadron mid-rapidity density [197]. 
However, this pQCD based argument can also not reproduce the magnitude of 
the effective transport coefficient (equation I5.12p . q^ff ss 10 — 15 GeV'^ /fm, 
shown in refs. |194l \19b\ I199| to be required by the large observed suppression 
(see Fig. [OD. 

It has been argued [1991 I211J that these models need refinement by introducing 
more realistic dynamics, and/or by completely abandoning the pQCD ansatz for 
hard parton in-medium transport J2121I213"] . We shall return to these novel sug- 
gestions, of how to treat dynamical, non equilibrium quantities of the "strongly 
coupled" parton plasma of non perturbative QCD (toward which equilibrium 
lattice theory can only give hints) , in our final conclusion (chapter [8]) . In the 
meanwhile, we note that the expected non abelian behaviour, AE oc L"^, could 
not be verified quantitatively, as of yet [STlj . because of the unexpectedly high 
opacity of the fireball interior sections at ^/s = 200 GeV, in combination with 
the limited jet energy range that is available at RHIC energy. It appears possi- 
ble, however, to extend the analysis of the "back side jet re-appearance" data 
[2081 1209] (Fig. [5TT]1 toward this goal. The situation should improve at LHC 
energy, y/s = 5.5 TeV, where primordial 100-200 GeV jets are abundant, such 
that the opposite side jet can be reconstructed with explicit use of the complete 
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Figure 5.12: The ratio of the hadron-triggered fragmentation function, equa- 
tion 15.171 in central Au+Au and in p+p collisions, for different 
values of p^ *^ [196] . 



Fermilab jet cone recognition algorithms |2U5j even if their in-medium Et loss 
ranges up to 50 GeV. 

A further prediction of the model employed in Fig. 15.121 has been confirmed by 
the RHIC experiments. The high medium opacity at top RHIC energy leads 
to an intriguing emission pattern of low pT opposite side hadrons. Clearly, the 
trigger-selected Et flux, of up to 20 GeV, toward the away-side, can not re- 
main unnoticeable. Inspection of the attenuated fragmentation functions |196j 
in Fig. 15.121 reveals an enhanced emission of bremsstrahlung gluon hadroniza- 
tion products at zy < 0.1 This fraction of in-medium jet-degradation products 
has in fact been observed, as is shown in Fig. 15.131 The left panel shows STAR 
results [214J for the di-hadron correlation in central Au-|-Au at ^/s = 200 GeV, 
with near-side hadron trigger 4 < p^ < 6 GeV/c, and opposite side observation 
extended to soft hadrons, 0.15 < pr < 4 GeV/c. A prominent double-peak 
structure is indicated, symmetric about A$ = vr. The right panel shows high 
resolution PHENIX results [215] for Au-|-Au at three centralities, from periph- 
eral to central collisions. For the former, the typical p-|-p-like away side peak 
(c.f. Fig. 15. 7|) is recovered, while a double peak appears in the latter case, 
shifted from A$ = vr by ib(5<l> ~ 70*^. A hypothetical mechanism comes to mind 
[213] ■ of sideward matter acceleration in a "Mach-cone" directed mechanism 
of compressional shock waves initiated by the in-medium energy loss of the 
opposite side leading jet parton, which traverses the medium at "super-sonic" 
velocity, i.e. at f > u^, the appropriate speed of sound in a parton plasma. 

If confirmed by pending studies of the away-side multi-hadron correlation, that 
might ascertain the implied conical shape of the soft hadron emission pattern 
(about the direction A<I> ~ vr of the leading parton), this mechanism might lead 
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Figure 5.13: Di-hadron correlation: 



to the determination of the sound (or shock wave) velocity of a strongly coupled 
parton plasma: a third characteristic QGP matter property, in addition to vis- 
cosity T] (from elliptic flow) and q (from high pT parton attenuation). We note 
that the implied concept, of measuring the shock wave transport velocity of a 
strongly interacting medium, dates back to the 1959 idea of Glassgold, Heck- 
rotte and Watson [216], to study shock wave emission in central p+A collisions 
at AGS and PS proton energy, of about 30 GeV . In hindsight we can under- 
stand the inconclusiveness of such searches: the "low energy" incident proton 
does not preserve its identity, breaking up into relatively soft partons in the 
course of traversing the target nucleus. Their coupling to the surrounding cold 
baryonic target matter is weak, and possible collective transverse compressional 
waves are swiftly dissipated by the cold, highly viscous [61] hadronic medium. 
The shock compression and Mach-cone emission mechanism was then revived 
by Greiner and collaborators [217j for central "^^Ne + U collisions at Bevalac 
energy |218| . The recent RHIC observations in Fig. 15. Rlj have been called "con- 
ical flow" [219j : they demonstrate, again, the strong coupling to the medium, 
of even a 20 GeV jc leading parton |220j . 
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6 Vector meson and direct photon 
production: penetrating probes 



This chapter is devoted to three observables that have been profoundly studied 
already at the SPS, as tracers of the created fireball medium: 

1. J/^ "charmonium" production, as suppressed in a high T QGP medium, 

2. "direct" photons as black body T sensors of the QGP state; 

3. p meson in medium production, studied by di-lepton decay. 

These three observables represent an internally connected set of ideas of high 
density QCD matter diagnostics: all serve as medium tracers but in a comple- 
mentary way. Matsui and Satz [31] realized that at the modest top SPS energy, 
17.3 < ^/s < 20 GeV, the production rate of cc pairs (that would in part de- 
velop toward charmonium production, J/^,^ , etc.) was so low that only the 
primordial, first generation nucleon-nucleon collisions in an A+A event would 
create a measurable yield (nothing after thermalization) . I.e. the primordial 
cc yield would be well-estimated by A'^'^ app{cc). Initially produced in a su- 
perposition of states including color triplets and octets |221) the emerging cc 
pairs thus co-travel with the developing high energy density fireball medium, 
as tracers, on their way toward J/^ or ^ and D, D formation. Attenuation, 
characteristic of medium properties, will break up leading cc pairs resulting 
in a suppression of the eventually observed J/^ yield: another "quenching" 
observable. 

As the suppression of J/ip is related to medium temperature and density, the 
extent of charmonium quenching should also be related to the rate of black 
body thermal fireball radiation, by photon emission via elementary qq — > 57 
and qg -^ q^ processes in the plasma |2221 1223] . Photons leave the interaction 
volume essentially un-rescattered, and their radiative intensity, proportional to 
T^, makes them an ideal probe of the initial fireball temperature. This thus 
could be an ideal diagnostics of the early deconfined matter, but it is difficult 
to disentangle from a multitude of concurrent low pT photon sources |224j . 
most prominently vr'^ -^ 27 decay, with cross sections higher by several orders 
of magnitude. The thermal photon signal thus becomes more promising the 
higher the initial temperature Tj which might reach up to 500-600 MeV at 
LHC energy. 

Similar to photons, in medium created lepton pairs |225j escape essentially 
un-attenuated as was shown for the Drell-Yan process, qq -^ LL by CERN 
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experiment NA38/NA50 |226j . Thermal di-lepton production in the mass re- 
gion < 1 GeV is largely mediated by light vector mesons. Among these, the p 
meson is of particular interest due to its short lifetime (1.3 fm/c), making it an 
ideal tracer (via its in-medium di-lepton decay) for the modification of hadrons 
composed of light quarks, in the vicinity of T = Tc. This modification signal is 
thus complementary to the heavy quark charmonium J/^ suppression (break 
up) effect that sets in at T > 1.5 — 2Tc (see below). Moreover, in addition 
to the deconfinement breakup mechanism by QCD plasma Debye screening of 
the color force potential acting on the cc pair [fflj, the QCD chiral symmetry 
restoration mechanism [227] can be studied via in-medium modification of the p 
spectral function as T —> Tc. Note that the in vacuum p mass and width prop- 
erties are owed to non-perturbative QCD condensate structures [U I228J which 
spontaneously break the chiral symmetry of the QCD Lagrangian, at T — > 0. 
These properties should change, in the vicinity of Tc, and be reflected in mod- 
ifications of the di-electron or di-myon decay spectra - unlike the suppression 
effect on the J/^ which simply dissolves primordial cc pairs before they can 
hadronize, a yes-no-effect whose onset with y/s or centrality serves as a plasma 
thermometer, by observing Raa or Rcp < 1- 

6.1 Charmonium Suppression 

Due to the high charm and bottom quark masses, the "quarkonium" states of 
cc and bb can be described in non-relativistic potential theory [229^1230] . using 

V{r) = ar-- (6.1) 

r 

as the confining potential [231j . with string tension a = 0.2 GeV^ and gauge 
coupling a = tt/12. We are interested in the states J/^ (3.097), Xc (3.53) and 
^ (3.685) which are the IS, IP and 2S levels. The decay of the latter two feeds 
into the J/^, accounting for about 40% of its yield. The radii derived from 
equation 16. II are 0.25 fm, 0.36 fm and 0.45 /m, respectively, well below hadron 
size at least for the J/^ and Xc states. 

With increasing temperature, (7{T) decreases, and at deconfinement cr(Tc) = 0. 
For T > Tc we thus expect 

o 

V(r) = exp l-r/rD (T)l (6.2) 

r 

where rj:,{T) is the QCD Debye screening radius. It was initially estimated 
from a SU (2) gauge theory of thermal gluons [H], to amount to about 0.2-0.3 
/to at T/Tc = 1.5. In this picture, the screened potential (equation 16. 2|) can 
still give rise to bound cc states provided their radius is smaller than tq. The 
pioneering study of Matsui and Satz [H] concluded that screening in a QGP 
system would dissolve the J/^, or its cc precursor, at T > 1.3 Tc whereas the 
Xc and ^ states would be suppressed already directly above Tc. 

The corresponding energy density for J/^ suppression, employing 

e/ec ^ {T/Tcf « 2.9 (6.3) 
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(obtained with ec ~ 1 GeV/fm^ from lattice QCD), would thus amount to 
about 2.9 GeV/fm^. This motivated an extensive experimental effort at the 
CERN SPS Pb beam at ^/s = 17.3 GeV. We have seen in chapter [2] that the 
Bjorken estimate |45] of average transverse energy density reached in central 
Pb+Pb colhsions [SI E] amounts to e = (3.0 ± 0.6) GeV/ fm?, with higher e 
to be expected in the interior fireball sections: encouraging conditions. 

However, the above approach assumes the validity of a two-body potential 
treatment at finite T, near a conjectured critical point of QCD. More recently 
the quarkonium spectrum was calculated directly in finite temperature lattice 
QCD |232j . with the striking result that the J/^ dissociation temperature in 
a realistic non-perturbative treatment of the QCD plasma state moves up to 
about T = 2 Tc, whereas Xc and ^ dissociation is expected |230j to occur at 
T= (l.l-1.2)rc. 

In addition to high T breakup of cc or J/^, we have to take account of the so- 
called "normal suppression" of charmonium yields, observed in proton- nucleus 
collisions [233J. This effect is due to a re-scattering dissociation of the pri- 
mordially produced, pre-hadronic cc system upon traversal of (cold) hadronic 
matter [234j . It can be studied in p+A collisions where the data on J/^ pro- 
duction relative to pp collisions can be described by the survival probability 

SpA = -T^— = cfb dzpA{b, z)exp \ -{A - 1) / dz pA{b, z)crabs \ (6.4) 

where Gabs is the effective cross section for the "absorption" (break-up) of the cc 
in cold nuclear matter, and pA is the transverse nuclear density profile. The data 
[233 J suggest (Tabs = 4.2 mh. The generalization of equation 16.41 to the nucleus- 
nucleus case [235J gives a good description of the J/^ suppression (relative 
to binary pp scaling) in S-|-U and peripheral Pb-|-Pb collisions at top SPS 
energy [226j . It has thus become customary to quantify the J/^ suppression 
in central A-|-A collisions by relating the observed yield, not directly to the 
scaled pp yield (thus obtaining Raa)-, but to a hypothetical "normal absorption" 
yield baseline, established by equation 16.41 All further absorption is called 
"anomalous suppression" . 

Fig. O shows the results gathered aX ^ = 17.3 GeV by the NA38 - NA50 
- NA60 di-muon spectrometer [236], for minimum bias S+U {^/s = 20 GeV), 
In+In and Pb-|-Pb. Up to Npart ~ 100 all yields gather at the "normal absorp- 
tion" expectation from p-|-A scaling. A plateau at 0.8 follows for intermediate 
^part values up to about 200 (which corresponds to central In-|-In collisions, 
so the NA60 data end here), and a final falloff toward 0.6 for central Pb-|-Pb 
collisions. It appears natural to interpret the former step as an indication of ^ 
suppression, the final step as Xc dissociation. No genuine J/^ suppression is 
indicated. The expectation from lattice calculations |232j that J/^ dissociation 
does not occur until T ~ 2 Tc (and, thus, e ~ 16 Cc) is thus compatible with 
Fig. 16.11 We know from equation 16.31 that T < 1.3 Tc at top SPS energy. The 
data are thus compatible with no break-up of the J/^ at the SPS, unlike at 
top RHIC energy where one expects T ss 2Tc [HU [96] . 
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Figure 6.1: J/^ production measured in minimum bias collisions of S+U at 
Vs = 20 GeV and Pb+Pb and In+In at ./J = 17.3 GeV. The yield 
is scaled by " normal nuclear absorption" , equation 16.4! |236j . 



The RHIC data obtained by PHENIX [237J are shown in Fig. O Minimum 
bias Au+Au collisions are covered at mid-rapidity, as well as at 1.2 < y < 2.2, 
and plotted vs. Npart- Due to a parallel measurement of J/^ production in p+p 
collisions [238) the PHENIX experiment is in the position to show Raa-, but this 
is done without re-normalization to p-A absorption. J/^ is suppressed in cen- 
tral collisions, Raa ^ 0.2. Note that Raa can not drop down below about 0.15, 
due to unsuppressed surface contributions. The suppression is thus stronger 
than at top SPS energ}o - in fact it is almost maximal. We conclude that 
in central Au+Au at ^/s = 200 GeV the charmonium signal gets significantly 
quenched, in accord with the inferences about the primordial temperature that 
we presented in section [231 and HI to amount to about 300 MeV, i.e. T/T^ ~ 2 
as implied for J/^ dissociation by the lattice calculations }232j . 

The above interpretation is still a matter of controversy [239]. It remains un- 
clear whether successive, distinct stages of normal and anomalous J/^ suppres- 
sion are compatible with the dynamical evolution in central A+A collisions. A 
further open question refers to the difference in Fig. 16.2! of the (y) = and 
(y) = 1.7 data at intermediate centrality [239j interpretation of which should 

^Dropping the unfortunate distinction between normal and anomalous absorption one gets 
Raa = 0.35 for the central Pb+Pb collisions in Fig. 16.11 almost a factor 2 above the RHIC 
value. 
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be suspended until the advent of higher statistics data. 
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Figure 6.2: Raa for J/^ production in minimum bias Au+Au colhsions at ^/s - 
200 GeV, by PHENIX ^237] at mid-rapidity and at 1.2 < y < 2.2. 



A different aspect of charmonium production in A+A cohisions requires at- 
tention: at the stage of universal hadronization a certain fraction of bound cc 
mesons results, as the outcome of the density of uncorrelated c and c quarks 
in the QGP medium as it enters hadronization |240| I241J . The stage of sta- 
tistical hadronization (recall chapter E]) is omitted in the charmonium suppres- 
sion models [411 1230| 12341 1235] . which proceed in two steps (only): primordial 
nucleon-nucleon collisions produce an initial input of cc pairs, proportional to 
a^ X N^-jJ (b). In vacuum, the low relative momentum fraction of these cc 
pairs (about 1%) would evolve into charmonium hadrons, after a formation 
time of several fm/c. In the concurrently developing QGP medium, however, 
the initial cc correlated pairs may dissociate owing to the absence (due to color 
screening ^T\) of vacuum-like bound states |232j . at high T. At T > 2Tc all 
fireball charmonium production would thus cease, all cc pairs ending up in singly 
charmed hadrons. This picture 1230^ is incomplete. 



Even if all cc pairs break up at RHIC during the early phase of central collisions, 
the single charm and anti-charm quarks flow along in the expanding medium 
of equilibrated light and strange quarks. This picture is supported [M] by the 
observed elliptic flow of charm merging into the universal trend. Charm can 
not, however, be chemically (flavor) equilibrated at the low plasma tempera- 
ture, where it is essentially neither newly created nor annihilated |241| . The 
initially produced c and c quarks (a few per central Au+Au collision at RHIC, 
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a few tens at LHC) finally undergo statistical hadronization sX T = T^ along 
with all other light and strange quarks. To deal with the non-equilibrium over- 
all charm abundance an extra charm fugacity factor 7c is introduced into the 
statistical model [108] calculation (for details see [241]). J/\l' and ^ are thus 
created in non-perturbative hadronization, with multiplicities proportional to 
7c and phase space weights, along with all other charmed hadrons. This "re- 
generation" model also agrees with the RHIC data of Fig. 16.21 albeit within a 
large systematic uncertainty |241j . 

We note that the term regeneration is, in fact, misleading. The statistical 
hadronization process does not recover the initial, small fraction of correlated 
cc pairs that would end up in J /^ in vacuum. It arises from the total density of 
primordially produced c and c, uncorrelated in the hadronizing fireball volume. 

The statistical hadronization J/^ production process, sketched above, thus has 
the unfortunate property of providing a trivial background charmonium yield, 
unrelated to the deconfinement signal [41j referring to the primordial J/^ yield. 
Only about 1% of the primordial cc yield results in charmonia, in vacuum. The 
in-medium deconfinement process breaking up the cc correlation on its way 
to charmonia, thus constitutes a mere 1% fraction of the total charmed quark 
and anti-quark number. The regeneration process is insensitive to this 1% 
fraction, deconfined or not. At Tc, charm hadronization reacts only to the 
total abundance of c and c, as imprinted into the dynamical evolution by the 
perturbative QCD cc production rate of initial nucleon-nucleon collisions. At 
RHIC, it turns out [ 241J that the c and c density is low, giving rise to substantial 
canonical suppression (recalling equations 13. 13113. 171 in chapter [3]) of the two 
charm quark charmonia, relative to D mesons, during hadronization. With 
a tenfold c, c density at LHC, grand canonical charmonium production will 
set in, thus probably overshooting the primordial yield reference, cr^^ x Ncoii- 
Thus we expect Raa > 1 at the LHC. The role of a critical deconfinement 
"thermometer" is lost for J/^ at LHC, but the bottonium Y states can take 
over, being deconfined well above T = 300 MeV |242j . 

The RHIC result [257] for J/* in central Au+Au collisions (Fig. E2D, namely 
that Raa -^ 0.2, represents the lucky coincidence that the initial temperature, 
T « 300 MeV, is high enough to dissolve the correlated cc charmonium pre- 
cursor states, while the J/\I' suppression is not yet overshadowed by the trivial 
hadronization yield of J/^. 



6.2 Direct Photons 

Photons are produced during all stages of the dynamical evolution in A+A 
collisions. About 98% stem from final electromagnetic hadron decays, not of 
interest in the present context, other then by noting that their rate has to 
be painstakingly measured experimentally, in order to obtain "direct" photon 
spectra at low px by subtracting the decay fraction from the total. This was 
done [224j by WA98 in central Pb+Pb collisions at top SPS energy; we show 
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Figure 6.3: The WA98 direct photon transverse momentum spectrum for central 
Pb+Pb collisions at -y/i = 17.3 GeV . Also indicated are scaled pA 
results above 2.0 GeV/c and pQCD estimates |243j . From 



Only upper limits could be obtained at pT < 1-5 GeV/c due to overwhelming 
background from vr^ and r/ decay, thus also obscuring the major spectral domain 
in which to look for a direct photon QCD plasma black body radiation source, 
i.e. for thermal photons of a fireball at T between Tc and about 250 MeV |223j . 
Several data from p+p and p+A collisions at nearby ^/s and scaled up to central 
Pb+Pb are included in Fig. 16.31 at px > 2 GeV/c, and one sees the Pb+Pb 
data in clear excess of such contributions from primordial bremsstrahlung and 
hard, pQCD initial partonic collisions |243j . This excess ranges up to about 
3.5 GeV/c, in Fig. 16. 3i Above, the hard pQCD initial collision yield dominates 
|244j over thermal production. 

In contrast to all other primordial high pT pQCD yields (e.g. J/^, charm, jet 
leading partons) this photon yield is not attenuated in the medium of A+A 
collisions. We have shown in Fig. 15.41 the Raa for the PHENIX central Au+Au 
direct photon results |181j at RHIC ^/s = 200 GeV, obtained in a background 
substraction procedure |245j similar to the one undertaken by WA98. This 
procedure gives reliable data at pT > 4.0 GeV/c, at RHIC, and we observe 
Raa = 1- Hard initial photons are not attenuated, and there is no sign of any 
other direct photon contribution besides the primordial pQCD yield which, in 
fact, is shown (by Raa = 1) to obey binary scaling. However, there is no hint 
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to plasma thermal radiation (except for a trend at the very limit of statistical 
significance, at pT < 4.5 GeV/c) in this high px window. 

The WA98 SPS data, with thermal radiation enhancement indicated in the in- 
terval 1.5 < pt < 3.5 GeV/c, thus remained as the sole evidence until, more 
recently, the PHENIX experiment gained low px data [ 246] exploiting the fact 
that any source of real photons emits also virtual photons 7* leading to internal 
conversion to an e^e~ pair (the Dalitz effect). To identify this yield the invari- 
ant mass distribution of e^e~ pairs is analyzed outside the phase space limits of 
-K^ Dalitz decay; the decay pairs of all remaining hadron sources (??, A) is sub- 
tracted as a "cocktail". The remaining pair yield is then converted assuming 
^lirh inclusive = Idir- ll inclusive (sec ref. [216] for detail), thus finally obtaining 
data representative oi^dir in this approach. Fig. l6.4l shows the corresponding pj^ 
distribution which covers the interval 1.3 < pt ^ 4.5 GeV/c, within which the 
conventional direct photon extraction method did not give significant results 

[m]. 

The PHENIX experiment has also obtained direct photon spectra in p-|-p and 
d-FAu at ^ = 200 GeV pi7] which are both weh accounted for pi6] by a 
next to leading order (NLO) pQCD photon production model [248j . These 
data were already employed in deriving Raa = 1 for central Au-|-Au collisions, 
as shown in Fig. 15.41 and referred to, above. The pQCD fits derived from 
p-|-p and d-|-A are shown in Fig. 16.41 after binary scaling to Au-|-Au (pQCD 
xTaa)- They merge with the yield at pT ^ 4 GeV/c but demonstrate a large 
excess yield below 3 GeV jc. That excess is well described by adding a thermal 
photon component resulting from the hydrodynamic model of d'Enterria and 
Peressounko |249j . It traces the dynamical evolution during the early stages 
of equilibrium attainment, in which the photon luminosity of the emerging 
QGP matter phase is maximal. The hydrodynamic model provides for the 
space-time evolution of the local photon emission rate |223j which includes hard 
thermal loop diagrams to all orders, and Landau- Migdal- Pomeranchuk (LPM) 
in-medium interference effects. Similar, in outline, to previous models that 
combined algorithms of plasma photon radiation luminosity with hydrodynamic 
expansion }250| . the model |249] fits the data in Fig. 16.41 It implies a picture 
in which the early stage of approach toward thermal equilibrium at RHIC is 
governed by a symbolic, initial, effective "temperature" of about 550 MeV 
which, after equilibration at i ~ 0.6 fm/c, corresponds to T k, 360 MeV in the 
primordial plasma [249j: close to the consensus about initial T as derived from 
J/^ suppression, jet attenuation, elliptic flow and transverse energy coupled 
with the 1-dimensional Bjorken expansion model. 

However, we note that the employed theoretical framework, hard thermal loop 
(HTL) QCD perturbation theory of a weakly coupled plasma state, as combined 
with hydrodynamics, has a tendency to call for rather high initial T values. This 
applies both to the above analysis [249J of RHIC data which is based on the 
thermal field theory model of Arnold, Moore and Yaffe [223], and to previous 
analysis [250j of the WA98 SPS data of Fig. 16.31 Direct photon production 
is even more strongly biased toward the primordial, high T evolution than jet 
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Figure 6.4: Internal conversion measurement of direct photons in central 
Au+Au collisions at 200 GeV [MS]- Predictions by pQCD piH] 
and thermal hydrodynamic |249j models are included. 



attenuation (that may be proportional to T^ |212j ). Thus, by implication of 
the model [219] that fits the low pr RHIC data of Fig. MM the yield is highly 
sensitive to the (pre-equilibrium) formation period, 0.15 <t< 0.6 fm/c, where 
the HTL model might not be fully applicable. This illustrates the present 
state of the art. The model(s) based on perturbative QCD require extreme 
initial "temperatures" to produce the high photon yield, indicated by the RHIC 
experiment. The strongly coupled nature j213] of the non perturbative local 
equilibrium QGP state at RHIC, T ~ 300 MeV, mayprovide for an alternative 
approach to plasma photon production. 



6.3 Low mass dilepton spectra: vector mesons 
in-medium 

We have dealt with dilepton spectra throughout the above discussion of J/^ 
and direct photon production, as messenger processes sensitive to the energy 
density prevailing (during and) at the end of the primordial equilibration phase. 
The third tracer observable, low mass vector meson dilepton decay in-medium, 
samples - on the contrary - the conditions and modalities of hadron deconfine- 
ment in the vicinity of T = Tc. SPS energy is ideally suited for such studies 
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as the QGP fireball is prepared near Tc whereas, at RHIC, it races through 
the Tc domain with developed expansion flow. The major relevant data thus 
stem from the CERN SPS experiments NA45 ^25Tj and NA60 [136j . which have 
analyzed e~^e~ and ^"'"/i" production at invariant mass from 0.2 to 1.4 GeV. 

Fig. 16.51 shows NA45 data [251J for e^e^ production in central Pb+Au colli- 
sions at ^/s = 17.3 GeV. In searching for modifications of g properties and 
of vr+vr^ annihilation via virtual intermediate g decay to e+e~, in the high 
density environment near the hadron-parton coexistence line at T = T^ the 
various background sources have to be under firm control. The Dalitz decays 
of T^^,!] and T] and the in vacuo e^e~ decays of g,uj and <1>, which occur after 
hadronic freeze-out to on-shell particles, form a hadronic "cocktail" (left panel 
in Fig. 16. 5p that is generated from yields provided by the grand canonical sta- 
tistical model |108j . Within detector resolution, tTjUJ and $ leave distinct peaks 
but the observed invariant mass distribution is not accounted for. 

One needs also to account for background from Drell-Yan lepton pair produc- 
tion and open charm decay, both scaling with "number of collisions" A^'^. The 
latter contribution arises from primordial cc charm production, (c) = (c) which 
leads to synchronous formation of (D) = (^D^ at hadronization; subsequent de- 
cays D — > lepton + Ji, D ^ antilepton -|- Y create LL pairs. This procedure is 
straight forward as no significant medium attenuation occurs besides the statis- 
tical charm redistribution conditions at hadronization (section 16. ID , governing 
D,^ production |241j . 
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Figure 6.5: Di-electron mass spectrum for central Pb-|-Au collisions at y/s = 
17.3 GeV with the hadron decay cocktail (left) and further in 
medium contributions (right) [251j : see text for detail. 



Onward to non-trivial backgrounds in the invariant mass plot of Fig. 16. 5^ we 
recall the presence of thermal lepton pairs from virtual photon production in 
the early plasma phase, in parallel to real "direct" photon emission |252j . The 
spectrum of such pairs contains an average decay factor exp{ — y^}, with T the 
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(initial) plasma temperature. With T > 220 MeV assumed for top SPS energy 
|431 H3] , this contribution is a background candidate over the entire invariant 
mass interval covered by the data. In general Drell-Yan, open charm decay 
and plasma radiation contributions are smooth, partially closing the "holes" in 
the hadronic cocktail undershoot of the data. This smoothing effect is helped, 
finally, by consideration of modifications concerning the g meson spectral func- 
tion near T = Tc, which [1331 I134t I135J both affects the immediate g -^ e+e" 
decay invariant mass region (through the fraction of in-medium g decays vs. 
the in vacuum decay fraction after hadronic freeze-out) and, even more im- 
portantly, the contribution of in-medium tt~^tt~ annihilation to dileptons. The 
latter contribution accounts for the most obvious deviation between cocktail 
and data at 0.3 < m^/ < 0.7 GeV in Fig. [63] (left panel). 

The right hand panel of Fig. 16.51 shows the results of various theoretical models 
which address the sources of the significant dilepton excess over hadronic (in 
vacuum) and Drell-Yan cocktails, labeled " Rapp- Wambach" |133t 1135] . "drop- 
ping mass" (Brown-Rho |134j ) and " Kaempfer" |252| . We shall return to these 
models below but note, for now, that the extra yield in central A+A colli- 
sions chiefly derives from tt+tt" annihilation via the (in medium modified) g 
resonance, and from modification of the g peak itself. 

With improved statistics and background control, the A+A specific extra dilep- 
ton yield below M ~ 1.4 GeV/(? can be represented by itself, after cocktail 
subtraction. This has been first accomplished by NA60 |136t 1253) and, more 
recently, also by NA45 [254j . We show the former results in Fig. 16.61 The 
left panel shows the di-muon invariant mass spectrum in semi-central Indium- 
Indium collisions at top SPS energy y/s = 17.3 GeV , compared to the hadronic 
cocktail, qualitatively in agreement with the left hand panel of Fig. I6.5l but with 
superior resolution and statistics. The cocktail subtraction procedure (see |136j 
for details) leads to an invariant mass spectrum of di-muon excess in In+In, 
shown in the right side panel of Fig. 16.61 an experimental landmark accom- 
plishment. The g vacuum decay contribution to the hadronic cocktail has been 
retained and is shown (thin solid line) to be a small fraction of the excess mass 
spectrum, which exhibits a broad distribution (that is almost structureless if the 
cocktail g is also subtracted out), widening with collision centrality [136] . The 
best theoretical representation of the excess yield again results (like in Fig. 16.51 
right panel) from the broadening model [1331 11351 I255J where the g spectral 
function is smeared due to various coupling mechanisms within the medium 
via the vector dominance model, prior to hadronic freeze-out. In the VDM g 
couples, both, to pion pair annihilation, and to excited baryon states like the 
N* (1520), via their A^vrvr decay branches. 

The above data still imply serious conceptual questions. It is unclear to exactly 
which stage of the dynamical evolution (in the vicinity of T = Tc) the excess 
dilepton yield should correspond. As the observations appear to be coupled to 
in-medium g meson "metabolism" we need to identify a period of temporal ex- 
tension above the (in vacuo) g half life (of 1.3 fmjc)^ safely 2 fmjc. This period 
should be located in the vicinity of hadro- chemical freeze-out. From top SPS to 
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Figure 6.6: (left) Di-muon invariant mass spectrum in semi-central In+In col- 
lisions at ^/s = 17.3 GeV , with hadron final state decay cocktail; 
(right) excess mass spectrum after cocktail subtraction, confronted 
with fits from the broadening [VS'6\ I135t I255| and the dropping mass 
[228j models. From [136]. 



RHIC energy, hadro-chemical freeze-out should closely coincide with hadron for- 
mation, i.e. it should occur near the parton-hadron coexistence line, at T = Tc. 
The microscopic parton cascade model of reference [85] implements the Web- 
ber [121J phenomenological, non-perturbative QCD hadronization model (recall 
chapter 13. 3p which proposes pre-hadronization clusters of color neutralization 
(Fig. 13. 7p as the central hadronization step. In it, so one might speculate, the 
transition from pQCD to non perturbative QCD creates the chiral condensates 
(qq), spontaneously breaking chiral symmetry [256] (see below), and creating 
hadronic mass. The overall process, from pQCD color neutralization to on-shell 
hadrons, takes about 2.5 fm/c [851 E^ ^^ top SPS energy. This could, thus, 
be the period of excess dilepton yield creation. However, the relation of the 
models employed above [1331 1134^ 11351 12551 1256] to this primordial spontaneous 
creation of chiral condensates is still essentially unknown ^256j . 

Thus, at present, the 1990's paradigm of a direct observation of the chiral 
phase transition in QCD has been lost. The Brown-Rho model [134J predicted 
the Q mass to drop to zero at T = T^ occuring as a certain power of the ra- 
tio {qq)^^ I {qqY"''^ of the chiral condensate in medium and in vacuum which 
approaches zero at the chiral phase transition temperature, then expected to 
coincidence with the deconfinement temperature. This "dropping mass" model 
is ruled out by the data in Fig. 16.51 and |6.6[ This is, perhaps, a further mani- 
festation of the fact that the deconfined QGP state at T > Tc is not a simple 
pQCD gas of quarks and gluons [213] . In fact, lattice calculations [2321 1257] 
find indications of surviving light qq pair correlations in the vector channel at 
T > Tc- Thus the two most prominent symmetries of the QCD Lagrangian, 
non abelian gauge invariance (related to confinement) and chiral invariance 
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(related to mass) might exhibit different critical patterns at T = Tc and low 
baryo-chemical potential. This conjecture is best illustrated by the observation 
that the broad, structureless NA60 excess dilepton spectrum of Fig. 16.61 (after 
cocktail g subtraction) is equally well reproduced by a T ~ 160 — 170 MeV 
calculation in hadronic (equilibrium) matter \l'6'6\ 12531 1254| , and by a thermal 
QGP fireball of qq annihilation at this average temperature |252j . as illustrated 
here by the model curve labeled "Kaempfer" in Fig. 16.51 (right panel). This ob- 
servation has invited the concept of parton-hadron duality near T^ [258], which 
might be provocatively translated as "the QCD chiral transition properties can 
not be unambiguously disentangled from the deconfinement transition effects 
at Tc" [256]. 

We may be looking at the wrong domain of [T, hb] space: too high ^/s and, thus, 
too high T, too low ^b- Already at top SPS energy the medium is dominated 
by the deconfinement, not by the chiral QCD transition. After hadronization 
the medium is still at T close to Tc but the density drops off within a few fm/c, 
not allowing for an equilibrium mean field state of chirally restored hadrons. It 
is thus perhaps not surprising that the data are seen to be dominated by simple 
broadening and lack of structure: perhaps that is all that happens to hadrons 
at Tc- 

The chiral restoration transition should thus be studied at higher ^b and 
lower T such that the dynamics achieves high baryon densities but still merely 
touches the critical (deconfinement) temperature. In fact, at hb -^ 1 GeV and 
T < 100 MeV the chiral transition should be of first order [2 55] . Here, in fact, 
the chiral condensate mass plays the role of the order parameter (in analogy 
to the magnetization in a spin system), which approaches zero as T ^ Tc- We 
might thus speculate that, unlike at top SPS to LHC energy (where decon- 
finement dominates), the chiral QCD first order phase transition will dominate 
the phenomena occuring near the hadron-parton borderline, in the vicinity of 
-y/s = 4 — 6 GeV [142j. This requires a new experimental program, with low 
energy running at the SPS [260], at RHIC [Ml] and at the GSI FAIR project 
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Fluctuation and correlation signals in A+A collisions can be evaluated in single 
events, due to the high multiplicity of produced particles. Depending on the 
physics context we may be interested to see either a small, or a large nonstatis- 
tical fluctuation effect. For example in view of the universality of hadronization 
(chapter [3j) we would have difficulty with an event by event pion to baryon ra- 
tio (essentially ^^ ) fluctuating by, say, 50%. Conversely, searching for critical 
fluctuations in the vicinity of a predicted critical point of QCD [146^ I147J we 
would be frustrated if event- wise (pt), -^ (low px pion) |263j or strange to non- 
strange ratios like K/n [148^ I264J would not exhibit any significant fluctuation 
beyond statistics. Overall, event by event fluctuation observables also carry a 
message concerning the robustness of our assumptions about equilibrium at- 
tainment. It turns out that equilibrium properties are not, merely, central limit 
consequences of ensemble averaging. Rather to the contrary, each A ~ 200 cen- 
tral collision event at y/s > WGeV appears to offer, to the dynamical evolution 
of bulk properties, a sufficiently complete macroscopic limit. Such that we can 
view the event- wise bulk dynamics as " self analyzing" . 



-E 7.5 
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Figure 7.1: Initial transverse energy density distribution of a single central 
Au-FAu colhsion event at ^ = 200 GeV [265] . 
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7.1 Elliptic flow fluctuation 

That this ascertation is far from trivial is illustrated in Fig. 17. 1[ It shows 
|265j the primordial transverse energy density projection (the stage basic to 
all primordial abservables) at i ~ 0.2 fm/c, in a central Au+Au collision at 
^/s = 200GeV , exhibiting an extremely clumpy and nonhomogeneous structure, 
apparently far away from equilibrium. 

The imaging of this initial geometrical pattern of the collision volume by a 
hydrodynamic evolution even applies at the level of individual events, such as 
illustrated in Fig. EH The PHOBOS Collaboration has shown [Ml [2671 EM] 
that an event by event analysis of the elliptic flow coefficient V2 is possible (see 
ref. |266j for detail), by means of a maximum likelyhood method. For initial- 
ization they sample the seemingly random initial spatial density distribution 
of single Au+Au collision events by the " participant excentricity" of individual 
Monte Carlo events, 



^\ 



2 - ^2^2 



+ 4ct2 , 



y ^xl ' ^" xy 
tpart = 2~~i 2 \'-^) 

a^ + ai 



where axy = {xy) — (x) (y). The average values of epart turn out to be sim- 
ilar to €x from equation 14. H as expected, but the relative fluctuation width 
(T(e)/ {e)part turns out to be considerable. It is the point of this investigation 
[2671 I268J to show that the observed relative event by event flow fluctuation 
equals the magnitude of the relative excentricity fluctuation. This is shown 
|268j in Fig. 17.21 The left panel demonstrates that the average (^2) obtained 
vs. Npart from the event-wise analysis agrees with the previously published 
|156j event-averaged PHOBOS data. The right panel shows that the event-wise 
relative fluctuation of V2 is large: it amounts to about 0.45 and is equal to the 
relative fluctuation of epart, i-e. 



o-(^2)/ (-"2) ~ Cr{epart)/ (epart) ■ (7.2) 



The initial geometry appears to drive the hydrodynamic evolution of the system, 
not only on average but event-by-event [268j . thus providing for an example of 
the self-analyzing property mentioned above. The f 2 signal thus emerges as the 
most sensitive and specific diagnostic instrument for the primordial conditions 
and their subsequent evolution: it reveals even the random (Fig. 17. ip initial 
fluctuations. In comparison the analysis with thermal photons is only sensitive 
to the primordial temperature j249j , and restricted by very small cross sections 
and significant background from other sources and evolution times. It also does 
not give viscosity information. 
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Figure 7.2: (left) Event by event elliptic flow analysis by PHOBOS gives an 
average (^2) that agrees with the result of ensemble analysis, in 
Au+Au at 200 GeV, for charged hadrons. (right) The event-wise 
relative V2 fluctuation vs. Npart, compared to the event-wise rela- 
tive fluctuation of the participant excentricity, a{epart)/ impart)- The 
closed line gives V2 variation due to Npart number fluctuation. From 



7.2 Critical point: fluctuations from diverging 
susceptibilities 



Recalling the goal defined in the introduction we seek observables that help to 
elaborate points or regions of the QCD phase diagram, Fig. 11.11 We have seen 
several observables that refer to the QCD plasma environment at T > 300MeV, 
H ^ (elliptic flow, jet attenuation, J/^ suppression, direct photon pro- 
duction), which could be first resolved within the primordial time interval 
r < 1 fm/c accessible at RHIC energy. The LHC will extend the reach of 
such observables toward T « 600 MeV, xp < 10~^, at fis = 0. On the other 
hand, relativistic A+A collisions at lower energy permit a focus on the hypo- 
thetical QCD parton-hadron coexistence line, T = Tc{^b), with the domain 
IXb -^ 500 MeV being accessible at the SPS. Characteristic observables are ra- 
dial flow, hadro-chemical freeze-out, and chiral symmetry restoration effects in 
dilepton vector meson spectra. Focusing on this domain, we discuss fluctuations 
potentially associated with the existence of a critical point [8l[9t[T0 t[Tnil46iri47j . 



At the end of chapter 16.31 we mentioned the conclusion from chiral symmetry 
restoration models |15[ I259J that at high fis the phase transformation occuring 
at Tc {^j-b) should be a chiral first order phase transition. On the other hand, 
lattice QCD has characterized [16J the phase transformation at ^b -^ 0, to be 
merely a rapid cross-over. Thus, the first order nature of the phase coexistence 
line in Fig. 11.11 has to end, with decreasing ^b, in a QCD critical point, ten- 
tatively located by recent lattice QCD calculations [U [lOl [H] in the interval 
fj-B = 300 — 500 MeV. The existence of such a point in the [T, /ib] plane would 
imply fluctuations analogous to critical opalescence in QED [1461 I147t l263j . Be- 
yond this second order phase transition point the coexistence line would be the 
site of a rapid cross-over p^. This overall theoretical proposal places poten- 
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tial observations related to the critical point itself, and/or to the onset of first 
order phase transition conditions at higher fiB , within the domain of the lower 
SPS energies, ^/s < 10 GeV. Note that, at such low energies, the initialization 
of thermal equilibrium conditions should occur in the vicinity of Tc, unlike at 
RHIC and LHC, and that the central fireball spends considerable time near the 
coexistence line, at 300 < /Ub < 500 MeV. 
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Figure 7.3: Quark number density susceptibility vs. temperature for light 
quarks in 2 flavor lattice QCD at flnite fiB- The calculation refers to 
Tc = IbOMeV and quark chemical potential fiq/Tc = 0, 0.5 and 1.0, 
respectively |270j . Smooth lines interpolate the calculated points; 
error bars indicate lattice statistics. 



To analyze potential observable effects of a critical point, we recall briefly the 
procedure in flnite fiB lattice theory that led to its discovery. One method to 
compute thermodynamic functions at /ib > from the grand canonical parti- 
tion function Z{V, T, Hg) at /Xg = is to employ a Taylor expansion with respect 
to the chemical quark potential [TOl [HI EZO] , defined by the derivatives of Z at 
/x = 0. Of particular interest is the quark number density susceptibility, 

Xu4 = T^ i^r^ ^^ (7.3) 



which can also be written as 



Xq 



+ 



Uu + Ud 



(7.4) 



with Xq = iXu + Xd)/2 and quark number densities n„, n^. We see that the 
susceptibility refers to the quark number density fiuctuation. The lattice result 



116 



7.2 Critical point: Euctuations from diverging susceptihilities 



[101 I270J is shown in Fig. 17.31 a calculation with two dynamical flavors assuming 
Tc = 150MeV and three choices of chemical quark potential, fiq = 0, 75 and 150 
MeV, respectively, corresponding to fiB = 3 //^ = 0, 225 and 450 MeV. These 
choices correspond to LHC/RHIC energy, top SPS energy and -y/i ss 6.5 GeV, 
respectively. At fiB = one sees a typical smooth cross-over transition at 
T = Tc whereas a steep maximum of susceptibility occurs with ^b = 4:50 M eV . 
This suggests the presence of a critical point in the (T, hb) plane [270] in 
the vicinity of (150 MeV, 450 MeV). For final confirmation one would like 
to see this maximum disappear again, toward fig > Tc, but this is beyond 
the convergence domain of the employed Taylor expansion (see ref. [9] for 
alternative approaches). 
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Figure 7.4: Event by event fluctuation of the K~ /tt^ ratio in central collisions 
of Pb+Pb at ^/s = 17.3 GeV, relative to mixed event background 
(histogram) [272j . 



From Fig. 17.31 we also expect a divergence of the strangeness susceptibility, for 
which no results from a 3 dynamical flavors calculation at finite fiB exist to date. 
A lattice calculation at /i^ = 0, T ps 1.5 Tc suggests |271j that the u, d, s quark 
flavors densities fluctuate uncorrelated (but we do not know whether that is also 
true at fj,B = M^**). This could thus be observed in event by event analysis, in 
particular as a fluctuation of the Wroblewski ratio A^ = 2(s + s)/(u + u + d + d) 
which is approximated by the event- wise ratio (K^ + K^)/{7r^ + 'it~). This was 
first measured by NA49 in central collisions of Pb+Pb at top SPS energy; the 
result ^272j is shown in Fig. 17.41 The data result from a maximum likelyhood 
analysis of track- wise specific ionization in the domain 3.5 < y < 5 slightly 
forward of mid-rapidity. The width adata is almost perfectly reproduced by the 
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mixed event reference, such that the difference, 



^dyn 



a 



data 



cr„ 



(7.5) 



amounts to about 3% of (Jdata only, at -y/i = 17.3 GeV. This analysis has 
more recently been extended to all energies available thus far, at the SPS [273] 
and at RHIC |274j . Fig. 17.51 shows that adyn stays constant from top SPS to 
top RHIC energy but exhibits a steep rise toward lower energies that persists 
down to the lowest SPS energy, ^/s = 6.2 GeV. Fig. 13.91 shows |107j that at 
this energy f^s = 450 MeV, thus corresponding to the susceptibility peak in 
Fig. 17.31 Again, as we noted about the peak in Fig. 17.31 if these data indicate 
a critical point effect in the vicinity of fis = 450 MeV the relative fluctuation 
should decrease again, toward yet higher hb and lower ^/s. These data will, 
hence, be re-measured and extended to lower ^/s by experiments in preparation 
[2601 ESU |262] at CERN, RHIC and GSI-FAIR. This will also help to evaluate 
alternative tentative explanations invoking fluctuating canonical suppression 
[120j . or strangeness trapping |275j . Finally, the position of the critical point 
needs to be ascertained by lattice theory. 
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Figure 7.5: The relative deviation of the event by event K/tt fluctuation width 
from the mixed event background width, adyn (see equation 17. 5p : at 
SPS [573] and RHIC [271] energies. 



7.3 Critical fluctuation of the sigma-field, and related 
pionic observables 

Earlier investigations of critical QCD phenomena that might occur in high 
energy nuclear collisions were based on QCD chiral field theory |276j . The QCD 
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critical point is associated with the chiral phase transition in so far as it appears 
as a remnant of a tri-critical point |147j corresponding to the " ideal" chiral limit 
that would occur if rriu = rud = 0. Therefore the existence of a second-order 
critical point, at /i^ > 0, is a fundamental property of QCD with small but 
non-zero quark masses j277| . The magnitude of the quark condensate, which 
plays the role of an order parameter of the spontaneously broken symmetry 
(generating hadronic mass), has the thermal expectation value 

im) = ^'^{n\qq\n) exp[-En/T) (7.6) 

n 

with the partition function of hadronic states En 

Z = Y,exp{-EjT). (7.7) 

n 

The low energy behaviour of the matrix elements {n\qq\ n) can be worked 
out in chiral perturbation theory |277j . At the QCD critical point the order 
parameter fluctuates strongly. Its magnitude (gg) is identified with an isoscalar 
quantity, the so-called a-field. The critical point communicates to the hadronic 
population via the a <-^ tttt reaction, generating fluctuating fractions of the 
direct pion yield present near T = Tc, which thus gets imprinted with a fluctu- 
ation of transverse momentum (in the low pT domain) stemming from a mass 
fluctuation, downward toward the critical point. At it the isoscalar field ide- 
ally approaches zero mass, in order to provide for the long wavelength mode 
required by the divergence of the correlation length |147j . 

Note the relatively fragile structure of the argument. In an ideal, stationary in- 
finite volume situation the sigma field would really become massless, or at least 
fall below the TT~^7r~ threshold; thus its coupling to tt~^tt~ becomes weak, and 
restricted to very small px- Furthermore, such primary soft pions, already small 
in number, are subject to intense subsequent re-absorption and re-scattering in 
the final hadronic cascade evolution |114j . In fact, experimental investigations 
of event by event pT fluctuations in central A+A collisions, covering the en- 
tire ^/s domain from low SPS to top RHIC energy have not found significant 
dynamic effects [Ml EZS IMl ESI]- Fig. ES illustrates the first such mea- 
surement by NA49 [278j in central Pb-|-Pb collisions at y/s = 17.3 GeV, at 
forward rapidity 4 < y < 5.5, showing the distribution of event-wise charged 
particle average transverse momentum, a perfect Gaussian. It is very closely 
approximated by the mixed event distribution, ruling out a significant value of 
adyn from equation 17.51 More sensitive measures of changes, event by event, 
in the parent distribution in transverse momentum space, have been developed 
[2821 EM]. NA49 has employed [STHJEZS] the measure $(pt), defined as ^ 



where Zi = pTi — Pj- for each particle, with pj. the overall inclusive average, 
and for each event Z = ^^ Zj is calculated. With the second term the trivial 
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independent particle emission fluctuation is subtracted out, i.e. ^ vanishes 
if this is all. Indeed, the data of Fig. 17.61 lead to <1> compatible with zero. 
Furthermore, a recent NA49 study j279j at mid-rapidity, covers the range from 
^/s = 17.3 to 6.3 GeV (where the K/n ratio fluctuation in Fig. 17.51 exhibits 
the much-discussed rise, and even the ensemble average in Fig. 13.41 shows the 
unexplained sharp peak) but finds no significant ^ signal. 
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Figure 7.6: Event by event fluctuation of average charged hadron pT in the in- 
terval 4.0 < y < 5.5, in central Pb+Pb collisions at ^/s = 17.3GeV. 
Mixed event background given by histogram |278] . 



Alternatively, one can base a signal of dynamical pT fluctuation on the binary 
correlation of particle transverse momenta in a given event, i.e. on the co- 
variance {pTiPTj) [281^ I283J of particles i,j in one event. Of course, the co- 
variance receives contributions from sources beyond our present concern, i.e. 
Bose-Einstein correlation, flow and jets (the jet activity becomes prominent 
at high ^/s, and will dominate the pT fluctuation signal at the LHC). In co- 
variance analysis, the dynamical pr fluctuation (of whatever origin) is recovered 
via its effect on correlations among the transverse momentum of particles. Such 
correlations can be quantified employing the two-particle pT correlator [28 1^ 



{ApTi ApTj) 



1 



„ Nik) N{k) 



M, 



pairs 



E E E ^PT^^PTj 



(7.9) 



fc=l i=l j=i+l 



where Mpairs is the total number of track pairs of the events k contained in 
the entire ensemble of n events, N{k) is the number of tracks in event k, and 
ApTi = PTi — Pt where p-p is the global ensemble mean px- The normalized 



120 



7.4 Bose-Einstein-correlation 



dynamical fluctuation is then expressed |281| as 



(y{pT)dyn = Ji^PTi ^PTj) / PT- 



(7.10) 



It is zero for uncorrelated particle emission. 



Fig. 17.71 shows the analysis of pT fluctuations based on the px correlator, for 
central Pb+Au SPS collisions by CERES [280J and for central Au+Au at four 
RHIC energies by STAR [281] . The signal is at the 1% level at all ^/s, with 
no hint at critical point phenomena. Its small but finite size could arise from 
a multitude of sources, e.g. Bose-Einstein statistics. Coulomb or flow effects, 
mini-jet-formation, but also from experimental conditions such as two-track 
resolution limits |284j . We note that even if a critical opalescence effect, related 
to a fluctuating chiral condensate at T = Tcrit, couples to the primordial, low 
Pt pion pair population ^1471 1285] . this signal might be dissipated away, and 
thus " thermalized" to the thermal freeze-out scale of about 90-110 MeV, as a 
pion experiences about 6 re-scatterings during the hadronic cascade [114] . On 
the other hand the hadro-chemical K/tt ratio fluctuation (Fig. 17. 5p would be 
preserved throughout the cascade (chapter [3]) . 
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Figure 7.7: Dynamical px event by event fluctuation analysis by cr{pT)dyn of 
equation 17101 vs. ^/s, showing SPS [280] and RHIC [28l] data. 



7.4 Bose-Einstein-correlation 

Identical boson pairs exhibit a positive correlation enhancement when Apij -^ 0, 
an intensity correlation analogous to the historical Hanbury-Brown and Twiss 
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effect (HBT) of two photon interferoinetry |286j employed in astrophysics. In 
nucleus-nucleus collisions this is an aspect of the symmetry of the N-pion 
wave function that describes pion pairs at the instant of their decoupling from 
strong interaction. The momentum space correlation function C{q,K), q = 
Pi—p2, K = ^{pi +P2) is the Fourier transform of the spatial emission function 
S{x, K) which can be viewed as the probability that a meson pair with mo- 
mentum K is emitted from the space-time point x in the freezing-out fireball 
density distribution [9611287] . 

The aim of HBT two particle interferometry is to extract from the measured 
correlator C(g, K) as much information about S{x, K) as possible. In addition 
to the traditional HBT focus on geometrical properties of the source distribution 
("HBT radii") there occurs information on time and duration of emission, as 
we have already employed in Fig. 12.181 illustrating tj and At in central Pb-|-Pb 
collisions at top SPS energy [90j. Moreover, even more of dynamical informa- 
tion is accessible via the effects of collective flow on the emission function. In 
this respect, the observations make close contact to the hydrodynamic model 
freeze-out hyper-surface (chapter 13. ip and the contained flow profiles. HBT 
thus helps to visualize the end of the collective radial and elliptic flow evolu- 
tion. This implies that we expect to gather evidence for the, perhaps, most 
consequential property of A-l-A collisions at high ^/s, namely the primordially 
imprinted Hubble expansion. 

We do not engage here in a detailed exposition of HBT formalism and results 
as comprehensive recent reviews are available [96^ ,287] . Briefly, the measured 
quantity in 7r"'"7r''', tt^vt^ or K^K'^ interferometry is the correlator 

r(, K^ <fN/dpldpl .7^^^ 

^^'^'^^- d^N/dp\d^N/dpl ^^-^^^ 

of pair yield normalized to the inclusive yield product. The correlator is related 
to the emission function which is interpreted as the Wigner phase space density 
of the emitting source |288j : 

I fd^xS(x,K)e''^^ |2 

^'''•^'"^+ ifrf%5(;.K)p ■ '"'' 

Due to the experimental on-shell requirement Kq = \JK^ -|- w? the 4-vector 
components of K are constrained on the left hand side. Thus, relation l7. 121 can 
not simply be inverted. 

To proceed one employs a Gaussian ansatz on either side of equation 17.121 The 
experimental data are parametrized by source "radius" parameters Rij{K), 

C{q,K) = 1 + X{K)exp [-Y,R%{K)mj ] (7.13) 



employing X{K) essentially as a fudge factor (related nominally to possible 
coherence of emission effects in the source, which would dilute or remove the 
Bose-Einstein statistics effect of C — > 2 for q ^ Q but have never been seen 
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in nuclear collisions). In equation 17.131 the sum runs over three of the four 
components of q, due again to the on-shell requirements [287] . For the emis- 
sion function S{x,K) a Gaussian profile is assumed about an "effective source 
center" x{K), thus 

S{x,K)^ S{x{K),K)xG (7.14) 

where G is a Gaussian in coordinates x^^(K) relative to the center coordinates 
x^{K). Inserting equation 17.141 into 17.121 one finally obtains 

C{q, K) = l + \{K) exp [-q^q, (x'^ i")] (7.15) 

where {x^ x'^) are the elements of the space-time variance of the correlation 
function, which re- interpret the "radii" Rf: in equation I7.13[ Assuming az- 
imuthal symmetry (central collisions), cartesian parametrizations of the pair 
relative momentum q coordinates (corresponding to fixation of the space-time 
variance in equation 17. 13p have been introduced by Yano, Koonin and Podgoret- 
skii [289] , and, alternatively, by Pratt |290j . The latter, out- side-longitudinal 
coordinate system has the "long" direction along the beam axis. In the trans- 
verse plane, the "out" direction is chosen parallel to Kt = (pit +P2r)/2, the 
transverse component of the pair momentum K. The "side" direction is then 
orthogonal to the out- and long-direction but, moreover, it has the simplest 
geometrical interpretation (see ref. |287j for detail), to essentially reflect the 
transverse system size |288j . The parameters of equation 17.151 are thus defined; 
as an example we quote, from identification of equation 17.131 with I7.15t the 
resulting geometrical definition of the "side" radius, 

Rl,,{K) = (yiKf) ■ (7.16) 

Overall, this model of Fourier related correlators, C{q, K) the experimentally 
accessible quantity (see equation 17. lip , and S{x.,K) the to-be-inferred spatial 
freeze-out fireball configuration, leads to the Gaussian ansatz }287j 



C{q,K) = 1 + \{K) exp[-Rl^,iK)ql, - Rlae{K)ql 



Rfong {K)qJong + cross terms] (7.17) 

which is fitted to the experimental correlation function (equation 17. lip . The 
experiment thus determines the variances Rout-, Rside^^'^^Riong- In the so-called 
"local co-moving system" (LCMS), defined as the frame in which p^^i = —Pz,2, 
i.e. (3iong = 0, we obtain in addition to equation 17.161 



RaatiKf = {{x{K) - (3Ti{K)f 

Hong 



RLJK) = (~z{Kf) (7.18) 



and finally, for azimuthal symmetry in central collisions with (5;^) ~ (y^) we 
find the " duration of emission" parameter illustrated in Fig. 12.181 [90j : 

{P)^^{Rl^^-Rl.^^). (7.19) 

The resulting reduction of the initial 8 coordinates of the meson pair stems, in 
summary, from the on-shell requirement, from azimuthal symmetry and from 
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approximate Bjorken invariance in the LCMS system |287l I288J . One has to 
be aware of the latter ideahzations. Note that in an expanding source all HBT 
parameters depend on K, the pair mean momentum (see below). 

We make first use of the above parametrization in Fig. 17.81 From the purely 
spatial radii Rside ^-nd Riong one can define a mid-rapidity volume at pionic 
decoupling, Vj = (27r)2/3 i?^.^^ ^^^^^ ^i^^^i^ jg shown ^mj at K = 0.2 GeV 
for central Pb+Pb and Au+Au colhsions from AGS [M! via SPS ^ [MS] 
to RHIC [294] energy. The upper panel shows the ^/s dependence, the lower 
illustrates the dependence on the charged particle rapidity density dNch/dy 
which one might intuitively expect to be related to the freeze-out volume [295j . 
We see, firstly, that the plot vs. y/s exhibits a non-monotonous pattern at the 
transition from AGS to SPS energies [287[ 1295] . whereas the plot vs. dN/dy 
rather features a rise toward a plateau that ends in a steep increase at RHIC 
energies. Second, the tenfold increase in charged particle rapidity density is 
reflected in only a doubling of the "volume" Vf. 
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Figure 7.8: Coherence freeze-out volume Vf from 7r~pair Bose-Einstein correla- 
tion analysis in central Au+Au and Pb+Pb collisions, (upper panel) 
plotted vs. ^/s, (lower panel) vs. mid-rapidity charged particle den- 
sity dN^^/dy [2M] . 



The latter observation reminds us of the fact that the geometrical parameters 
obtained via the above analysis do not refer to the global source volume if that 
volume undergoes a collective Hubble expansion [96l I287J . A pion pair emitted 
with small relative momentum into the azimuthal direction Kt is likely to 
stem (only) from the fraction of the expanding source that also moves into 
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this direction. This coupling between position and nionientuni in the source 
becomes more pronounced, both, with increasing Kt and increasing sources 
transverse velocity I3t from radial expansion. We have seen in Fig. 12.231 that 
the latter increases dramatically with ^/s, such that the coherence volume Vf 
comprises a decreasing fraction of the total fireball. It should thus rise much 
more slowly than proportional to the global dN/dy |96 tl287j . 

A striking experimental confirmation of the Hubble expansion pattern in central 
A+A collisions is shown in Fig. 17.91 The illustrated HBT analysis of NA49 
[90] at Vs = 17.3 GeV, and of PHOBOS [29i] at ^ = 200 GeV , employs 
the alternative parametrization of the correlation function C(g, K) introduced 
|289j by Yano, Koonin and Podgoretskii (YKP). Without describing the detail 
we note that the YKP correlation function contains the "YK velocity" I3yk 
describing the source's longitudinal collective motion in each interval of pion 
pair rapidity. 

Defining the "YKP rapidity" corresponding to (Syk by 

Yykp = 2 ^'^TZTr ^ ^'^"^ i^-'^^) 

leads to the results in Fig. 17. 9^ indicating a strong correlation of the collective 
longitudinal source rapidity Yykp with the position of the pion pair in rapidity 
space. 

We are interested in a possible non-monotonous ^/s dependence of freeze-out 
parameters such as Vf because of recent predictions of the relativistic hydrody- 
namic model that a QCD critical point should act as an attractor of the isen- 
tropic expansion trajectories in the [T, ^b] plane [lOlj 1296] . Fig. 17.101 shows 
such trajectories characterized by the ratio of entropy to net baryon number, 
S/ub, which is conserved in each hydro- fluid cell throughout the expansion. 
Note that the S/ub ratio is determined from the EOS during the primordial 
initialization stage in the partonic phase; the relationship between S/ub and ^/s 
is thus not affected by the presence or absence of a critical end point (CEP) at 
T = Tc which, however, drastically influences the trajectory at later times as is 
obvious from comparing the left and right panels, the latter obtained with a first 
order phase transition at all fiB but no CEP. In this case, the S/ub = const, 
trajectories in the partonic phase all point to the origin in the [T, fig] plane 
because //g/T oc In(^) in an ideal parton gas; whereas they point oppositely 
below T = Tc because T^''^/q oc ln(-^) in a hadron gas. This organization is 
dramatically upset for the cases S'/n=100, 50, and 33 by the assumption of a 
CEP, tentatively placed here at T = 155 MeV, ^b = 368 MeV. 

A first conclusion from this model is that it should not be essential to fine-tune 
^/s to make the system pass near the CEP because the attractor character 
would dominate over a substantial domain, e.g. at hadro-chemical freeze-out 
conditions 250 < fiB < 500 MeV in the example of Fig. 17.101 Please note that 
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Figure 7.9: Emitting source rapidity (Yykp) as a function of pion pair rapidity 
(I^tt)- From 7r~pair correlation analysis in central Pb+Pb collisions 
at ^/s = 17.3 GeV ^90j, and in central Au+Au collisions at ^/s = 
200 GeV [291]. 
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Figure 7.10: (left) Influence of a critical point on hydrodynamic model isen- 
tropic expansion trajectories characterized by various values of en- 
tropy to net baryon number s/ub- (right) The same trajectories 
without a critical point but a first order transition all along the 
phase boundary 



hadro-chemical freeze-out is not treated correctly in this model, to occur at Th 
which is 160 < Th < 130 MeV from Figs. 11.11 and 13.91 The trajectories shown 
here below Th are thus not very realistic. The expected pattern could cause 
the "plateau" behaviour of several observables at hadronic freeze-out over the 
range of SPS energies that corresponds to the above interval of ^b and Th, e.g. 
{ruT) and T in Figs. 12.201 and 12.211 elliptic flow V2 at mid-rapidity in Fig. 14.31 
and coherent hadronic freeze-out volume Vf from HBT in Fig. | 
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A consistent search for critical point effects implies, at first, a correct treat- 
ment of the hadronic expansion phase in hydro- models as above |10H I296J . 
properly distinguishing chemical composition freeze-out, and eventual "ther- 
mal" freeze-out at hadronic decoupling. From such a model predictions for the 
^/s or S/ns systematics could be provided for the HBT source parametrization 
implied by equations I7.13ll7.15l I.e. the hydro dynamic model provides the cor- 
relator S{x, K) in cases with, and without a critical point which, as Fig. I7.1UI 
shows, leads to considerable changes in the system trajectory in the domain 
near hadronization and onward to hadronic thermal freeze-out, at each given 
S/n or y/s. On the experimental side, this trajectory is documented by [T, ^b\ 
at hadronic freeze-out from grand canonical analysis |140j (chapter [3]) which 
also yields S/ub- Furthermore, as we shall show below, HBT in combination 
with the analysis of pT or m^ spectra will describe the situation at thermal 
hadron freeze-out yielding Tf, (5t at the surface and the "true" transverse ra- 
dius Rgeom of the source, in addition to the coherence volume Vj illustrated 
in Fig. 17.81 (which also documents the present, insufficient data situation). Of 
course, we have to note that the lattice results concerning the critical point are 
not yet final |27U| , and neither is the hydrodynamic treatment p96j concerning 
the EOS in the vicinity of the CEP. 



200 




Figure 7.11: Allowed regions of freeze-out temperature vs. radial expansion ve- 
locity for central Pb+Pb collisions at ^/s = 17.3 GeV and mid- 
rapidity, combining negative hadron and deuterium spectral data 
analysis with BE vr^ correlation results on Kt dependence of 
Ra_ ^ Rs [90]. 



We turn to combined analysis of the two final processes of expansive evolu- 
tion, formation of pT spectra and decoupling to Bose-Einstein pair correlation, 
in order to show how the experimental information mentioned above can be 
gathered. At the level of " hydrodynamically inspired" analytic models for the 
emission function S{x, K) several approaches [1061 I297J have established the 
connection, arising from radial flow, between hadronic rriT spectra and the Kt 
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dependence of the HBT radii Rside, Rout and Riong, equation 17. 16^ I7.2H which 
fall down steeply with increasing Kt |29H 12921 12931 I294J . A combined analy- 
sis lifts continuous ambiguities between the thermal freeze-out temperature Tf 
and the radial expansion velocity /3t (for which a radial profile Pt = w' — /^o is 
assumed), that exist in the blast wave model derivation of both the pT spectra, 
equation 12.251 and the Kt dependence of Rs- 

This was first demonstrated in a NA49 study of pion correlation in central 
Pb-fPb collisions at ^ = 17.3 GeV [90] that is shown in Fig. \7JT[ The 
ambiguous anti-correlation of fit parameters Tf and /3y can be somewhat con- 
strained if several hadronic species are fit concurrently. This is obvious in 
Fig. 17.111 from the overlap of parametrization regions for negative hadron, and 
for deuterium nirp spectra. An orthogonal constraint on P^/Tj results from a 
blast wave model fit of the HBT Kt dependence of Rt ~ Rside [ 298j . employed 
here. We see that the three independent 1 a fit domains pin down Tf and 
Pt rather sharply, to [115 MeV,0.55]. A relativistic correction |297j leads to 
Tf = 105 MeV, Pt = 0.60. The ^/s dependence of Pt at the freeze-out surface, 
from such an analysis [2911 I297J is shown in Fig. 17.121 The data again exhibit a 
plateau at SPS energies, which remains to be understood [299J. In the light of 
the considerations above, the plateau might turn out to reflect a critical point 
focussing effect. 
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We have seen that many physics observables in relativistic nucleus-nucleus col- 
lisions can, at RHIC energy y/s = 200Gey, be related to the primordial dy- 
namical phase, from initial QCD parton shower formation to local momentum 
space equilibration. The time interval is 0.35 to 0.65 fm/c. This domain can 
be investigated at RHIC due to the short interpenetration time, of 0.15fm/c. 
From among the bulk hadron production signals, total and midrapidity charged 
particle rapidity densities, dNch/dy, reflect the primordial parton saturation 
phenomenon [72]. It leads to an unexpectedly slow increase of multiplicity with 
fireball participant number Npart, and with y^. This observation signals the 
onset of nonperturbative QCD, a coherent shower multiplication by multigluon 
coherence (HU EH [701 EH ES]- It is expected to be even more dominant in 
LHC Pb+Pb collisions at 5.5TeV. Furthermore, elliptic flow, a collective bulk 
hadron emission anisotropy, also originates from the primordial, nonisotropic 
spatial density profile of shower-produced partons [93l |96] . A hydrodynamic 
evolution sets in at t < Ifm/c, implying the existence of a primordial equa- 
tion of state (EOS) for partonic matter in local equilibrium. Moreover, the 
experimental elliptic flow data at RHIC are well described by "ideal fluid" hy- 
drodynamics, i.e. by a very small shear viscosity ry and a small mean free path 
A [Ml [95]. 

These observations indicate the existence of a (local) equilibrium state at very 
early times in the bulk parton matter stage of the dynamical evolution. This 
quark-gluon plasma (QGP) state of nonperturbative QCD was predicted by 
QCD lattice theory [IT]. In Au+Au collisions at RHIC energy this state appears 
to be realized under the following conditions |61| . 

The energy density e amounts to 6 — 15GeV/fm'^, far above the parton-hadron 
QCD confinement region at e = IGeV/fm^, and at about 55 times the density 
of nuclear ground state matter, po = 0.14GeV/fm^. Translating to macro- 
scopic units we are dealing with a matter density of about 1.3 • 10^'^ kg /m^, 
the density prevailing in the picosecond era of the cosmological evolution. The 
corresponding temperature amounts to T=300-330 MeV (about 3.6 • 10^^ de- 
grees Kelvin), far above the "Hagedorn limit" ^38j for the hadronic phase 
(Th = 165MeV). From an analysis of the production ratios of the various 
hadronic species (from pions to Omega hyperons) arising from this primordial 
partonic state at hadronization, the statistical hadronization model (SHM) de- 
termines its baryo-chemical potential |108j , hb = 20MeV at RHIC energy. This 
value indicates that one is already close to the near-zero net baryon number 
conditions, hb ~ 0, prevailing in the corresponding big bang evolution, where 
the density of particles exceeds that of antiparticles by a fraction of about 10~^ 
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only. 

Overall we thus obtain entries in the QCD phase diagram of Fig. II. 1[ RHIC 
creates a parton plasma at about T = SOOMeV and fiB = 20MeV. It turns 
out to behave like an ideal fluid and features an extremely short mean free path 
A, thus inviting a description by relativistic hydrodynamics. The small shear 
viscosity 77 (or the small viscosity to entropy ratio r]/s) are highlighted by the 
striking observation that even the fluctuations of primordial parton density in 
individual, single events, are preserved in the single event variation of elliptic 
flow [267[ I268J . Moreover, the observed scaling of elliptic flow with hadron va- 
lence quark number [9l] - and thus not with hadronic mass as in radial flow |103j 
- conflrms the implied origin of the elliptic flow signal from the partonic phase. 

At the LHC the phase of early QCD plasma formation is expected to shift to 
yet higher energy density, in the vicinity of T = 600MeV and fj,b = 5MeV.0iie 
is thus getting nearer to the domain of QCD asymptotic freedom, and might ex- 
pect a falloff of the partonic cross section which is extremely high at RHIC [61] , 
as reflected by the small t]/s and A values. 

The observed features of the QCD plasma produced at RHIC energy have in- 
vited the terminology of a "strongly coupled" quark-gluon plasma (sQGP [l_65j). 
Further evidence for the strongly coupled, non perturbative nature of the pri- 
mordial partonic state stems from the various observed, strong in-medium at- 
tenuation effects on initially produced high pj- partons. In Au-|-Au collisions at 
^/s = 200GeV this high medium opacity leads to a universal quenching of the 
high pt hadron yield |175j including, most remarkably, heavy charm quark prop- 
agation to D mesons [175', '178] . We have shown in chapter [5] that the interior 
of the collisional flreball at midrapidity is almost "black" at t < 1/m/c.This is 
also reflected in a strong suppression of the back-to-back correlation of hadrons 
from primordially produced di-jets [207^1208] . and in a similarly strong suppres- 
sion of the J/^ yield [237] which we have shown in chapter [6] to be ascribed to 
an in-medium dissolution of the primordially produced cc pairs |41j at T about 
SOOMeV. 



The underlying high medium opacity can be formally expressed [1881 11931 I195J 
by an effective parton transport coefficient q (equations 15.81 and I5.12p which 
quantifies the medium induced transverse momentum squared per unit mean 
free path A. The value of q derived from analysis of the various attenuation 
phenomena turns out to be almost an order of magnitude higher than what 
was expected from former, perturbative QCD based models [196] . Analogously, 
rj/s has turned out to be much smaller than the previous perturbative QCD 
expectation |92]. The two quantities may be related |300j via the heuristic 
expression 

-^ ^ 3.75 C — (8.1) 

s q 

with C a to be determined constant; C = 1/3 from |300j . This relation shows 
that a larger value of q implies a small value for the ratio rj/s. The latter has a 
lower bound by the general quantum gauge field theory limit rj/s > {Att)~^ [171j . 
a value not too far from the estimate rj/s = 0.09 it 0.02 derived from relativistic 
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hydrodynamics applied to elliptic flow V2 [9HI170J . As a consequence, q can not 
grow beyond a certain upper bound that should be established at LHC energy. 

These considerations are an example of the recent intense theoretical search 
for alternative methods of real-time strong coupling calculations, complemen- 
tary to lattice QCD. In this regime, lattice QCD has to date been the prime 
non-perturbative calculational tool. However, understanding collective flow, 
jet quenching and primordial photon radiation requires real time dynamics, on 
which lattice QCD information is at present both scarce and indirect. Com- 
plementary methods for real-time strong coupling calculations at high temper- 
ature are therefore being explored. For a class of non-abelian thermal gauge 
field theories, the conjecture of a correspondence between anti-de Sitter space- 
time theory and conformal field theory (the so-called AdS/CFT conjecture) 
has been shown |301j to present such an alternative. It maps nonperturba- 
tive problems at strong coupling onto calculable problems of classical gravity 
in a five-dimensional anti-de Sitter (ADS5) black hole space-time theory |302j . 
In fact, this formalism has been recently applied [212J in a calculation of the 
transport coefficient q that governs in-medium jet attenuation, resulting in an 
effective, expansion time averaged q^ff = 5 GeV^/fm at T = 300 MeV corre- 
sponding to top RHIC energy, rather close to the experimental estimates (c.f. 
Figs. O and EHD- 

Thus, it does not seem to be too far-fetched to imagine |301j that the quark- 
gluon plasma of QCD, as explored at RHIC, and soon at the LHC (and the- 
oretically in lattice QCD), and the thermal plasma of certain supersymmetric 
conformal gauge field theories (for example N = 4 " Super- Yang-Mills" (SYM) 
theory as employed in |212tl30l] ) share certain fundamental common properties. 

The focus on early time in the dynamical evolution of matter in nucleus-nucleus 
collisions is specific to RHIC energy as the initial interpenetration period of 
two Lorentz contracted mass 200 nuclei amounts to 0.15fm/c only. The sub- 
sequent evolution is thus reflected in a multitude of observables. It is, at first, 
described as a one-dimensional Hubble expansion [[6lJ , setting the stage for the 
emergence of the medium specific quantities addressed above ( gluon satura- 
tion, direct photon production, hydrodynamic elliptic flow, jet quenching and 
J/^ suppression). These observables tend to settle toward their eventually ob- 
served patterns at t < 1.0 — 1.5fm/c, owing to the fact that they are most 
prominently determined under primordial conditions of high temperature and 
density. For example, photon production reflects T^, and the transport coeffi- 
cient q falls with T^ [212j . 

On the contrary, at the energy of former studies at the CERN SPS, 6 < ^/s < 
20GeV, such early times stay essentially unresolved as the initial interpenetra- 
tion takes upward of 1.5fm/c. A fireball system in local (or global) equilibrium 
thus develops toward t = 3fm/c, at T about 220MeV, closer to the onset of 
hadronization |85t I86j. Also the baryo-chemical potential is much higher than 
in RHIC collisions, 250 < fiB ^ 4:50MeV . However, we thus gain insight into 
the QCD physics of the hadronization, and high fis domain of the phase di- 
agram sketched in Fig. II. H in the vicinity of the conjectured parton-hadron 
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coexistence line of QCD. 

For reference of such data, e.g. statistical species equilibrium (chapter [3]) , dilep- 
ton studies of light vector meson "melting" at the phase boundary (chapter l6.3|) . 
and hadronic event-by-event fluctuations (chapters 17.21 and I7.3p , to theoretical 
QCD predictions, a recent progress of lattice QCD [Q [DJ |TU] is of foremost 
importance. The technical limitation of lattice QCD to the singular case of 
vanishing chemical potential, hb = (which arises from the Fermion determi- 
nant in the lattice formulation of the grand canonical partition function), has 
been overcome recently. Three different approaches have been outlined, the re- 
spective inherent approximation schemes touching upon the limits of both the 
mathematical and conceptual framework of lattice theory, and of present day 
computation power even with multi-teraflop machines. First results comprise 
the prediction of the parton-hadron phase boundary line, which interpolates 
between the well studied limits of the crossover-region at /ig — > 0, T > Tc and 
the high net baryon density, low T region for which a first order character of 
the phase transition has been predicted by chiral QCD models [15J. We have 
illustrated this line in Fig. 11.11 and we have shown in chapter [3] that hadronic 
freeze-out occurs at, or near this line at ^/s > 17.3GeV (top SPS energy). The 
coexistence line includes an intermediate (T, fis) domain featuring a critical 
point of QCD at which the first order line at higher fis terminates, in a critical 
domain of (T, ^b) in which the transition is of second order. One thus expects 
the nature of the confining hadronization transition — an open QCD question 
— to change from a crossover to a second order, and onward to a first order 
characteristics in a relatively small intervall of fiB that is accessible to nuclear 
collision dynamics at the relatively modest ^/s of about 5 to 15 GeV. This do- 
main has as of yet only received a first experimental glance, but the top AGS 
and low SPS energy experiments exhibit phenomena that can be connected to 
the occurrence of a critical point and/or a first order phase transition, notably 
the "SPS plateau" in (rn-j'), the non-monotoneous K~^ /tt~^ excitation function, 
and the eventwise fluctuations of this ratio (chapters 12.61 El 17.21 and 17. 3p . A 
renewed effort is underway at RHIC, at the CERN SPS and at the future GSI 
FAIR facility to study hadronization, in-medium meson modification induced 
by the onset of QCD chiral restoration, as well as critical fluctuations and flow, 
in the low ^/s domain. 
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